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 = 105;
const int base = 1e4;

struct Bign {
    int len,N[maxn];
    int& operator[](const int& i){
        return N[i];
    }
    void init() {
        len=1;
        memset(N,0,sizeof(N));
    }
};
Bign tmp; char s[maxn];
void read(Bign& A)
{
    A.init(),tmp.init();
    scanf("%s",s+1);
    int i,j;
    tmp.len=strlen(s+1);
    for(int i=1;i<=tmp.len;i++)
        tmp[i]=s[tmp.len+1-i]-'0';
    A.len=(tmp.len-1)/4 +1;
    for(i=1;i<=A.len;i++) {
        for(j=1;j<=4;j++)
            A[i]=A[i]*10+tmp[(i-1)*4+(5-j)];
    }
}
void print(Bign A)
{
    printf("%d",A[A.len]);
    for(int i=A.len-1;i;i--)
        printf("%04d",A[i]);
    puts("");
}
Bign operator + (Bign A,Bign B)
{
    A.len=max(B.len,A.len);
    for(int i=1;i<=A.len;i++) {
        A[i]+=B[i];
        A[i+1]+=A[i]/base;
        A[i]%=base;
    }
    while(A[A.len+1]) A.len++;
    return A;
}
Bign operator + (Bign A,int x)
{
    Bign B; B.init(); B[1]=x;
    return A+B;
}
Bign operator / (Bign A,int p)
{
    for(int i=1;i<=A.len;i++) {
        if(A[i]&1) A[i-1]+=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==1 && x[1]==1) {
        p=x; q.init(); return ;
    }
    F((x+1)/2);
    if(x[1]&1) 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 0;
}

 

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