筛素数

Posted 2020-09-19 bennettz

普通筛

const int MAXN = 50000；  
void Prime()    // 素数表p[i]=1为素数 
{  
    int i, j;  
    for (i=0; i<MAXN; i++) prime[i] = 1;  
    prime[0] = prime[1] = 0;  
    for (i=2; i<MAXN; i++)  
    {  
        if (!prime[i]) continue;  
        for (j=i*2; j<MAXN; j+=i) prime[ j ] = 0; //将每个素数的倍数都删去（有重复）
    }  
}  

线性筛

每个合数除1外最小的因数一定是素数，只利用这个素数筛没有重复

int prim[50000],prime_n; 
bool is_pr[500000]={1,1};
void prime(){
    int m=500000;
    for(int i=2;i<m;i++){
        if(!is_pr[i])prim[prime_n++]=i;//没有筛掉的为素数
        for(int j=0;j<prime_n&&i*prim[j]<m;j++){
            is_pr[prim[j]*i]=1;//筛掉
            if(i%prim[j]==0)break;//后面的数最小因数一定小于等于prim[j]
        }
    }
}