HDU-1568 Fibonacci【数学公式题】

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说的是否正确。

解题报告

由于数据范围很大，所以应该不能用递推法，可以考虑利用通项公式。

我们要求的是某一项的前四位数字，这就是要去一个大实数末尾的数字，可以用对数，由于$lg(n)$的整数部分就是这个数字的位数，$10^x$(x为小数部分)就是一个$(1,10)$的实数，可以看成第n项的科学计数法。

这样，我们只要快速算出第n项的$lg$值就可以了，利用对数的运算法则，很快可以得到

$$lg(f_n)=-lg5+nlg(\frac{1+\sqrt5}{1+\sqrt5})+lg[1-\big(\frac{1-\sqrt5}{1+\sqrt5}\big)^5]$$
这样，问题就解决了。

代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define g5 sqrt(5)
#define LL long long
#define lg(x) log10(x)
using namespace std;
int n,f[21];
double power(double x,int y){
    if(!y)return 1;
    if(y==1)return x;
    double c=power(x,y>>1);
    if(y&1)return x*c*c;
      else return c*c;
}
int main(){
    freopen("calc.in","r",stdin);
    freopen("calc.out","w",stdout);
    f[1]=1;
    for(int i=2;i<=20;i++)f[i]=f[i-1]+f[i-2];
    while(scanf("%d",&n)==1){
        if(n<=20){
            printf("%d\n",f[n]);
            continue;
        }
        double t=lg((double)1/g5)+(double)n*lg((1+g5)/2)+lg(1-power((1-g5)/(1+g5),n));LL nn=t;
        if(nn>3)t-=nn;
        t=pow(10,t);int ans=t*10000;ans/=10;
        printf("%d\n",ans);
    }
    return 0;
}