素数的三种筛法比较

最新推荐文章于 2022-06-06 15:48:42 发布

SF-_-

最新推荐文章于 2022-06-06 15:48:42 发布

阅读量1.8k

点赞数

分类专栏：数学 ACM 文章标签：素数

本文链接：https://blog.csdn.net/sf____/article/details/24844545

版权

ACM 同时被 2 个专栏收录

224 篇文章 0 订阅

订阅专栏

数学

96 篇文章 0 订阅

订阅专栏

ACM比赛中经常能遇到关于素数的题目，可能需要我们预处理求出100W，1000W以内的所有素数。笔者在这里介绍共三种筛素数的方法，代码如下：

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

const int maxn=100000000;
int prime[maxn];
int primeNum;
bool h[maxn];

// 普通版
void calPrime()
{
    for(int i=2;i<maxn;i++) if(!h[i])
    {
        prime[primeNum++]=i;
        for(int j=i+i;j<maxn;j+=i)
            h[j]=true;
    }
    printf("%d\n", primeNum);
}

// 优化版
void calPrime2()
{
    for(int i=2;i<maxn;i++) if(!h[i])
    {
        prime[primeNum++]=i;
        if(maxn>(long long)i*i) for(int j=i*i;j<maxn;j+=i)
            h[j]=true;
    }
    printf("%d\n", primeNum);
}

// 欧拉筛法
void calPrime3()
{
    for(int i=2;i<maxn;i++)
    {
        if(!h[i]) prime[primeNum++]=i;
        for(int j=0;j<primeNum;j++)
        {
            if(i*prime[j]>maxn) break;
            h[i*prime[j]]=true;
            if(i%prime[j]==0) break;
        }
    }
    printf("%d\n", primeNum);
}

int main()
{
    calPrime();
//    calPrime2();
//    calPrime3();
}

在笔者的电脑上，求1亿以内的所有素数，普通方法大概4.4S，优化版大概2.9S，欧拉筛法大概1.5S，优化程度还是很高的。

当然，只求100W以内的所有素数三种方法基本无差别。