bzoj 2656: [Zjoi2012]数列(sequence)（简单高精度模板2.0）

2656: [Zjoi2012]数列(sequence)

Time Limit: 2 Sec Memory Limit: 128 MB
Submit: 1663 Solved: 860
[Submit][Status][Discuss]

Description

小白和小蓝在一起上数学课，下课后老师留了一道作业，求下面这个数列的通项公式：

小白作为一个数学爱好者，很快就计算出了这个数列的通项公式。于是，小白告诉小蓝自己已经做出来了，但为了防止小蓝抄作业，小白并不想把公式公布出来。于是小白为了向小蓝证明自己的确做出来了此题以达到其炫耀的目的，想出了一个绝妙的方法：即让小蓝说一个正整数N，小白则说出的值，如果当N很大时小白仍能很快的说出正确答案，这就说明小白的确得到了公式。但这个方法有一个很大的漏洞：小蓝自己不会做，没法验证小白的答案是否正确。作为小蓝的好友，你能帮帮小蓝吗？

Input

输入文件第一行有且只有一个正整数T，表示测试数据的组数。

第2～T+1行，每行一个非负整数N。

Output

输出文件共包含T行。

第i行应包含一个不含多余前缀0的数，它的值应等于An(n为输入数据中第i+1行被读入的整数)

【样例输入】

Sample Input

Sample Output

假设n是奇数

那么有F(n) = F(n/2)+F(n/2+1) = 2F(n/4)+F(n/4+1) = …

可以logn的复杂度递归

n是大数，需要高精度

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
char str[105];
typedef struct Bnum
{int l;int s[108];Bnum(){l = 0;memset(s, 0, sizeof(s));}
}Bnum;
Bnum operator + (Bnum x, Bnum y)
{int i;Bnum z;z.l = max(x.l, y.l);for(i=1;i<=z.l;i++){z.s[i] += x.s[i]+y.s[i];if(z.s[i]>9){z.s[i+1] += 1;z.s[i] %= 10;}}while(z.s[z.l+1])z.l++;return z;
}
Bnum operator - (Bnum x, Bnum y)
{int i;Bnum z;z.l = max(x.l, y.l);for(i=1;i<=z.l;i++){if(x.s[i]<y.s[i]){x.s[i] += 10;x.s[i+1]--;}z.s[i] += x.s[i]-y.s[i];if(z.s[i]>9){z.s[i+1] += z.s[i]/10;z.s[i] %= 10;}}while(z.s[z.l]==0)z.l--;return z;
}
Bnum operator * (Bnum x, Bnum y)
{int i, j;Bnum z;z.l = x.l+y.l-1;for(i=1;i<=x.l;i++){for(j=1;j<=y.l;j++){z.s[i+j-1] += x.s[i]*y.s[j];if(z.s[i+j-1]>9){z.s[i+j] += z.s[i+j-1]/10;z.s[i+j-1] %= 10;}}}while(z.s[z.l+1])z.l++;return z;
}
Bnum operator ^ (Bnum x, int y)
{Bnum z, u;z.l = z.s[1] = 1, u = x;while(y){if(y%2==1)z = z*u;u = u*u;y /= 2;}return z;
}
Bnum Half(Bnum x)
{int i, now;Bnum z;z.l = x.l, now = 0;for(i=z.l;i>=1;i--){z.s[i] = (now*10+x.s[i])/2;if(x.s[i]%2)now = 1;elsenow = 0;}while(z.s[z.l]==0)z.l--;return z;
}
void Print(Bnum x)
{int i;for(i=x.l;i>=1;i--)printf("%d", x.s[i]);printf("\n");
}
int main(void)
{Bnum x, y, e, a, b;int T, i;scanf("%d", &T);while(T--){scanf("%s", str+1);x.l = strlen(str+1);if(x.l==1 && str[1]=='0'){printf("0\n");continue;}for(i=x.l;i>=1;i--)x.s[i] = str[x.l-i+1]-'0';while(x.s[1]%2==0)x = Half(x);e.l = e.s[1] = 1;x = Half(x), y = x+e;a.l = a.s[1] = b.l = b.s[1] = 1;while(1){//Print(a); Print(x);//Print(b); Print(y);if(x.l==1 && x.s[1]==1){if(y.l==1 && y.s[1]==2)Print(a+b);elsePrint(a);break;}if(y.l==1 && y.s[1]==1){if(x.l==1 && x.s[1]==2)Print(a+b);elsePrint(b);break;}if(x.s[1]%2==0){a = a+b;x = Half(x);y = x+e;}else{b = a+b;y = Half(y);x = y-e;}}}return 0;
}