HDU 1568 Fibonacci

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

Input 输入若干数字n(0 <= n <= 100000000)，每个数字一行。读到文件尾。

Output             输出f[n]的前4个数字（若不足4个数字，就全部输出）。

Sample Input 0 1 2 3 4 5 35 36 37 38 39 40

Sample Output 0 1 1 2 3 5 9227 1493 2415 3908 6324 1023

Author daringQQ

Source Happy 2007

Recommend 8600

思路：n(0 <= n <= 100000000)的数据范围，递推公式不适合

考虑O（1）公式

先对10取对数，可以等到一个数numdf,由于log10(F[n])=numdf,所以F[n]=pow(10.0,numdf);

由于要求的是前四位数，先把F[n]处理成为小数；

参考：log10(12345)=log10(1.2345*10^4)=log10(1.2345)+4=n(4.09149）;

那么10^n=12345,但我们需要的是1.2345，那么只需要求出n的小数部分n1（0.09149）=(double)n-(int)n就可以了；

注：10^(4.09149)=10^4*10^0.019149=10000*1.2345;

下面是代码

#include<algorithm>
#include<iostream>
#include<cmath>
using namespace std;
const double fi=sqrt(5.0);
int Fibon[21];//F[20]=6765,F[21]=10946>1000
//Fibon公式：F(n)=[((1+√5)/2)^n-((1-√5)/2)^n]/√5
int main()
{
	double numdf;
	int numif;
	Fibon[0]=0;
	Fibon[1]=1;
	for(int i=2;i<21;i++){
		Fibon[i]=Fibon[i-1]+Fibon[i-2];
	}
	int n;
	while(cin>>n){
		if(n<=20){
			cout<<Fibon[n]<<endl;
			continue;
		}
		numdf=n*log10((1+fi)/2.0)-log10(fi);//由于对数运算中 n-((1-√5)/2)^n十分小，忽略不计 
		numif=(int)numdf;
		double numf=numdf-numif;
		double num=pow(10.0,numf);
//		double k=pow(10.0,numdf);
//		cout<<numdf<<"  "<<k<<endl;
//		cout<<numf<<"  "<<num<<endl;
		while(num<1000){
			num*=10;
		}
		cout<<(int)num<<endl;
	}
}