BZOJ2656 [Zjoi2012]数列(sequence)[模拟]

这个递推式子可以发现$i$是偶数下标可以缩一半，是奇数下标就可以拆成两个下标，$\lfloor \frac{i}{2} \rfloor$以及$\lfloor \frac{i}{2}+1 \rfloor$。

然后发现每次奇数下标可以拆成一奇一偶，而拆出来的奇再拆和偶的减半得到的下标只有两种。比如：

$A_{13}=1*A_6+1*A_7$

　　 $=1*A_3+1*A_3+1*A_4$

　　 $=2*A_3+1*A_4$

　　 $=2*A_1+2*A_2+1*A_2$

　　 $=2*A_1+3*A_2$

大概就是这个亚子。于是只要记录每次拆出来的两个数$a,b$各需要多少个加起来才能表示原数，

再拆时讨论一下是$b$奇$a$偶，还是$a$奇$b$偶，分两种情况将这个系数累加即可。一直折半到$A_1$停下输出。

---

高精又WA了一发。。我好菜。WA:line31是b.n。

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #include<cmath>
 #define dbg(x) cerr << #x << " = " << x <<endl
 using namespace std;
 typedef long long ll;
 typedef double db;
 typedef pair<int,int> pii;
 template<typename T>inline T _min(T A,T B){return A<B?A:B;}
 template<typename T>inline T _max(T A,T B){return A>B?A:B;}
 template<typename T>inline char MIN(T&A,T B){return A>B?(A=B,):;}
 template<typename T>inline char MAX(T&A,T B){return A<B?(A=B,):;}
 template<typename T>inline void _swap(T&A,T&B){A^=B^=A^=B;}
 template<typename T>inline T read(T&x){
     x=;int f=;char c;while(!isdigit(c=getchar()))if(c=='-')f=;
     while(isdigit(c))x=x*+(c&),c=getchar();return f?x=-x:x;
 }
 const int N=+,base=1e8,dig=;
 struct thx{
     int a[N],n;
     inline void mst(){memset(a,,sizeof a),n=;}
     thx(){mst();}
     inline int &operator [](int x){return a[x];}
     inline void get_num(int x){mst(),n=;while(x)a[++n]=x%base,x/=base;}
     inline void operator =(int x){get_num(x);}
     inline thx operator +(thx&orz){
         thx b;b.n=_max(n,orz.n);
         for(register int i=;i<=b.n;++i)b[i]+=a[i]+orz[i],b[i+]=b[i]/base,b[i]%=base;
         while(b[b.n+]&&++b.n);//mistake!
         return b;
     }
     inline void operator +=(thx&orz){*this=*this+orz;}
     inline void operator ++(){
         ++a[];int tmp=;
         while(a[tmp]==base)a[tmp]=,++a[++tmp];
         while(a[n+]&&++n);
     }
     inline void div2(){
         for(register int i=n;i;--i)a[i-]+=(a[i]&)*base,a[i]/=;
         while(!a[n]&&--n);if(!n)n=;
     }
     inline void read(){
         mst();
         char s[+];scanf("%s",s+);int len=strlen(s+);
         for(register int i=;dig*(i-)<len;++i)
             for(register int j=,bas=;j<=_min(dig,len-dig*(i-));++j,bas*=)a[i]+=(s[len-dig*(i-)-j+]-'')*bas;
         n=len/dig+(len%dig!=);
     }
     inline void print(){
         printf("%d",a[n]);
         for(register int i=n-;i;--i)printf("%08d",a[i]);
         puts("");
     }
 }val,a,b,x,y,tmp;
 int T;

 int main(){//freopen("9.in","r",stdin);freopen("9.ans","w",stdout);
     read(T);while(T--){
         val.read();
         if(val.n==&&val[]==){puts("");continue;}
         while(!(val[]&))val.div2();
         if(val.n==&&val[]==){puts("");continue;}
         val.div2();b=a=val;++b;x=,y=;
         while(!(a.n==&&a[]==)){//printf("x="),x.print(),printf("y="),y.print();
             if(a[]&)y+=x,a.div2(),b=a,++b;
             else x+=y,a.div2(),b=a,++b;
         }//puts("done");
         x+=y;x.print();
     }
     return ;
 }