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说的是否正确。
(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
斐波那契数列公式取对数最后一项 log10(1-((1-√5)/(1+√5))^n)的值趋近于0,可忽略
所以log10(f[n])=log10(1/sqrt(5))+n*log10((1+sqrt(5))/2=n*log10((1+sqrt(5))/2)-log10(sqrt(5));
令log10(f(n))=m; m为小数 10^m=f(n);
小数部分为m=m-(int)m;
假设给出一个数10234432,
那么log10(10234432)=log10(1.0234432*10^7)=log10(1.0234432)+7;
log10(1.0234432)就是log10(10234432)的小数部分.
log10(1.0234432)=0.010063744
10^0.010063744=1.023443198
#include<stdio.h> #include<math.h> int main() { int n,a[21],i; double m; a[0]=0; a[1]=1; for(i=2;i<=20;i++) a[i]=a[i-1]+a[i-2]; while(scanf("%d",&n)!=EOF) { if(n<21) printf("%d\n",a[n]); else { m=n*log10((1+sqrt(5))/2)-log10(sqrt(5)); //公式 m=m-(int)m; n=pow(10,m)*1000; printf("%d\n",n); } } return 0; }