ACM刷题之HDU————Fibonacci

Fibonacci

Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)

Total Submission(s): 3846 Accepted Submission(s): 1747

Problem Description 2007年到来了。经过2006年一年的修炼，数学神童zouyu终于把0到100000000的Fibonacci数列 (f[0]=0,f[1]=1;f[i] = f[i-1]+f[i-2](i>=2))的值全部给背了下来。 接下来，CodeStar决定要考考他，于是每问他一个数字，他就要把答案说出来，不过有的数字太长了。所以规定超过4位的只要说出前4位就可以了，可是CodeStar自己又记不住。于是他决定编写一个程序来测验zouyu说的是否正确。

Input 输入若干数字n(0 <= n <= 100000000)，每个数字一行。读到文件尾。

Output             输出f[n]的前4个数字（若不足4个数字，就全部输出）。

Sample Input 0 1 2 3 4 5 35 36 37 38 39 40

Sample Output 0 1 1 2 3 5 9227 1493 2415 3908 6324 1023   这题用平常的方法，肯定会爆数据或者超时。一开始我想着取余，等数字大起来的时候舍弃后面的数字，但是会wa的。 百度了下，发现原来fibonacci数列还有一个通项公式的。 （图片来自百度百科） 但是如果直接这么算的话，还是会爆数据的。 所以我们就要用到log10了 下面是ac代码 #include<stdio.h> #include<math.h> #include<stdlib.h> int main() { double x,n,bit,b,c,d; int y; int q[100]={0,1,1,2,3,5,8,13,21,34,55,89,144,233,377,610,987,1597,2584,4181,6765}; while(scanf("%lf",&n)!=EOF) { if(n>20) { b=(1+sqrt(5.0))/2.0; c=(1-sqrt(5.0))/(b*2.0); x=n*log10(b)-log10(sqrt(5.0)); y=x; bit=x-(double)y; bit=pow(10,bit); while(bit<1000) { bit*=10; } y=bit; printf("%d\n",y); } else { y=(int)n; printf("%d\n",q[y]); } } }