埃氏筛法及其改进

#埃氏筛法
原理:如果一个数是质数，那么他的倍数必是合数，所以我们在碰到质数的时候把他的倍数全都标记为合数，那么没有被标记到的就是质数
下面看代码（判定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;
            }
        }
    }
}