#埃氏筛法
原理:如果一个数是质数,那么他的倍数必是合数,所以我们在碰到质数的时候把他的倍数全都标记为合数,那么没有被标记到的就是质数
下面看代码(判定n以内的素数):
const int N=1e7+5;
bool ok[N];
int prim[N];
int pos=0;
void p(int n)
{
memset(ok,true,sizeof(ok));
ok[0]=ok[1]=false;
for(int i=2;i<=n;i++)
{
if(ok[i])
{
prim[pos++]=i;
for(int j=2*i;j<=n;j+=i)
{
ok[j]=false;
}
}
}
}
注意到:这里的“倍数”是从2 * i开始的。其实从i * i开始就行了,因为小于i * i的素数已经被[2,i-1]区间的某个数标记过了,比如(i-2)*i已经被i-2标记过,所以不需要用i再去跑一遍。
下面是优化代码(其实就是把2 * i 换成了i * i,这样更快一些):
const int N=1e7+5;
bool ok[N];
int prim[N];
int pos=0;
void p(int n)
{
memset(ok,true,sizeof(ok));
ok[0]=ok[1]=false;
for(int i=2;i<=n;i++)
{
if(ok[i])
{
prim[pos++]=i;
for(int j=i*i;j<=n;j+=i)
{
ok[j]=false;
}
}
}
}