素数的线性欧拉筛

先上代码:


#define maxx 400005
using namespace std;
int prime[100005];
bool p[maxx];
int cnt=0;
void init()
{
    for(int i=2;i<maxx;i++)
    {
        if(!p[i])
            prime[cnt++]=i;
        for(int j=0;j<cnt;j++)
        {
            if(i*prime[j]>=maxx)break;
            p[i*prime[j]]=true;
            if(i%prime[j]==0)break;
        }
    }
}

如何理解这个代码只需要知道一点,每个合数只会被他的最小的质因子筛去。
证明:
我们设有合数写成x=pqx=pq,pp为质数,且为x的最小质因子。假设当i在遍历到xx之前肯定已经遍历到q,从而x=pqx=p∗q必然已经被打上了标记。每个数的值只会被修改一次(即被最小的质因子筛去),原因是判断了q%p==0q%p==0,因为我们每次是从已经筛出的质因子里从小到大考虑,所以遇到的第一个pp,若满足q%p==0,则一定是qq的最小质因子。
假设没有这个条件,循环没有break,那么假设当我们遇到了第二质数p满足q%p==0q%p′==0,则x=pqx′=p′q会被修改一次,但是pp′并不是xx′的最小质数,所以就重复修改了。

xx′会被筛去么?会的,因为p|xp|x′,设q=xpq′=x′p,当i遍历到qq′时,便会把x=pqx′=pq′筛去,这在逻辑上是能自恰的。所以这样证明了每个合数只会被他最小的质因子筛去。

因为数组的每个数只会被修改一次,所以复杂度就是O(n)O(n)啦。

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值