bzoj 2656 [Zjoi2012]数列(sequence)（高精度）

【题目链接】

　　http://www.lydsy.com/JudgeOnline/problem.php?id=2656

【题意】

　　计算大数递推式

【思路】

　　高精度

【代码】

 #include<cstdio>
 #include<cstring>
 #include<iostream>
 using namespace std;

 const int maxn = ;
 const int base = 1e4;

 struct Bign {
     int len,N[maxn];
     int& operator[](const int& i){
         return N[i];
     }
     void init() {
         len=;
         memset(N,,sizeof(N));
     }
 };
 Bign tmp; char s[maxn];
 void read(Bign& A)
 {
     A.init(),tmp.init();
     scanf("%s",s+);
     int i,j;
     tmp.len=strlen(s+);
     for(int i=;i<=tmp.len;i++)
         tmp[i]=s[tmp.len+-i]-'';
     A.len=(tmp.len-)/ +;
     for(i=;i<=A.len;i++) {
         for(j=;j<=;j++)
             A[i]=A[i]*+tmp[(i-)*+(-j)];
     }
 }
 void print(Bign A)
 {
     printf("%d",A[A.len]);
     for(int i=A.len-;i;i--)
         printf("%04d",A[i]);
     puts("");
 }
 Bign operator + (Bign A,Bign B)
 {
     A.len=max(B.len,A.len);
     for(int i=;i<=A.len;i++) {
         A[i]+=B[i];
         A[i+]+=A[i]/base;
         A[i]%=base;
     }
     while(A[A.len+]) A.len++;
     return A;
 }
 Bign operator + (Bign A,int x)
 {
     Bign B; B.init(); B[]=x;
     return A+B;
 }
 Bign operator / (Bign A,int p)
 {
     for(int i=;i<=A.len;i++) {
         if(A[i]&) A[i-]+=base/p;
         A[i]/=p;
     }
     while(!A[A.len]&&A.len) A.len--;
     return A;
 }

 Bign p,q,n;

 void F(Bign x)
 {
     if(x.len== && x[]==) {
         p=x; q.init(); return ;
     }
     F((x+)/);
     if(x[]&) p=p+q;
     else q=p+q;
 }

 int main()
 {
     int T; scanf("%d",&T);
     while(T--) {
         p.init(),q.init();
         read(n);
         F(n);
         print(p);
     }
     return ;
 }

Ps：每次做高精度都会涨姿势 >_<