欧拉筛法和埃氏筛法

最新推荐文章于 2022-08-02 12:08:10 发布

猫的玖月

最新推荐文章于 2022-08-02 12:08:10 发布

阅读量164

点赞数 1

分类专栏： ACM竞赛（C++）文章标签：算法

本文链接：https://blog.csdn.net/qq_44061866/article/details/105322740

版权

ACM竞赛（C++）专栏收录该内容

126 篇文章 4 订阅

订阅专栏

在这里插入图片描述

void FilterPrime(int n){
    bool* isPrimes = new bool[n + 1];//动态数组；释放空间：delete[] isPrimes;
    for(int i = 2; i <= n; i++){
        isPrimes[i] = true;
    }
    for(int j = 2; j <= n; j++){
        if(isPrimes[j] == true){
            for(int k = 2; k * j <= n; k++){
                isPrimes[k * j] = false;
            }
        }
    }

}

问题：埃氏筛法复杂度为：O(nloglogn)，而在一些数据范围达到1e7这样的题目中，也很难让人满意，于是我便学习了欧拉筛法，也即 O(n)的线性筛法

欧拉筛法的基本思想：在埃氏筛法的基础上，它不会重复标记一个数是不是素数的倍数（所有合数都可以进行质因子分解，因此只要用已知素数筛选掉合数，剩下的就是素数）

#include <bits/stdc++.h>
using namespace std;
int prime[maxn]; //记录素数(和埃式筛法不同，下标不是素数，只是存了素数表)
bool visited[maxn]; //下标对应数字是否访问过
int Prime(int n){
    int cnt = 0;//小于n的素数
    memset(visited, false , sizeof(visited));
    memset(prime, 0 , sizeof(prime));
    for(int i = 2; i < n; i++){
        if(!visited[i]){
            prime[cnt++] = i;
        }
    
        for(int j = 0; j < cnt&& i * prime[j] < n; j++){
            visited[i*prime[j]] = true;
            if(i % prime[j] == 0) break; 
            /*
            证明：
            1.所有合数均可以表示成素数的积
            2. i % prime[j] == 0 可以化成 i == k * prime[j]
            3. 判断i是不是合数（素数倍数）
            */
        }
    }
    return cnt;//小于n的素数
}