顾名思义,在线性的时间里筛选出质数,核心在于每次用合数的最小质因子筛去这个合数,时间复杂度为 O ( N ) O(N) O(N)。
void xxs()
{
int cnt=0;
for(int i=2;i<=3e6;i++)
{
if(!pd[i]) prime[++cnt]=i;
for(int j=1;j<=cnt&&i*prime[j]<=3e6;j++)
{
pd[prime[j]*i]=1;
if(i%prime[j]==0) break; // 核心
}
}
}
设 x = i ∗ p r i m e [ j ] ~x=i*prime[j] x=i∗prime[j],若 i % p r i m e [ j ] = 0 ~i\%prime[j]=0 i%prime[j]=0,说明 i i i 有这个质因子,而 j ~j~ j 从小到大遍历,因此 p r i m e [ j ] ~prime[j]~ prime[j] 就是 i ~i~ i 的最小质因子。
但此时若不跳出, j + 1 j+1 j+1 后,一定有 p r i m e [ j + 1 ] > p r i m e [ j ] prime[j+1]>prime[j] prime[j+1]>prime[j] ,对于 i ∗ p r i m e [ j + 1 ] i*prime[j+1] i∗prime[j+1] , i ~i~ i 中含有 p r i m e [ j ] prime[j] prime[j] 这个质因子了,因此 p r i m e [ j + 1 ] prime[j+1] prime[j+1] 并不是 i ∗ p r i m e [ j + 1 ] i*prime[j+1] i∗prime[j+1] 的最小质因子,所以必须跳出,这保证了每一个合数都是由其最小质因子筛掉。