【HDU1568】fib数列公式的使用

1.题目链接这个题目的意思就是输出fib数列每一项的前四位数字。这个题应该没有什么好的做法，只能用fib数列的通项公式了。

2.fib的通项公式：

 

                                             

至于这个公式是怎么来的呢？三言两语也说不清。可能有很多别的推导办法吧，但是我所知道的时使用数列递推的特征方程，本质上就是矩阵的特征方程和特征根来解决这个问题。有兴趣的可以去看一下（其实我觉得大家应该都会吧，高中讲数列的时候老师一般都会把这个通项公式告诉大家，因为他是在时太特殊了，虽然每一项都是无理数，但是计算的结果却是一个有理数）。OK，知道了这个就很好办了，我们直接根据这个公示变一下形，两边同时取对数：

     

这个式子在fib比较大的时候时成立的，因为n比较大的时候我们知道后边的一项是很小的一个数，对我们找前四项不会产生大的影响。所以问题解决，代码如下：

#include<bits/stdc++.h>
using namespace std;
#pragma warning(disable:4996)
#define ll long long
const int N = 20;//这里为什么设置成20，因为第20项已经不再是4位数，并且这个时候使用公式的精度已经足够
ll f[N];
int main()
{
	int n, i;
	double s;
	f[0] = 0, f[1] = 1;
	for (i = 2; i<N; i++)
		f[i] = f[i - 1] + f[i - 2];
	while (~scanf("%d", &n))
	{
		if (n<N)
		{
			printf("%d\n", f[n]);
			continue;
		}
		s = log10(1.0 / sqrt(5.0)) + n * log10((1 + sqrt(5.0)) / 2);
		s = s - (int)s;
		s = pow(10, s);
		while (s<1000)
			s *= 10;
		printf("%d\n", (int)s);
	}
	return 0;
}