【转载】 判断质数/素数——我知道的最快的方法

转载自:https://blog.csdn.net/songyunli1111/article/details/78690447

标准版:大部分人都知道的比较快的方法:判断从 2 到 sqrt(n) 是否存在其约数,时间复杂度 O(sqrt(n))

高配版:判断 2 之后,就可以判断从 3 到 sqrt(n) 之间的奇数了,无需再判断之间的偶数,时间复杂度 O(sqrt(n)/2)

尊享版:

首先看一个关于质数分布的规律:大于等于 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 的相邻两侧。因此在 5 到 sqrt(n) 中每 6 个数只判断 2 个,时间复杂度 O(sqrt(n)/3)。

(有同学对这里的原理不理解的话,可以去参考我的另一篇原创《我理解的素数——约数求法》)

在高配版和尊享版中,都是一个剪枝的思想,高配版中裁剪了不必要的偶数,尊享版中裁剪了不和 6 的倍数相邻的数,虽然都没有降低时间复杂度的阶数,但都一定程度上加快了判断的速度。

在此给出尊享版 C++ 代码:

#include <iostream>  
#include <math.h>  
using namespace std;  
int isPrime(int n)  
{   //返回1表示判断为质数,0为非质数,在此没有进行输入异常检测  
    float n_sqrt;  
    if(n==2 || n==3) return 1;  
    if(n%6!=1 && n%6!=5) return 0;  
    n_sqrt=floor(sqrt((float)n));  
    for(int i=5;i<=n_sqrt;i+=6)  
    {  
        if(n%(i)==0 | n%(i+2)==0) return 0;  
    }
    return 1;  
}   
int main()  
{  
    int flag;  
    flag=isPrime(101);  
    cout<<flag<<endl;  
    return 0;  
}  
  • 2
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值