素数线性筛——欧拉Euler筛

时间复杂度O(n)

prime[]数组中的素数是递增的,当i能整除prime[j]，那么i*prime[j+1]这个合数肯定被prime[j]乘以某个数筛掉。因为i中含有prime[j],prime[j]比prime[j+1]小，即i=k*prime[j]，那么i*prime[j+1]=(k*prime[j])*prime [j+1]=k’*prime[j]，接下去的素数同理。所以不用筛下去了。因此，在满足i%prime[j]==0这个条件之前以及第一次满足改条件时,prime[j]必定是prime[j]*i的最小因子。

bool vis[maxn];
int prime[maxn];
void isprime()
{
    memset(vis,true,sizeof(vis));
    int c=0;
    for(int i=2;i<=maxn;i++)
    {
        if(vis[i])
        {
            prime[c++]=i;
        }
        for(int j=0;j<c&&i*prime[j]<=maxn;j++)
        {
            vis[i*prime[j]]=false;
            if(i%prime[j]==0)//不中断会将prime[j+1]的倍数筛掉，造成重复
            {
                break;
            }
        }
    }
}

对于 i%prime[j] == 0 就break的解释 ：当 i是prime[j]的倍数时，i = kprime[j]，如果继续运算 j+1，i * prime[j+1] = prime[j] * k prime[j+1]，这里prime[j]是最小的素因子，当i = k * prime[j+1]时会重复，所以才跳出循环。
举个例子 ：i = 8 ，j = 1，prime[j] = 2，如果不跳出循环，prime[j+1] = 3，8 * 3 = 2 * 4 * 3 = 2 * 12，在i = 12时会计算。因为欧拉筛法的原理便是通过最小素因子来消除。