PJU2480-Longge's problem—算法10.6 计算5000以内的素数，筛法

zjutsunny

于 2019-11-12 22:14:18 发布

阅读量167

点赞数

分类专栏：算法分析

本文链接：https://blog.youkuaiyun.com/zjutsunny/article/details/103039719

版权

算法分析专栏收录该内容

149 篇文章

订阅专栏

本文详细介绍了使用素数筛算法生成大量素数的方法，并通过实例展示了如何利用生成的素数表计算乘法逆元，这对于解决数学问题及提高编程竞赛中的算法效率具有重要意义。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include <stdio.h>

#define MAX 50000
__int64 prime[MAX];
int pnum = 0;

void primtable()
{
	bool flag[MAX] = {0};
	for (int i=2; i<MAX; i++)
	{
		if (!flag[i])
		{
			prime[pnum++] = i;
			for (long j = i+i; j<MAX; j += i) 
				flag[j] = 1;
		}
	}
}

int main()
{
	primtable();
	int n;
	while (scanf("%d", &n) != EOF) 
	{
		int i;
		__int64 ans = n;
		int p, a;
		for (i=0; prime[i]*prime[i]<=n; ++i) 
		{
			p = prime[i];
			if (n%p == 0) 
			{
				a = 0;
				while (n%p == 0) 
				{
					a++;
					n /= p;
				}
				ans = ans + ans*a*(p-1)/p;
			}
		}
		if (n!=1) ans = ans + ans*(n-1)/n;
		printf("%I64d\n", ans);
	}
	return 0;
}