Fibonacci ：斯特林公式

#include<iostream>
#include<cmath>
using namespace std;
int a[21]={0,1};
int main()
{
    int n;
    double m,l,des;
    while(~scanf("%d",&n))
    {
        for(int i=2;i<=20;i++)
        {
            a[i]=a[i-1]+a[i-2];
        }
        if(n<=20)printf("%d\n",a[n]);
        else 
        {
            m=(1.0+sqrt(5.0))/2.0;
            l=-log10(5.0)*0.5+n*1.0*log10(m);
            l-=(int)l;
            des=pow(10.0,l);
            while(des<1000)
            {
                des*=10;
            }
            printf("%d\n",(int)des);
        }  
    }
    return 0;
}