从埃筛到欧拉筛

埃式筛法是第一个比较先进的筛法,以后的所有筛法大多都是以埃筛为基础推进的
埃筛的思想也很简单,筛掉所以素数的倍数,最后还剩的就是素数

bool isprime[maxn];
void Eratosthenes_prime(){
    memset(isprime,true,sizeof isprime);
    isprime[0]=isprime[1]=false;
    for(int i=2;i<maxn;i++){
        if(isprime[i])
            for(int j=i+i;j<maxn;j+=i)
                isprime[j]=false;
    }
}

这个筛法的唯一缺陷是被重复筛去的元素太多,导致算法不够精妙
比如6会被2筛一次,还会被3筛一次,这就导致埃筛的时间只能达到O(nlogn),这样的复杂度对付小数据,要是来个1e6是绝对不可能过的
再这之上,欧拉筛法就出现了,它避免了重复筛这一情况

bool isprime[maxn];
int prime[maxn/10],tot;
void Euler_prime(){
    memset(isprime,true,sizeof isprime);
    isprime[0]=isprime[1]=false;
    for(int i=2;i<maxn;i++){
        if(isprime[i])prime[tot++]=i;
        for(int j=0;j<tot&&prime[j]<=maxn/i;j++){
            isprime[i*prime[j]]=false;
            if(i%prime[j]==0)break;
        }
    }
}

欧拉筛唯一亮点所在在于if(i%prime[j]==0)break;这句,是从埃筛提升到欧拉筛的关键
欧拉筛的时间是O(n)。

上面这个欧筛中的isprime数组其实是作用不大的,我们可以把它换成欧拉函数,适当的我们也可以添加莫比乌斯函数

int prime[maxn/10],phi[maxn],mu[maxn],tot;
void Euler(){
    phi[1]=1;
    for(int i=2;i<maxn;i++){
        if(!phi[i]){
            phi[i]=i-1;
            prime[tot++]=i;
            mu[i]=-1;
        }
        for(int j=0;j<tot&&prime[j]*i<maxn;j++){
            if(i%prime[j]==0){
                phi[i*prime[j]]=phi[i]*prime[j];
                mu[i*prime[j]]=-mu[i];
                break;
            }
            phi[i*prime[j]]=phi[i]*(prime[j]-1);
            mu[i*prime[j]]=0;
        }
    }
}

除此之外,比较先进的素筛算法还有Min_25筛,杜教筛,洲阁筛
三种筛法时间都在O(n)之下,有些出题人就喜欢把数据故意整到1e10,就看你会不会这三种先进筛法

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Macarons_i

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值