UVa 10236 The Fibonacci Primes (斐波那契素数)

最新推荐文章于 2024-11-11 20:45:30 发布

synapse7

最新推荐文章于 2024-11-11 20:45:30 发布

阅读量1.7k

点赞数 1

分类专栏： UVa acm之路--数学数论文章标签： acm UVa c++

本文链接：https://blog.csdn.net/synapse7/article/details/16901299

版权

UVa 同时被 3 个专栏收录

407 篇文章 0 订阅

订阅专栏

acm之路--数学

262 篇文章 1 订阅

订阅专栏

数论

100 篇文章 0 订阅

订阅专栏

10236 - The Fibonacci Primes

Time limit: 3.000 seconds

http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&category=24&page=show_problem&problem=1177

注意此题的描述和维基百科上Fibonacci Prime的描述不同。

上面加着重符的文字说明：可以用类似素数筛的方式筛出Fibonacci Prime。

但是，定理“若b|a，则fb|fa”并不能推出“若p是质数，则fp是Fibonacci Prime”。

要证明这一点，我们可以证明“若fb|fa，则b|a”（详见《具体数学》P247及《Proofs that Really Count: The Art of Combinatorial Proof》P12）

完整代码：

/*0.022s*/

#include<cstdio>
const int maxn = 249450;

int fp[22001];
bool vis[maxn];

int main()
{
	long long a = 0, b = 1, tmp, i, j;
	int n, cp;
	for (cp = 1, i = 2; cp < 22001; ++i)
	{
		tmp = a + b, a = b, b = tmp;
		if (b >= 1e18) b /= 10, a /= 10;
		if (!vis[i])
		{
			tmp = b;
			while (tmp >= 1e9) tmp /= 10;
			fp[cp++] = tmp;
			for (j = i * i; j < maxn; j += i) vis[j] = true;
		}
	}
	//printf("%lld\n",i);
	fp[1] = 2, fp[2] = 3;
	while (~scanf("%d", &n))
		printf("%d\n", fp[n]);
	return 0;
}