首页 > 代码库 > 【BZOJ 2656】2656: [Zjoi2012]数列(sequence) (高精度)
【BZOJ 2656】2656: [Zjoi2012]数列(sequence) (高精度)
2656: [Zjoi2012]数列(sequence)
Time Limit: 2 Sec Memory Limit: 128 MB
Submit: 1499 Solved: 786Description
小白和小蓝在一起上数学课,下课后老师留了一道作业,求下面这个数列的通项公式:
小白作为一个数学爱好者,很快就计算出了这个数列的通项公式。于是,小白告诉小蓝自己已经做出来了,但为了防止小蓝抄作业,小白并不想把公式公布出来。于是小白为了向小蓝证明自己的确做出来了此题以达到其炫耀的目的,想出了一个绝妙的方法:即让小蓝说一个正整数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
Source
【分析】
就是可以直接高精度暴力的。因为每次最多分出两个数而已。
比如:
$25→12,13→6,7→3,4→1,2→0,1$
1 #include<cstdio> 2 #include<cstdlib> 3 #include<cstring> 4 #include<iostream> 5 #include<algorithm> 6 using namespace std; 7 8 int mymax(int x,int y) {return x>y?x:y;} 9 10 struct hugeint 11 { 12 int a[110],l; 13 void clear() 14 { 15 memset(a,0,sizeof(a)); 16 l=1; 17 } 18 friend hugeint operator + (hugeint x,hugeint y) 19 { 20 int ll=mymax(x.l,y.l); 21 for(int i=1;i<=ll;i++) x.a[i]=x.a[i]+y.a[i]; 22 for(int i=1;i<=ll;i++) x.a[i+1]+=x.a[i]/10,x.a[i]%=10; 23 while(x.a[ll+1]!=0) x.a[ll+2]+=x.a[ll+1]/10,x.a[ll+1]%=10,ll++; 24 x.l=ll; 25 return x; 26 } 27 friend hugeint operator + (hugeint x,int y) 28 { 29 int ll=x.l; x.a[1]+=y; 30 for(int i=1;i<=ll;i++) x.a[i+1]+=x.a[i]/10,x.a[i]%=10; 31 while(x.a[ll+1]!=0) x.a[ll+2]+=x.a[ll+1]/10,x.a[ll+1]%=10,ll++; 32 x.l=ll; 33 return x; 34 } 35 friend hugeint operator / (hugeint x,int y) 36 { 37 int nw=0; 38 for(int i=x.l;i>=1;i--) 39 { 40 nw=nw*10+x.a[i]; 41 x.a[i]=nw/y; 42 nw%=y; 43 } 44 while(x.a[x.l]==0&&x.l>1) x.l--; 45 return x; 46 } 47 }; 48 49 hugeint k1,k2; 50 void dfs(hugeint xx) 51 { 52 if(xx.l==1&&xx.a[1]==1) 53 { 54 k1=xx;k2.clear(); 55 return; 56 } 57 dfs((xx+1)/2); 58 if(xx.a[1]&1) k1=k1+k2; 59 else k2=k1+k2; 60 } 61 62 char s[110]; 63 64 int main() 65 { 66 int T; 67 scanf("%d",&T); 68 while(T--) 69 { 70 scanf("%s",s+1); 71 int l=strlen(s+1); 72 hugeint xx;xx.clear(); 73 for(int i=1;i<=l;i++) xx.a[l-i+1]=s[i]-‘0‘; 74 xx.l=l; 75 dfs(xx); 76 for(int i=k1.l;i>=1;i--) printf("%d",k1.a[i]); 77 printf("\n"); 78 } 79 return 0; 80 }
【BZOJ 2656】2656: [Zjoi2012]数列(sequence) (高精度)
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。