nyoj 461 Fibonacci数列(四)

原创 2012年03月27日 20:36:13

Fibonacci数列是个神奇的数列,它包含了太多的知识,有些你甚至想都想不到……

这个题目的意思很简单了,但是这个题目绝对不是让你循环算出的,或许你做过fibonacci数列(二)的话,你也许在想使用矩阵二分幂的方法,但是这个方法只适合找F(n)的低位数字,至于这个题目,则是另一种思路,使用通项公式。

实际上我们可以把这一类题目归类:求最Fibonacci数列F(n)的高位,求n^m的高位,求n!的高位……,求n!有多少位数……,这类为题是要找一个很大的数的高位,或者估计出这个数有多大。区别于以前的求低位,这些问题实际上是一种近似计算,我们一般的解法是找出欲求的表达式an的同阶的一个表达式bn;

那么n很大的时候我们就可以用bn来近似代替an了,这里我们一般要求bn是比较容易计算出来的表达式。

具体到这个题目,我们知道Fibonacci数列有通项公为:

 

 

 

 

上面通项公式中间应该为减号,copy过来的图片,改不了,,,无语

那么问题至此就很简单了,我们最后可以计算出log(an),保留它的的小数部分为ans即可。

#include<stdio.h>
#include<math.h>
int main()
{
	int n,m;
	int p[]={0,1,1,2,3,5,8,13,21,34,55,89,144,233,377,610,987,1597,2584,4181,6765};
	while(scanf("%d",&n)!=EOF)
	{
		double sum,a=-0.5*log(5)/log(10),b=log( ( 1+sqrt(5.0) )/2 )/log(10);
		if(n<=20)
			printf("%d\n",p[n]);
		else
		{
			double sum;
			sum=a+n*b;
			sum=sum-int(sum);
			sum=pow(10,sum);
			while(sum<1000)
				sum=sum*10;
			printf("%d\n",int(sum));

		}
	}
}




相关文章推荐

nyoj461 Fibonacci数列(4)解通项公式

Fibonacci数列(四) 时间限制:1000 ms  |  内存限制:65535 KB 难度:4 描述 数学神童小明终于把0到100000000的Fibonacci数列(...

NYOJ109 数列转换 【守恒法】

数列转换 时间限制:3000 ms  |  内存限制:65535 KB 难度:3 描述 有一个数列a1,a2,a3...an,每次可以从中任意选三个相邻的数ai-1 ,...

NYOJ 461-Fibonacci数列(四)(求斐波那契数列前4位)

题目地址:NYOJ 461 思路:斐波那契数列的通项公式为 然后下一步考虑如何产生前4位: 先看对数的性质,loga(b^c)=c*loga(b),loga(b*c)=loga(b)+loga(c);...

NYOJ 461 又见Fibonacci数列

原题链接:http://acm.nyist.net/JudgeOnline/problem.php?pid=461 这个题一看貌似和以前的斐波那契数列一样  就是数据可能大些。。。于是就认为...

nyoj 148 fibonacci数列(二)【矩阵】

fibonacci数列(二) 时间限制:1000 ms  |  内存限制:65535 KB 难度:3 描述 In the Fibonacci integer sequence, F0...

nyoj 148 fibonacci数列(二)

题目链接:http://acm.nyist.net/JudgeOnline/problem.php?pid=148 矩阵快速幂入门题,身为连线代都没学过的渣渣,为了学矩阵快速幂,还得百度矩阵学习矩阵...
  • AcmLzq
  • AcmLzq
  • 2016年05月07日 17:13
  • 455

nyoj fibonacci数列(二) 矩阵乘法

http://acm.nyist.net/JudgeOnline/problem.php?pid=148 fibonacci数列(二) 时间限制:1000 ms  |  内存限制:65...

NYOJ 148 fibonacci数列(二)

fibonacci数列(二)

nyoj 301 递推求值 和 nyoj 148 fibonacci数列(二) 【矩阵】

原题链接: nyoj 301 :http://acm.nyist.net/JudgeOnline/problem.php?pid=301 nyoj 148 :http://acm.nyist.ne...

acm-Fibonacci数列(四)

Fibonacci数列(四) 时间限制:1000 ms  |  内存限制:65535 KB 难度:4 描述  数学神童小明终于把0到100000000的Fibonacci数列(f[0]=...
  • zzyy17
  • zzyy17
  • 2013年11月12日 16:43
  • 447
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:nyoj 461 Fibonacci数列(四)
举报原因:
原因补充:

(最多只允许输入30个字)