质数筛的Blog

提示：文章写完后，目录可以自动生成，如何生成可参考右边的帮助文档

---

前言

最近忘了质数筛原理了写个blog复习一下

一、如何判定质数？

1.试除法：设置$i$从2->$n$分别去除此数，若不会被整除，则为素数；
2.注意应为比如一个数$15$有$<3,5>,<1,15>$质因数；所以假如，
$$n|i$$ $$则n|n/i$$所以我们只需要验证$$i:i<=n/i(注意！必须有”=”)$$的部分就可以了（代码如下）；

bool is_prime(int x)
{
    if(x<2)return false;
    for(int i=2;i<=x/i;i++)
    {
        if(x%i==0)
            return false;
    }
    return true;
}

二、分解质因数

1.还是根据我们的算数基本定理：

$$\mathbf{m}={\mathbf{P}}_1^{{\mathbf{\alpha }}_1}{\mathbf{P}}_2^{{\mathbf{\alpha }}_2}{\mathbf{P}}_3^{{\mathbf{\alpha }}_3}{\mathbf{P}}_4^{{\mathbf{\alpha }}_4}.....{\mathbf{P}}_{\mathbf{k}}^{{\mathbf{\alpha }}_{\mathbf{k}}}$$
其中$P为质因数$：所以一个数可以拆成许多质因数的K次幂乘积；
如$$12=2^2\times 3^1$$
所以我们也只需要从小到大把能见到的质数全除就可：
比如：
$$12/2=6/2=3/2(\times )-->3/(2+1)=1$$
这样子的计算我们就已经完成了质因数的分解：
既有2个2，和1个3；
代码如下：

#include<iostream>
using namespace std;

void divide(int x)
{
    for(int i=2;i<=x/i;i++)
    {
        if(x%i==0)
        {
            int cnt=0;
            while(x%i==0)
            {
                x/=i;
                cnt++;
            }
        cout<<i<<" "<<cnt<<endl;
            
        }
    }
    if(x>1)
    cout<<x<<" "<<1<<endl;
}

三、质数筛

1：朴素筛法的原理：

比如$24$;
我们一个for循环从小到大穷举质数数，碰到质数就把他的倍数标记上；

$i$	标记
2	$4,6,8,10,12,14,16,,18,20,22,24$
3	$6,9,12,15,18,21,24$
4	‘4’已被标记跳过
5	$10,15,20$
6	‘6’已被标记跳过
7	…
…

这样子我们就可以把1~n素数找出来(标记的倍数都被排除了，剩下的就都是最初的质数)不信你看表格$i$列的4，6是不是因为之前的标记给排除了？

2：线性筛法

不难发现我们上面的做法有好多重复操作！
我们在观察一下表格

$i$	标记
2	$4,$$6$$,8,10,$$12$$,14,16,,18,20,22,24$
3	$6$ $,9,$$12$$,15,18,21,24$
4	‘4’已被标记跳过
5	$10,15,20$
6	‘6’已被标记跳过
7	…
…

上面的6，12，24…早在第一行就已经被标记过没必要再标记一次，少标记一次就减少一次循环那么每一个数都有且仅有一次被访问；

3.那我们重新定义一个筛法规则：

$i$是质数就筛掉$i$与$primes$数组中所有质数的乘积；i是非质数就筛掉$i$与$primes$中所有<=最小质因子的乘积，并且若if (i % primes[j] == 0)break;
primes是存质数的数组

那我们依照这个规则再列一遍：24

$i$	标记	primes
2	4	2
3	$6,9$	2,3
4	被标记但不跳过$8$，但4%primes[0]=0;break;	2,3
5	$10，15$	2,3,5
6	被标记但不跳过$12$，但6%primes[0]=0;break;	2,3,5
7	$14,21$	2,3,5,7
8	被标记 但不跳过$16$	2,3,5,7
9	被标记 但不跳过 $18,27>24$;但9%primes[1]=0;break;
10	被标记 但不跳过 $20$;但10%primes[0]=0;break
…
11	$22$	2,3,5,7,11
12	被标记 但不跳过但是if (i % primes[j] == 0)break;	2,3,5,7,11
13	26>24
…
17	34>24	2,3,5,7,11,13,17
…
19	38>24	2,3,5,7,11,13,17,19
…
23

如上通过这个方法，每个数都只筛了一遍，大减少了循环次数；
时间复杂度是O($n$)的！！！

4.完整代码如下：

#include<iostream>

using namespace std;

const int N=1e6+6;

int primes[N];bool st[N];
int cnt;

void get_primes(int x)
{
    for(int i=2;i<=x;i++)
    {
        if(!st[i])primes[cnt++]=i;
        for(int j=0;primes[j]<=x/i;j++)
        {
            st[primes[j]*i]=true;
            if(i%primes[j]==0)break;
        }
    }
}
int main()
{
    int n;
    cin>>n;
    get_primes(n);
    cout<<cnt<<endl;
}

一键三连