欧拉筛模板
#define N 1000
bool check[N];//数组表示第i个数是否为素数,是为true,不是为false;
int prime[N];//数组存放素数序列,从0开始都是素数(递增序列)
int min_prime[N];//数组存放每个数的最小素因数
int Prime(int n){//返回n中素数的个数
memset(check,true,sizeof(check));//初始化,所有都为true,筛过后标记为false
int sum=0,i,j;//sum为素数序列中素数的个数。
for(i=2;i<=n;i++){
if(check[i])//当i为素数时
prime[sum++]=i;//将i存放到素数序列中
for(j=0;j<sum&&i*prime[j]<=n;j++){
check[i*prime[j]]=false;//把所有素数的i倍全部标记为false
min_prime[i*prime[j]]=prime[j];//存放最小素因数
if(i%prime[j]==0) //当i是prime[j]的倍数时,跳出循环
break;
}
}
return sum;
}
欧拉筛与普通筛法不一样的是所有合数都被其最小素因子筛去且仅筛一次。
关键代码就是
if(i%prime[j]==0) //当i是prime[j]的倍数时,跳出循环
break;
如果不跳出循环,那么下一步就是将 check[i * prime[j+1]] 变成 false,但是,如果 i 为 prime[j] 的倍数,即 i = k * prime[j].
那么,i * prime[j+1]=k * prime[j] * prime[j+1] = ( k * prime[j+1] ) * prime[j] = k' * prime[j]。当外层循环 i 增长到 k' 时,会将 k' * prime[j] 筛去(因为 prime[j]< prime[j+1])。