HDacm212

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说的是否正确。

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

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

公式为：F(n)=（√5/5)*{[(1+√5)/2]^n - [(1-√5)/2]^n}

 

n比较小时n<=30，用f(n)=f(n-1)+f(n-2),n>=2求解

n比较大时使用进似计算的通式f(n)=x^n  ,x=(1+sqrt(5))/2  ,  (n<10^10)

巧妙之处在于计算x^n   

 

#include <stdio.h>
#include <math.h>
int main()
{
    int i,f[31]={0,1},n;
    double n1,n2,t,s;
    for(i=2;i<=30;i++)
    f[i]=f[i-1]+f[i-2];
    for(i=2;i<=30;i++)
    while(f[i]>=10000)
    f[i]=f[i]/10;
    while(scanf("%d",&n)!=EOF)
    {
        if(n<=30)
        {
            printf("%d\n",f[n]);
        }
        else
        {
            n1=log10((double)sqrt(5.0));
            n2=n*log10((double)sqrt(5.0)/2+0.5);
            t=n2-n1;
            s=t-(int)t;
            s=pow(10.0,s)*1000;
            printf("%d\n",(int)s);
 
        }
 
    }
    return 0;
 
}