线性筛

Posted 2021-01-03 wtsruvf

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的最小因子。
代码是斌神的  解释是这个大佬的：https://www.cnblogs.com/Miroerwf/p/7776390.html

const int maxn = 10010, INF = 0x7fffffff;
int prime[maxn+1];
void get_prime()
{
    mem(prime, 0);
    for(int i=2; i<=maxn; i++)
    {
        if(!prime[i]) prime[++prime[0] = i;
        for(int j=1; j<=prime[0] && prime[j] <= maxn/i; j++)
        {
            prime[prime[j]*i] = 1;
            if(i % prime[j] == 0) break;
        }
    }
}