质数的判定和应用

一，定义

• 质数：>1 的自然数中，仅有1和本身的为因数的数，就是质数

二，判定

1，试除法

预备性质：如果$x$是$a$的约数，那么$x/a$ 也是$x$的约数；已知$n|d$ 即有$n/d|d$
枚举所有$1-sqrt(n)$的数字，判断可否被整除即可

2，线性筛

• 用每个数的最小质因子筛去和数
• 算数基本定理的应用

void prime()
{
    for(int i=2;i<=n;i++)
    {
        if(!st[i])pri[++cnt]=i;
        for(int j=1;j<=cnt;j++)
        {
            st[pri[j]*i]=1;
            if(i%pri[j]==0)break;
        }
    }
}

参考：判断一个数是否为质数/素数——从普通判断算法到高效判断算法思路

3，六单元快进判法（最快）

• 首先看一个关于质数分布的规律：大于等于5的质数一定和6的倍数相邻。例如5和7，11和13,17和19等等；

证明：令x≥1，将大于等于5的自然数表示如下：
······ 6x-1，6x，6x+1，6x+2，6x+3，6x+4，6x+5，6(x+1），6(x+1)+1 ······
可以看到，不在6的倍数两侧，即6x两侧的数为6x+2，6x+3，6x+4，由于2(3x+1)，3(2x+1)，2(3x+2)，所以它们一定不是素数，再除去6x本身，显然，素数要出现只可能出现在6x的相邻两侧。这里有个题外话，关于孪生素数，有兴趣的道友可以再另行了解一下，由于与我们主题无关，暂且跳过。这里要注意的一点是，在6的倍数相邻两侧并不是一定就是质数。

此时判断质数可以6个为单元快进，即将方法（2）循环中i++的一步步长加大为六步，加快判断速度。
原因是：假如要判定的数为n，则n必定是6x-1或6x+1的形式，对于循环中6i-1，6i，6i+1,6i+2，6i+3，6i+4，其中如果n能被 6i，6i+2，6i+4整除，则n至少得是一个偶数，但是6x-1或6x+1的形式明显是一个奇数，故不成立；另外，如果n能被6i+3整除，则n至少能被3整除，但是6x能被3整除，故6x-1或6x+1（即n）不可能被3整除，故不成立。综上，循环中只需要考虑6i-1和6i+1的情况，即循环的步长可以定为6，每次判断循环变量k和k+2的情况即可

bool su(int a){
    if(a==1) return 0;
    if(a==2||a==3) return 1;
    if(a%6!=1&&a%6!=5) return 0;
    int temp=sqrt(a);
    for(int i=5;i<=temp;i+=6) if(a%i==0||a%(i+2)==0) return 0;
    return 1;
}

所以