HDU 3117 Fibonacci Numbers

 

http://acm.hdu.edu.cn/showproblem.php?pid=3117

 

题目大意:给定一个K（K很大很大），请你输出第K个Fibonacci数的前四位和后四位。

 

解题思路：这道题一开始后四位我直接就是想着去打表，我只是猜测说后四位的数字可能会形成一个周期，所以我就写了个测试代码，来测试是否存在周期，结果发现果然存在周期15000，所以就打表打出了后四位的值。前四位的求法也是参考了别人的思路后才了解的。

 

我们要搬出Fibonacci的通项公式——
f(n)=1/sqrt(5)(((1+sqrt(5))/2)^n+((1-sqrt(5))/2)^n)

假设F[n]可以表示成 t * 10^k（t是一个小数），那么对于F[n]取对数log10，答案就为log10 t + K，此时很明显log10 t<1，于是我们去除整数部分，就得到了log10 t ，再用pow（10，log10 t)我们就还原回了t。将t×1000就得到了F[n]的前四位。

具体实现的时候Log10 F[n]约等于((1+sqrt(5))/2)^n/sqrt(5)，这里我们把((1-sqrt(5))/2)^n这一项忽略了，因为当N>=40时，这个数已经小的可以忽略。于是log10 F[n]就可以化简成log10 1/sqrt(5) + n*log10 (1+sqrt(5))/2

 

#include <stdio.h> #include <math.h> #define aa (sqrt(5.0)+1.0)/2 int last[15000]={0,1,1,2,3,5}; int main() { int i,n,Fibonacci[40]={0,1,1,2,3,5}; double ans; /* 测试周期 int T = 0; for (i=3;i<=120000;i++) { last[i] = (last[i-2]+last[i-1])%10000; if(last[i] == last[2]&&last[i-1] == last[1]) { T = i; break; } } printf("%d/n",T-2); */ /*后4位打表*/ for (i=3;i<15000;i++) last[i] = (last[i-2]+last[i-1])%10000; for(i=6;i<40;i++) Fibonacci[i] = Fibonacci[i-1]+Fibonacci[i-2]; while (scanf("%d",&n)!=EOF) { if(n<40) printf("%d/n",Fibonacci[n]); else { ans=-0.5*(log10(5.0))+n*log10(aa); ans-=(int)ans; ans=pow(10.0,ans); while(ans<1000) ans*=10; printf("%d...",(int)ans); printf("%4.4d/n",last[n%15000]); } } return 0; }