HDOJ1568 Fibonacci[公式求前四位数]

Fibonacci

Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 2113    Accepted Submission(s): 979

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

 

Author

daringQQ

 

Source

Happy 2007

 

Recommend

8600

 

 

 

 

 

标准的Fibonacci数列F(0) = 0, F(1) = 1, F(n) = F(n-1) + F(n - 2)。通项公式为F(n) = (((1+Sqrt(5))/2)^n - ((1-Sqrt(5))/2)^n)*1/Sqrt(5)，因为(1-sqrt(5))/2的绝对值小于1，所以当i较大的时候，往往可以忽略掉这一项F(n)≈((1+Sqrt(5))/2)^n/sqrt(5)。所以Fibonacci前30项可直接求出。后面的大数取log10，取整后求出左边4位，方法和HDOJ1060类似http://hi.baidu.com/racebug/blog/item/e3b5cdefeb485af2b2fb9537.html

这个是对Fibonacci介绍很全的一个论文http://blog.163.com/baobao_zhang@126/blog/static/482523672008826103832368/

转自：http://hi.baidu.com/racebug/blog/item/54ad34ed3c3dfc3426979191.html

code ：

#include<iostream>
#include<cmath>
using namespace std;
int main()
{
    int n;
    int i;
    int f[31];
    double temp1=log10((1.0+sqrt(5.0))/2.0),temp2=log10(sqrt(5.0));
    f[0]=0;
    f[1]=1;
    for(i=2;i<31;i++)
        f[i]=f[i-1]+f[i-2];
    while(~scanf("%d",&n))
    {
        if(n<31)
        {
            int temp=f[n];
            while(temp>=10000)
                temp/=10;
            printf("%d\n",temp);
        }
        else
        {
            double t;
            int ans;
            t=(n*temp1-temp2);
            t-=(int)t;
            ans=pow(10.0,t)*1000.0;
            printf("%d\n",ans);
        }
    }
    return 0;
}