HDOJ 题目1568 Fibonacci（数学，公式）

Fibonacci

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

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   |   We have carefully selected several similar problems for you:   1788  1787  1573  3641  3123 

 思路：

1. f[n]=f[n-1]+f[n-2]; 
2.     在n<=20,即f[n]<10000时，用数组保存f[n],输入判断n<=20则直接用输出。 
3.     当n>20时，需要用通项公式求； 
4.      
5.     斐波那契通项公式：f[n]=1/sqrt(5)*(((1+sqrt(5))/2)^n-((1-sqrt(5))/2)^n) 
6.                     变化： =1/sqrt(5)*((1+sqrt(5))/2)^n*(1-((1-sqrt(5))/(1+sqrt(5)))^n); 
7.     令p=(1+sqrt(5))/2； 
8.     求对数:log(f[n])=-0.5*log(5.0)+n*log(p)+log(1-((1-sqrt(5))/(1+sqrt(5)))^n); 
9.     其中log(1-((1-sqrt(5))/(1+sqrt(5)))^n)趋于0，当n>20时 
10.      
11.     所以log(f[n])=-0.5*log(5.0)+n*log(p); 
12.     其整数部分为a，和小数部分b。 
13.     a表示f[n]的长度，则(int)(10^b*1000)就表示f[n]的前4位。(f[n]=10^(a+b),则10^b=f[n]/10^a.例:f[n]=12345,10^b=1.2345)
14. ac代码

15. #include<stdio.h>
    #include<string.h>
    #include<math.h>
    int a[25];
    int main()
    {
    	int n,i;
    	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)
    	{
    		double p,ans;
    		if(n<=20)
    		{
    			printf("%d\n",a[n]);
    			continue;
    		}
    		p=(1+sqrt(5.0))/2;
    		ans=-0.5*log10(5.0)+n*log10(p);
    		ans-=(int)ans;
    		printf("%d\n",(int)(pow(10.0,ans)*1000));
    	}
    }