快速找素数的方法

原理是合数一定可以被分解成两个质数相乘，所以我们只要去掉质数的倍数，那么剩下的一定是质数（假设p最后被留了下来，那么p一定不是2->p-1的倍数，那么p就是质数了）。

• 埃氏筛法(依据倍数)

int n,primes[1000010],cnt;
bool symbol[1000010];
int get_primes(int n){
    for(int i=2;i<=n;i++){
        if(!symbol[i]){//如果没有被标记的话
            primes[cnt++]=i;
            for(int j=i+i;j<=n;j+=i)symbol[j]=true;
        }
    }
    return cnt;
}

• 欧拉线性筛法（依据最小质因子）
  欧拉线性筛法的核心是用每个合数的最小质因子去筛，让每个合数都只被自己的最小质因子筛掉（因为一个数只有一个最小质因子），减少筛的次数。

int n,primes[1000010],cnt;
bool sl[1000010];//记录是否已被除掉
int get_primes(int n){
    for(int i=2;i<=n;i++){
        if(!sl[i])primes[cnt++]=i;//选素数
        for(int j=0;primes[j]<=n/i;j++){//筛数
            sl[primes[j]*i]=true;//让每个数都被其最小质因数筛掉
            if(i%primes[j]==0)break;//是i最小质因子时停止
        }
     }
     return cnt;
}

疑问1

为什么这段代码就说明合数是被自己的最小质因子筛掉的？

  for(int j=0;primes[j]<=n/i;j++){
            sl[primes[j]*i]=true;
            if(i%primes[j]==0)break;
        }

• 首先，因为我们的primes[j]是从小到大枚举的，第一次出现i%primes[j]==0时，说明primes[j]就是i的最小质因子，那么也是primes[j]*i（要筛掉的合数）的最小质因子。
• 如果i%primes[j]!=0时，说明pj小于i的所有质因子。那么pji的最小质因子就是pj（某数i自己的质因子的最小），它也就被pj干掉了。

疑问2

这样子真的能筛掉所有的合数吗？（仅仅靠一个合数的最小质因子来筛）
答案是可以的。

• 假设我们合数x的最小质因子是d，那么当我们i模拟到到x/d时（d *i(x/d)），就会在那一轮被筛掉。合数是任意的，所以所有合数都会被筛掉。

疑问3

所有的合数都会只筛一次吗？
答案是肯定的。我们从循环知道，当数是合数的最小质数时，我们枚举数就会停止。同时，每一个数都只有一个最小质数，说明我们这个数是第一次遇到被标记。第一次遇到被标记就马上停止了，所以每个数都只被筛了一次。这也是为什么被称作线性筛法。