【递推】【高精度】[ZJOI2012]数列 sequence

[Zjoi2012]数列(sequence)

Time Limit: 2 Sec  Memory Limit: 128 MB
Submit: 421  Solved: 216
[ Submit][ Status][ Discuss]

Description

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

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

Input

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

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

Output

      输出文件共包含T行。

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

【样例输入】

Sample Input

3

1

3

10

Sample Output

1
2
3

HINT

T<=20,N<=10^100

题目不能再简单了，明显的高精度

递归写完，发现超时。。。。。。

尝试记忆化，数据太大存不下。。。。。。

只能另想办法

比如

   9
4    5
2  2   3
1  1  1  2
         1

有很多重复调用，但是我们发现，算完1，可以算出2,3，然后可以算出4,5，然后就是9,10

利用这一特性，我们设 p 为 x 的答案，然后 q 为 (x-1) 的答案，y =2*x ，即 x 和 (x-1) 的父亲

每次只需要在p和q之间转换就可以得出答案

测评情况（Cena）

C++ AC Code

/*http://blog.csdn.net/jiangzh7
By Jiangzh*/
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

struct BigInt{//压位5为
	int len,c[50];
	BigInt(){len=1;memset(c,0,sizeof(c));}
	void Clear(){len=1;memset(c,0,sizeof(c));}
	void Zero(){while(c[len]==0 && len>1) len--;}
	void Read()
	{
		memset(c,0,sizeof(c));len=1;
		char s[200];
		memset(s,0,sizeof(s)); fgets(s,190,stdin);
		int slen=strlen(s);
		if(s[slen-1]=='\n') {s[slen]='\0';slen--;}
		int k=1;
		for(int i=slen-1;i>=0;i--)
		{
			c[len]+=(s[i]-'0')*k;
			k*=10;
			if(k==100000) {len++;k=1;}
		}
		Zero();
	}
	void Print()
	{
		printf("%d",c[len]);
		for(int i=len-1;i>=1;i--) printf("%05d",c[i]);
		puts("");
	}
	inline bool operator == (int a)
	{
		if(len==1 && c[1]==a) return true;
		return false;
	}
	inline int operator & (int a)
	{
		return c[1]&a;
	}
	BigInt operator + (int a)
	{
		BigInt res;res.len=len+1;
		for(int i=1;i<=len;i++)
		{
			res.c[i]+=c[i]+a;a=0;
			res.c[i+1]=res.c[i]/100000;
			res.c[i]%=100000;
		}
		res.Zero();
		return res;
	}
	BigInt operator + (BigInt a)
	{
		BigInt res;res.len=max(len,a.len)+1;
		for(int i=1;i<=res.len;i++)
		{
			res.c[i]+=c[i]+a.c[i];
			res.c[i+1]=res.c[i]/100000;
			res.c[i]%=100000;
		}
		res.Zero();
		return res;
	}
	BigInt operator / (int a)
	{
		BigInt res;res.len=len;
		for(int i=res.len;i>=1;i--)
		{
			res.c[i]=(res.c[i]+c[i])/a;
			if(i>1) res.c[i-1]+=c[i]%a*100000;
		}
		res.Zero();
		return res;
	}
};

BigInt n;
BigInt p,q;


void calc(BigInt x)
{
	//printf("---> ");x.Print();
	if(x==1)
	{
		p=x; q.Clear();
		return;
	}
	calc((x+1)/2);
	if(x&1) p=p+q;
	else q=q+p;
}

int main()
{
	freopen("sequence.in","r",stdin);
	freopen("sequence.out","w",stdout);
	int grp;scanf("%d\n",&grp);
	while(grp--)
	{
		n.Read();
		calc(n);
		p.Print();
	}
	return 0;
}