埃式筛法使用的场景:
例如输入一个数n,判断他是不是素数(有多组测试数据)n的范围是 (n<1e6)
我们知道用普通的判断一个素是否为素数的方法,时间复杂度为O(sqrt(n))次,如果测试数据有k组,时间复杂度为O( )
用埃式筛法可以先预处理出来 1e6以内的所有素数,用prime数组记录,需要O(nlog(log n))次,后面判断一个数是否是素数只需要查数组,只需要O(1)的复杂度
所以当测试数据k很小时,用普通方法;当测试数据很多时,用埃式筛法。
例如求1到20的素数个数,从2开始
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
用一个prime[21]数组,如果prime[i]=1表示i不是素数,prime[i]=0表示i是素数,一开始初始化prime为全0,筛选n轮后prime[i]=0即为素数
由于prime[2]=0,从2开始,所以筛掉2的所有倍数并prime置1,prime为0的剩下以下数
3 5 7 9 11 13 15 17 19
prime[3]=0,所以第二次从3开始,所以筛掉3的所有倍数并prime置1,prime为0的剩下以下数
5 7 11 13 17 19
prime[5]=0,所以第三次从5开始,所以筛掉5的所有倍数并prime置1,prime为0的剩下以下数
7 11 13 17 19
…
最后等于0的有:prime[2]=0,prime[3]=0, prime[5]=0, prime[7]=0, prime[11]=0,prime[13]=0,prime[17]=0,prime[19]=0
关键代码实现
for(i=2;i<N;i++){
if(prime[i]==0){ //只有prime等于0才要进入循环
int j=2*i; //一开始为i的两倍
for(j=i*2;j<N;j=j+i) //j每次加i,即3倍,4倍,5倍...
prime[j]=1;
}
}
for(i=1;i<N;i++){ //求1到20的素数个数
int ans;
if(prime[i]==0)
ans++;
}