线性筛素数

给定一个范围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;
}
评论 3
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值