二逼 普通 文艺的素数筛选法

二逼的素数筛选法

2b的筛选法，不多说。

const int N=1000000;
bool pr[N];
int p[N],lp;
bool isprime(int x){
    for(int i=2;i*i<=x;i++)
        if(x%i==0) return 0;
    return 1;
}
void initp(){
    for(int i=2;i<N;i++)
        if(isprime(i)){
            pr[i]=1;
            p[lp++]=i;
        }
}

普通的素数筛选法

用到了素数筛选的思想，但是很明显一些数被重复筛选。

普通筛选法1

const int N=100000;
bool pr[N];
int p[N],lp;
void initp()
{
    for(int i=2;i<N;i++)
     {
		if(!pr[i])
		{
            p[lp++]=i;
     
             for(int j=i+i;j<N;j+=i)
                  pr[j]=1; 
		}
      }
} 

普通筛选法2

const int N=100000;

bool pr[N];
int p[N],lp;
void initp()
{
	for(int i=4;i<N;i+=2)
		pr[i]=1;
    for(i=2;i<N;i++)
        if(!pr[i])
		{
            p[lp++]=i;
            if(N/i>=i)
                for(int j=i*i;j<N;j+=i*2)
                    pr[j]=1;
        }
}

普通筛选法3

void calprime(){
    long long i,j;
    memset(vis,0,sizeof(vis));
    priLen = 0;
    for(i = 2; i <= M; i++) 
        if(!vis[i]){
            prime[priLen++] = i;
        for(j = i*i; j <=M; j+=i)
            vis[j] = 1;
    }
}

文艺的素数筛选法

保证了每个数只被筛选一次，且是被最小质因数筛选，时间复杂度达 到o（n）。

const int N=1000000;
int pr[N],p[N/10],lp;

void gp()
{
    for(int i=2;i<N;i++)
    {
        if(!pr[i])
		p[lp++]=pr[i]=i;
        for(int j=0;j<lp && i*p[j]<N;j++)
	{
            pr[i*p[j]]=p[j];
            if(i%p[j]==0)
		break;
        }
    }
}