hdu1568（运用对数log10来解决斐波纳契数组【即(f[0]=0,f[1]=1;f[i] = f[i-1]+f[i-2](i>=2))】问题）

题目链接http://acm.hdu.edu.cn/showproblem.php?pid=1568
题目大意：
本文参考博文为：https://blog.csdn.net/queuelovestack/article/details/47859387
题目思路：

这道题，不能直接硬算，不是爆内存就是爆时间。
这道题需要用到网上的斐波纳契数组的公式：

将上面的式子进行如下的化简
首先说一下为什么要取对数log10，以及这个操作的好处。
比如我们需要求一个aⁿ的数的前4位,a为常数，而当a较大时，我们无法直接将其数值求出来。
那我们可以令M=aⁿ，
两边取对数log10,就化成了log10M=nlog10a（a为常数）
再将左边的log10移到右边去，得M=10^n*log10a,
然后我们可以注意到，nlog10a算出分为整数部分A和小数部分B，
而整数部分呢，10^A无论怎么算都是首位为1，无法影响最后的值。（即不可能让数首位变成2，3，4，5…等等）
那么也就是说只有小数部分会影响这个数的位数上的值。
所以，M=aⁿ这个数要求其个数位的值，
M=10^{(取小数部位的)nlog10a}【注意,这样求出来的值一定是X.XXX的，也就是整数位数只有一个数】
如果需要求其前四位，就相当于是将上面的结果乘上若干个10将整数部位拉成4位数记为ans。
然后最后输出的时候将ans取整即可得到aⁿ的前四位数了
【注意：为什么上面是乘上若干个10，而不是直接乘以1000】
因为有可能这个数是0.111或者是0.0011等之类的，单单是往前面移动3个数是不够的，所以得是若干个乘以10，相当于是移动若干个位

所以按照上面的思路来看我们需要将f[n]式子左右取对数log10，然后进行化简处理。

然后上面那个式子最右边为什么取约等于呢？
因为将下面这个式子给删掉了。所以要取约等于。

而删掉这个式子的原因就是，当n趋向比较大的时候，我们会发现((1-√5)/(1+√5))ⁿ这个式子是趋向0的。所以我们可以给他去掉，方便计算。同时这个也会给我们带来一个负面影响，就是当n取值不大的时候就会导致精度不够导致算错。这也就是我们为什么f[]数组前20个要自己算的原因。
同时前20个自己算，正好符合小于等于4位数直接输出。如下所示：
f[0]:0
f[1]:1
f[2]:1
f[3]:2
f[4]:3
f[5]:5
f[6]:8
f[7]:13
f[8]:21
f[9]:34
f[10]:55
f[11]:89
f[12]:144
f[13]:233
f[14]:377
f[15]:610
f[16]:987
f[17]:1597
f[18]:2584
f[19]:4181
f[20]:6765
==========在这儿之前都是小于等于四位的
f[21]:10946
f[22]:17711
f[23]:28657
f[24]:46368
f[25]:75025
f[26]:121393
…

ac代码：

#include<iostream>
#include<stdio.h>
#include<algorithm>

using namespace std;

#define N 100000005
typedef long long ll;

int f[30];


int main(){
	
	f[0] = 0;
	f[1] = 1;
	for(ll i = 2;i <= 20;i++){//直接算出前20个的值 
		f[i] = f[i-1] + f[i-2];
	}
	ll n;
	double temp;
	int ans;
	
	while(scanf("%lld",&n) != EOF){
		if(n<=20){
			printf("%d\n",f[n]);
			continue;
		}
		temp = log10(1.0/sqrt(5.0))+n*log10((1.0+sqrt(5.0))/2.0);
		temp = temp - (int)temp;//得到小数部分，因为其的前面位数，
								//只与小数部分有关
		temp = pow(10,temp);//这个算出来的首位肯定只有一位整数 
		
		ans = (int)temp;//将temp的整数部分给ans 
		
		while(ans<1000){//如果不是四位就将位数前移 
			temp = temp*10;
			ans = (int)temp;
		}
		
		printf("%d\n",ans);
	}
	
	
	return 0;
	
	
}