Fibonacci

  
  

   
   

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

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

   
   输出f
   
   
    
    
     
     n
    
    
     
     
      
      
        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

  
  


  
  


  
  


  
  
题解：
  
  
斐波那契通项　f[n]=(1/sqrt(5)) * { [(1 + sqrt(5)) / 2] ^ n - [(1 - sqrt(5)) / 2]^ n }

  
  
首先通过该通过该通项可以得出第n项b,由于b太大只能对其取对数：ans=lg(b)=lg(10^(k-3)*(x1*1000+x2*100+x3*10+x4+0.x5+........))
  
  
得：ans=k-3+lg(x1*1000+x2*100+x3*10+x4+0.x5+....)括号内取整即为所求的答案ret
  
  
则ret=int10^(ans-k+3)=int 10^(ans-int(ans)+3).
  
  


  
  
#include"stdio.h"
#include"cstdio"
#include"algorithm"
#include"string.h"
#include"math.h"
using namespace std;
int main()
{
	int A[50];
	A[0]=0;A[1]=1;
	for(int i=2;i<25;i++)
	A[i]=A[i-1]+A[i-2];
	int n;
	while(scanf("%d",&n)!=EOF)
	{
		int ret;
		if(n<=20)
		ret=A[n];
		else
		{
			double ans=log10(1.0/sqrt(5))+(double)n*log10((1+sqrt(5))/2.0);
			ret=pow(10.0,(ans-int(ans)+3));
		}
		printf("%d\n",ret);
	}
}