先说说埃氏筛
bool numlist[100000001];
int prime[20000001], cnt;
void work(int n){
for(int i=2; i<=n; i++){
if(numlist[i]==false){
prime[++cnt] = i ;
for(int j=i; i*j<=n; j++)
numlist[i*j] = true;
}
}
return;
}
埃氏筛的思想就是把不大于根号n的所有素数的倍数剔除,剩下的就是素数。
就像2这个素数,把不包括2在内的2的倍数全部剔,然后再向后枚举,查到一个未被剔除的数,再进行上述操作,比如2下一个数就是3,然后再把3的倍数剔除。但是这样的做法会造成重复剔除,比如2剔除了12,然后3又剔除一次12,接着6又来一次,就降低了效率。
所以我们要想避免重复剔除,将要用到另一种筛选方法-欧拉筛。
bool numlist[100000001];
int prime[20000001], cnt;
void work(int n){
for(int i=2; i<=n; i++){
if(numlist[i]==false)
prime[++cnt] = i ;
for(int j=1; j<=cnt&&i*prime[j]<=n; j++){
numlist[i*prime[j]] = true ;
if(i%prime[j]==0)
break;
}
}
return;
}
避免重复筛,应找到筛合数的一种原则:这个合数只会被它的最大非自身因数(对应最小质因数)筛。这样能保证每个合数只会被筛一次。
(运行过程)从 2 开始,2 加入 prime 数组,再从小到大枚举质数(现在只有 2),筛掉质数与 2 的乘积(4 被筛掉)。
#### 到了 3,3 加入 prime 数组,从小到大枚举质数(此时有 2,3),筛掉质数与 3 的乘积(6,9 被筛掉)。 到了 4,4 没加入 prime 数组,枚举质数(有2,3),筛掉 8 后,因为 4 mod 2=0,触发退出条件。(不触发,就会筛掉 12,而 12又会被2和6筛一次)