欧拉筛-质数筛的神c++

szm111213

于 2024-09-08 10:48:34 发布

阅读量730

点赞数 31

文章标签：算法 c++ 数据结构

本文链接：https://blog.csdn.net/szm111213/article/details/142023414

版权

欧拉筛（Euler's Sieve），又称线性筛法或欧拉线性筛，是一种高效筛选素数的方法。其核心思想是从小到大遍历每个数，同时标记其倍数为合数，但每个合数只被其最小的质因数标记一次，从而避免了重复标记，实现了线性时间复杂度的素数筛选。

欧拉筛的基本步骤

初始化：
- 初始化一个布尔数组isprime，其中isprime[i]表示自然数i是否为素数。
- 初始化一个整数数组prime，用来保存所有的素数。
- 初始化整数n，表示筛选范围n以内的所有素数。
- 初始化整数cnt，表示当前筛选的素数的个数，同时作为prime数组的下标进行填充。
- 将isprime[1]设置为false，因为1不是素数。
遍历与筛选：
- 从2开始循环遍历到n，对于每一次循环：
  - 如果isprime[i]为true，说明i为素数，将i添加到prime数组中。
  - 接下来，使用一个内部循环，从j=1到cnt（即当前prime数组的大小），同时要求i*prime[j] <= n。将isprime[i*prime[j]]设置为false，标记i*prime[j]为合数。
  - 在内部循环中，如果i % prime[j] == 0，则立即退出该内部循环。这是因为prime[j]已经是i的最小质因数，继续筛选会造成重复，并增加不必要的计算。

欧拉筛的关键点

每个合数只被筛一次：这是欧拉筛算法的核心优势，它避免了像埃拉托斯特尼筛法（Sieve of Eratosthenes）那样的重复筛选，从而实现了线性时间复杂度。
最小质因数原则：在内部循环中，通过检查i % prime[j] == 0来确保每个合数只被其最小的质因数筛掉。

欧拉筛的用途

欧拉筛法不仅限于筛选素数，它还可以用于计算其他与素数相关的数论函数，如欧拉函数φ(x)（表示小于x的、与x互质的数的个数）、1n之间每个数的因子个数、1n之间每个数的因数和等。

示例代码（C++风格）

#include <vector>  
#include <iostream>  
using namespace std;  
  
void EulerSieve(int n) {  
    vector<bool> isprime(n + 1, true);  
    vector<int> prime;  
    isprime[1] = false;  
  
    for (int i = 2; i <= n; i++) {  
        if (isprime[i]) {  
            prime.push_back(i);  
        }  
        for (int j = 0; j < prime.size() && i * prime[j] <= n; j++) {  
            isprime[i * prime[j]] = false;  
            if (i % prime[j] == 0) break;  
        }  
    }  
  
    // 输出素数  
    for (int p : prime) {  
        cout << p << " ";  
    }  
}  
  
int main() {  
    int n = 100;  
    EulerSieve(n);  
    return 0;  
}

这段代码展示了如何使用欧拉筛算法来筛选并输出1到100之间的所有素数。

欧拉筛的优化

实际上，欧拉筛算法本身已经非常高效，其时间复杂度为O(n)，空间复杂度也为O(n)（考虑到需要存储所有的素数）。然而，在实际应用中，我们可以根据具体需求对算法进行微调。

例如，如果只需要知道小于等于n的素数个数，而不需要知道具体的素数列表，我们可以使用一个计数器来代替prime数组。

欧拉筛的变种

欧拉筛的思想可以扩展到计算其他数论函数。以下是一些可能的应用：

计算欧拉函数φ(n)：
欧拉函数φ(n)定义为小于等于n的正整数中与n互质的数的个数。利用欧拉筛，我们可以在筛选素数的同时计算出每个数的欧拉函数值。
计算每个数的最小质因数：
在筛选过程中，我们可以额外维护一个数组来记录每个数的最小质因数。这对于后续的一些数论计算非常有用。
计算每个数的质因数分解：
虽然欧拉筛本身不直接用于计算质因数分解，但它可以为质因数分解提供一个高效的素数列表，从而加速质因数分解的过程。

欧拉筛与其他筛法的比较

埃拉托斯特尼筛法（Sieve of Eratosthenes）：
埃拉托斯特尼筛法是最直观的素数筛选方法，但其时间复杂度为O(n log log n)，且会重复标记合数。
埃拉托斯特尼筛法的优化版本（也称为“分段筛法”或“区间筛法”）：
当需要筛选的范围非常大，但只关心某个区间内的素数时，分段筛法可以显著提高效率。然而，它仍然不如欧拉筛那样简洁和高效。
轮询筛法（Wheel Sieve）：
轮询筛法是一种基于埃拉托斯特尼筛法的优化方法，它通过跳过一些显然不是素数的数（如偶数、3的倍数等）来减少不必要的检查。但它仍然无法避免重复标记的问题。