筛法求素数&线性筛法求素数

2021年更新版：
筛法求素数&线性筛法求素数

要理解筛法求素数首先要知道一个定理，
整数唯一分解定理：任意大于等于2的正整数都有且只有一种方式写出其质因子的乘积表达式。
A ＝p₁p₂p₃p₄…p_n(p_i是素数且p_i<=p_j)
eg:2=2
4=22
12=223
36=2233
也就是说任意一个合数都能分成一个素数和另一个数（另一个数可以是素数也可以是合数）的乘积
eg:4=22
12=26
15=3*5
即任意一个合数都有一个比它本身小的素因子。
这样的话，如果要求1－n的所有素数，只要把1－n内所有素数的倍数删掉，剩下的就都是素数了。
先看代码：
普通筛法求素数：

#include<iostream>
#include<cstring>
#include<cmath>
using namespace std;
#define ll long long
#define N 100
bool isprime[N];／／1
int pri[N];／／2
void prime()
{
    int cnt=1;
    memset(pri,0,sizeof(pri));
    isprime[1]=1;
    for(int i=2;i<=sqrt(N);i++)／／3
    {
        if(!isprime[i])
        {
            pri[cnt++]=i;／／4
            for(int k=i*i;k<=N;k+=i)／／5
            {
                isprime[k]=1;／／表示k是一个合数
            }
        }
    }
    return;
}

上边代码里：

1. isprime数组的作用是，判断一个数是否是素数，我们可以规定
   if(isprime[i]==0)代表i是素数
   else 代表i是合数
   一开始把isprime数组的值全部置0，假设1－n内所有数都是素数，然后把ispime［1］特判掉

2.pri数组用来存储已经确定的素数

3.一个循环从2开始到sqrt（n）结束，循环内如果发现isprime［i］＝＝0，那么把i存到pri数组内（即上面的4那一行），接着再来一个循环把i这个素数的倍数全部删掉（上面的5），这个循环从i＊i开始，每次加i，到n结束，这样就可以把i＊i，i＊(i+1)，i＊(i+2)…i＊(i+k)，把n以内i的倍数全部删掉。

接着说1 .为什么到sqrt(n)就可以结束了：
假设一个小于n的合数 x=ab ,设 a<=b;
那么a必定小于sqrt(n)，也就是说这个合数必定有一个小于 sqrt(n)的素因子，那么这个数一定能在2－sqrt（n）内被某个数删掉
2.为什么从i ＊i开始而不是从2＊i开始
假设一个合数x=ij （j < i ）
那么这个数在遍历 j 时就已经被删掉了

然后看线性筛法
仔细分析上面的筛法其实是会造成重复筛合数，比如 12 在遍历2的时候筛到2＊6 筛掉一次，在遍历3的时候筛到3＊4 又筛了一次
快速筛法就解决了这个问题
先看代码：

#include<iostream>
#include<cstring>
#include<cmath>
using namespace std;
#define ll long long
#define N 100
bool isprime[N];
int pri[N];
void prime2()
{
    int cnt=1;
    memset(pri,0,sizeof(pri));
    isprime[1]=1;
    for(int i=2;i<=N;i++)
    {
        if(!isprime[i])
        {
            pri[cnt++]=i;
        }
        
        for(int j=1;j<cnt&&i*pri[j]<=N;j++)／／1
        {
            isprime[i*pri[j]]=1;
            if(i%pri[j]==0)／／2
                break;
        }
    }
    return;
}

上边代码：
1。不管i是素数还是合数，都会执行到“1“这个地方。
这个循环的作用是

1. 如果i是素数，会筛掉所有的i＊k（其中k是所有小于等于i的素数）
   看代码j从1开始，到cnt-1结束，每次把i＊pri［j］标记为合数，即是说把从pri［1］＝2到pri[cnt-1]，与i的乘积筛掉
   eg：i＝5 时，j从1开始，此时cnt等于3 ，即j到3结束，会把5＊2，5＊3，筛掉
2. 如果i是合数，会筛掉所有的i＊k（其中k是所有小于等于i的最小素因子的素数）
   看代码j从1开始，到cnt-1结束，每次把i＊pri［j］标记为合数，然后“／／2“这一行的作用就是在pri［j］等于i的最小素因子时跳出循环。因为pri［j］是从小到大枚举的，所以遇见的第一个i％pri［j］＝＝0就是i的最小素因子

现在证明这个筛法

1. 不会造成重复筛
   对于任意一个 x ＝p₁p₂p₃p₄…p_n(p_i是素数且p_i<=p_j)
   这个x能且只能被y＝p₂p₃p₄…p_n(p_i是素数且p_i<=p_j)筛掉。
   eg：42＝2＊3＊7，42只能被3＊7＝21筛掉。
2. 可以把所有合数都筛掉
   对于任意一个 x ＝p₁p₂p₃p₄…p_n(p_i是素数且p_i<=p_j)
   这个数一定能被y＝p₂p₃p₄…p_n(p_i是素数且p_i<=p_j)筛掉。
   所以一定可以把所有合数都筛掉