NYOJ Fibonacci数列（四）

Fibonacci数列有通项公式为：
化简log10(f_n)=n*log10((1+sqrt(5))/2)-log10(sqrt(5))+log10(1-((1-sqrt(5))/(1+sqrt(5)))^n)后面红色部分是无穷小量，可以省略。
于是log10(f_n)=n*log10((1+sqrt(5))/2)-log10(sqrt(5))；

最后可以计算出log10(f_n)，保留它的的小数部分即可。

因为之前对fn求log10的对数，整数部分相当于(fn的位数-1)，把整数部分减去后，再求的10^t即为原来的fn除以10^(fn的位数-1)【比如：如果某一项为123456，则10^t=1.23456=123456/100000】，然后乘以1000之后取整即为所求

#include<cstdio>
#include<cmath>
#include<iostream>
using namespace std;
int f[23]={0,1};
double x=log((sqrt(5.0)+1.0)/2.0)/log(10.0);
double y=-0.5*log(5.0)/log(10.0);
int main(){
    int n,i;
    for(i=2;i<21;++i)
        f[i]=f[i-1]+f[i-2];//算出前20项，小于10000的
    while(scanf("%d",&n)!=EOF){
        if(n<=20) printf("%d\n",f[n]);//前二十个直接输出
        else {
            double t=n*x+y;//算出log10(fn)
            t-=floor(t);//减去整数部分
            t=pow(10.0,t);//因为之前对fn求log10的对数，整数部分相当于fn的位数-1，把整数部分减去后，再求的10^t即为原来的fn除以10^(fn的位数-1)【比如：如果某一项为123456，则10^t=1.23456=123456/100000】，然后乘以1000之后取整即为所求
            t*=1000;
            printf("%d\n",(int)t);
        }
    }
}