素数筛法详解（欧拉筛&埃氏筛）

素数的判断有一个最为朴素的算法思想就是从定义出发的试除法，也就是一些教科书上的版本，这里就不再介绍。

这里介绍的是 埃拉托斯特尼筛法 和 欧拉筛法，欧拉筛法可以说是埃氏筛法的升级版，因为少了一些重复的操作而更加高效，先来看看埃氏筛法。

埃氏筛法的思想也比较好理解，就是判断当前位置的数是否为素数，如果是，则以该数为因子的所有数均标记为合数（筛除掉 ），判断到√n的位置即可（稍加思考即可），n为数组总长度。最终剩下的未筛除的数即为素数。

埃拉托斯特尼筛法代码：

#include<iostream>
#include<cstdio>
#include<vector>
#include<cstring>
using namespace std;

const int maxn =1e7+10;
vector<int> primes;//存放素数的不定长数组
bool judge[maxn]; //筛除判断

int main()
{
    int i;
    memset(judge,true,sizeof(judge));
    judge[0]=false;
    judge[1]=false;     //0和1都不是素数
    for(i=2;i<maxn;i++)
    {
        if(i*i>maxn)   //i判断到√n即可
            break;
        if(judge[i])
        {
            primes.push_back(i);
            for(long long j=i*i;j<maxn;j+=i) //j从i²开始判断更高效，因为之前2*i，3*i等等已经被筛除了，j要用long long型来防止溢出
            {
                judge[j]=false;
            }
        }
    }
    vector<int>::iterator it;
    for (it = primes.begin(); it != primes.end(); it++)
    cout << *it << endl;
    return 0;

}

现在介绍第二种也是笔者已知的最高效的筛法：欧拉筛（线性筛）

欧拉筛思想的核心是要保证的是每个合数只被这个合数最小的质因子筛除，而且只筛一次，没有重复筛除，这里就需要用到这行神秘代码if(i%prime[j]==0) break;。

欧拉筛的难点就在于对if (i % prime[j] == 0)这步的理解，当i是prime[j]的整数倍时，记 m = i / prime[j]，那么 i * prime[j+1] 就可以变为 (m * prime[j+1]) * prime[j]，这说明 i * prime[j+1] 是 prime[j] 的整数倍，不需要现在筛出，因为在之后筛除过程中i * prime[j+1] 这个合数一定会被prime[j]筛除，prime[j]之后的所有素数同理，所以break跳出循环。

这样就保证了每个合数仅被它最小的质因子筛除一次，提高效率。（如果不理解可以多读几遍哈）

欧拉筛代码：

#include<iostream>
#include<cstdio>
#include<vector>
#include<cstring>
using namespace std;

const int maxn = 1e6+10;
int cnt=0; //记录已知素数数量
vector<int> primes; //存放素数的不定长数组
bool judge[maxn]; //筛除标记

int main()
{
    int i,j;
    memset(judge,true,sizeof(judge));
    judge[0]=false;
    judge[1]=false;// 0和1都不是素数
    for(i=2;i<maxn;i++)
    {
        if(judge[i])
        {
            primes.push_back(i);
            cnt++;
        }
        for(j=0;j<cnt && i*primes[j]<=maxn;j++)
        {
            judge[i*primes[j]]=false; //筛除
            if(i%primes[j]==0) //关键代码
                break;
        }
    }
    vector<int>::iterator it;
    for(it=primes.begin();it!=primes.end();it++)//遍历
    {
        printf("%d\n",*it);
    }
    return 0;
}

经过测验，在10^7数量级数据下，埃氏筛法平均需要177ms，而欧拉筛平均需要156ms。

好像也没差多少，掩面哭泣    {{{(>_<)}}}，可能数据还不够大？

以上即为个人对于素数筛法的理解，如有不足，请多多指教！