素数筛法（欧拉筛、埃氏筛）

最新推荐文章于 2023-09-18 16:28:37 发布

朝花夕拾喵

最新推荐文章于 2023-09-18 16:28:37 发布

阅读量262

点赞数

分类专栏：算法学习以及OJ题目原码文章标签：算法数学

本文链接：https://blog.csdn.net/Azxz0_0/article/details/120603481

版权

算法学习以及OJ题目原码专栏收录该内容

9 篇文章 0 订阅

订阅专栏

本文详细介绍了两种常用的素数筛法——埃氏筛和欧拉筛。埃氏筛基于质数的倍数不是质数的原理，而欧拉筛则优化了这一过程，确保每个合数仅被其最小质因子筛去，提高了效率。两种筛法的时间复杂度分别为O(nlog^2_2n)和O(n)，并提供了C语言的实现代码。此外，文章还提到可以使用std::bitset进行空间优化，并鼓励读者将筛法应用于其他问题的求解。

摘要由CSDN通过智能技术生成

一、埃氏筛

1.原理：

质数的倍数不是质数。

时间复杂度： $O(nlog^2_2n)$
埃氏筛的时间复杂度不会证明，可戳别处

2.代码

int prm[255], prm_cnt;//存放素数的数组、素数个数
bool iscompos[255];//true代表为合数，false代表素数数
void make_prime(int n) {
	if (n <= 1)return;
	for (int i = 2; i <= n; i++) {
		if (iscompos[i] == 0) {
			prm[++prm_cnt] = i;
			for (int j = i * i; j <= n; j += i) {
				iscompos[j] = 1;
			}
		}
	}
	return;
}

二、欧拉筛

1.前言

埃氏筛中，例如21会被3、7同时筛去。

我们希望每个合数仅被筛去一次以提高效率。于是有了欧拉筛。

2.原理

每个合数仅被最小质因子/最大因数筛去

时间复杂度： $O (n)$
每个数仅被处理一次故时间复杂度为 $O (n)$

3.代码

//时间复杂度：O(n)
int prm2[maxn], cnt2, tab2[maxn];
void make_prime2(int x) {
    if (x <= 1)return;
    for (int i = 2; i <= x; i++) {
        if (!tab2[i])
            prm2[++cnt2] = i;
        for (int j = 1; j <= cnt2 && prm2[j] * i < x; j++) { // 枚举素数，筛去倍数
            tab2[i * prm2[j]] = 1;
            if (!(i % prm2[j])) // 当前数为该素数的倍数，即 prm[j] 为最小质因子
                break;    
        }
    }
}