数论（fabonacci数列） hdu-1568-Fibonacci

题目链接：

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

题目意思：

给一个i，让你求出第i个febonacci数的前四位，不足四位的直接输出。

解题思路：

由feibonacci数列的通项公式an=1/√5*(((1+√5)/2)^n+((√5-1)/2)^n)

当n比较大的时候可以舍掉（√5-1)/2)^n 因为越来越小，不到一，可以舍去

log10(an)=log10(1/√5)+n*log10((1+√5)/2) =p;

令p=p1+p2(其中p1为p的整数部分，p2为p的小数部分),则10^(p1+p2)=10^p1*10^p2=an 其中把10^p1为an的尾部的零，去掉后不影响高位的计算。

详细处理见代码。

代码：

#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<string>
#include<cstring>
#include<algorithm>
#include<vector>
#include<map>
#include<stack>
#include<list>
#include<queue>
#define eps 1e-6
#define INF (1<<30)
#define PI acos(-1.0)
using namespace std;
#define ba sqrt(5.0)

/*
freopen("data.in","r",stdin);
freopen("data.out","w",stdout);
*/

//fabonacci 通项公式an=1/√5*(((1+√5)/2)^n+((√5-1)/2)^n)
//当n比较大的时候可以舍掉（√5-1)/2)^n 因为越来越小，不到一

int save[30];
int main()
{
save[0]=0;
save[1]=1;
for(int i=2;i<20;i++)
{
save[i]=save[i-1]+save[i-2];
// printf("%d\n",save[i]);
} //前二十位都是小于5位的，直接算出来就行了
int n;

while(scanf("%d",&n)!=EOF)
{
if(n<20)
{
printf("%d\n",save
);
continue;
}
double temp=log10(1.0/ba)+n*log10((ba+1.0)/2.0); //以十为底求出对数

temp-=floor(temp);//除掉整数部分，也就是save
的尾部的零
temp=pow(10.0,temp); //算出高位非零的部分

while(temp<1000) //不足四位补零凑齐四位
temp*=10;
while(temp>=10000) //超过五位的除掉低位
temp/=10;
printf("%d\n",(int)temp);

}

return 0;
}