Eratosthenes筛法(素数筛)

最一般的素数筛思想很简单,对于不超过maxx的每个非负整数p,删除2p,3p,4p,……,然后剩下的就是素数,复杂度O(nlogn),因为对内层循环n/2+n/3+......+n/n  小于  1+1/2+1/3+...+1/n=ln(n+1)+γ,其中γ为欧拉常数≈0.577218

应当注意,1不是素数哦~

这样已经不慢,但由于所有非素数都能由素数与素数的乘积表示,又可以优化如下

int maxx=1000000;
long long prime[maxx];
bool isPrime[maxx];
void Eratosthenes()
{
    memset(isPrime,1,sizeof(isPrime));
    isPrime[1]=0;
    for(int i=2;i<=maxx;++i)
    {
        if(isPrime[i])
        {
            for(int j=2*i;j<=maxx;j+=i)
                isPrime[j]=0;
        }
    }
}


时间复杂度是O(nloglogn)

但是还能再优化,内层循环不必从i*2开始,它已经在i=2时候就被筛掉了。优化如下:

int maxx=1000000;
long long prime[maxx];
bool isPrime[maxx];
void Eratosthenes()
{
    memset(isPrime,1,sizeof(isPrime));
    isPrime[1]=0;
    int m=sqrt(maxx+0.5)
    for(int i=2;i<=m;++i)
    {
        if(isPrime[i])
        {
            for(int j=i*i;j<=maxx;j+=i)
                isPrime[j]=0;
        }
    }
}




评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值