质数的判定和应用

一,定义

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

二,判定

1,试除法

预备性质:如果 x x x a a a的约数,那么 x / a x/a x/a 也是 x x x的约数;已知 n ∣ d n|d nd 即有 n / d ∣ d n/d|d n/dd
枚举所有 1 − s q r t ( n ) 1-sqrt(n) 1sqrt(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;
}

所以

  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

流苏贺风

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值