线性筛素数

最新推荐文章于 2024-07-05 00:10:08 发布

zzkya

最新推荐文章于 2024-07-05 00:10:08 发布

阅读量546

点赞数 6

分类专栏：数学

本文链接：https://blog.csdn.net/loi__zzk/article/details/52985738

版权

数学专栏收录该内容

1 篇文章 0 订阅

订阅专栏

给定一个范围N，你需要处理M个某数字是否为质数的询问（每个数字均在范围1-N内）

i%prime[j]==0是O(n)筛素数的关键
为什么这里要break？prime[j]为（i*prime[j]）的最小质因数，比一个合数数大的质数和该合数的乘积可用一个更大的合数和比其小的质数相乘得到，所以break。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
using namespace std;
const int maxn=1000005;
int prime[maxn];
bool not_prime[maxn];
int main()
{
    int n,cnt=0;
    scanf("%d",&n);
    memset(not_prime,0,sizeof(not_prime));
    not_prime[1]=true;
    for(int i=1;i<=n;i++)
    {
        if(!not_prime[i])
            prime[++cnt]=i;
        for(int j=1;j<=cnt;j++)
        {
            if(prime[j]*i>n) break;
            not_prime[prime[j]*i]=true;
            if(i%prime[j]==0) break;
        }
    }
    for(int i=1;i<=cnt;i++)
        printf("%d ",prime[i]);

    return 0;
}