首页 > 代码库 > 筛素数

筛素数

普通筛

const int MAXN = 50000void 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]        }    }}

 

筛素数