#include <iostream>
#include <cstring>
#include <cmath>
#include <cstdio>
using namespace std;
int f[101];
int main()
{
int n,i,j,k;
f[0]=0;f[1]=1;
for(i=2;i<=20;i++)
f[i]=f[i-1]+f[i-2];
while(cin>>n)
{
if(n<=20){cout<<f[n]<<endl;continue;}
double ans,p;
int s;
p=(1+sqrt(5.0))/2;
ans=-0.5*log10(5.0)+n*log10(p);
ans=ans-(int)ans;
s=(int)(pow(10.0,ans)*1000);
printf("%d\n",s);
}
return 0;
}
/*
f[n]=f[n-1]+f[n-2];
在n<=20,即f[n]<10000时,用数组保存f[n],输入判断n<=20则直接用输出。
当n>20时,需要用通项公式求;
斐波那契通项公式:f[n]=1/sqrt(5)*(((1+sqrt(5))/2)^n-((1-sqrt(5))/2)^n)
变化: =1/sqrt(5)*((1+sqrt(5))/2)^n*(1-((1-sqrt(5))/(1+sqrt(5)))^n);
令p=(1+sqrt(5))/2;
求对数:log(f[n])=-0.5*log(5.0)+n*log(p)+log(1-((1-sqrt(5))/(1+sqrt(5)))^n);
其中log(1-((1-sqrt(5))/(1+sqrt(5)))^n)趋于0,当n>20时
所以log(f[n])=-0.5*log(5.0)+n*log(p);
其整数部分为a,和小数部分b。
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)
*/
类似的有hdu2855 ,这题除了需要求前4位,还要求后4位。题解代码链接:http://blog.csdn.net/a601025382s/article/details/10195939