素性测试是数论题中比较常用的一个技巧。它可以很基础,也可以很高级(哲学)。这次主要要介绍一下有关素数判断的奇技淫巧

素数的判断主要分为两种:范围筛选型&&单个判断型

我们先从范围筛选型这种常用的开始讲起,这里采用模板题Luogu P3383 【模板】线性筛素数来进行测试

1.埃氏筛

这是最常用的筛法了,思路也很简单:任何一个素数的倍数都是合数

然后我们O(n)扫一遍,同时筛去素数的倍数

但是有一些数如6,会被2和3都筛去一次,就造成了效率上的浪费,所以复杂度经证明为**O(n log log n)

CODE

#include<cstdio>
using namespace std;
const int N=10000005;
bool vis[N];
int n,m,x;
inline char tc(void)
{
static char fl[100000],*A=fl,*B=fl;
return A==B&&(B=(A=fl)+fread(fl,1,100000,stdin),A==B)?EOF:*A++;
}
inline void read(int &x)
{
x=0; char ch=tc();
while (ch<'0'||ch>'9') ch=tc();
while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=tc();
}
inline void get_prime(int m)
{
register int i,j;
for (vis[1]=1,i=2;i<=m;++i)
if (!vis[i]) for (j=i<<1;j<=m;j+=i) vis[j]=1;
}
int main()
{
//freopen("CODE.in","r",stdin); freopen("CODE.out","w",stdout);
read(n); read(m); get_prime(n);
while (m--)
{
read(x);
puts(vis[x]?"No":"Yes");
}
return 0;
}

2.线性筛(欧拉筛)

这其实是对上者的优化,我们意识到一个数应该只有它的最小质因数删去,所以我们可以一边筛数的同时一边记录素数,这就是真正的O(n)复杂度

CODE

#include<cstdio>
using namespace std;
const int N=10000005;
int prime[N],n,m,x,cnt;
bool vis[N];
inline char tc(void)
{
static char fl[100000],*A=fl,*B=fl;
return A==B&&(B=(A=fl)+fread(fl,1,100000,stdin),A==B)?EOF:*A++;
}
inline void read(int &x)
{
x=0; char ch=tc();
while (ch<'0'||ch>'9') ch=tc();
while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=tc();
}
inline void Euler(int n)
{
register int i,j;
for (vis[1]=1,i=2;i<=n;++i)
{
if (!vis[i]) prime[++cnt]=i;
for (j=1;j<=cnt&&i*prime[j]<=n;++j)
{
vis[i*prime[j]]=1;
if (!(i%prime[j])) break;
}
}
}
int main()
{
//freopen("CODE.in","r",stdin); freopen("CODE.out","w",stdout);
read(n); read(m);
Euler(n);
while (m--)
{
read(x);
puts(vis[x]?"No":"Yes");
}
return 0;
}

注意上面的那句话:

if (!(i%prime[j])) break;

这保证了线性筛的效率,不会产生重复,因为当i%prime[j]==0时这个数就是让后面的数删去。

3.基础素性测试

这是最基本的素数判定法了吧。从2到sqrt(x)枚举是否有数能够整除x

证明的话很简单,因为如果这个数是素数,那么它的因数必定为1和x,若其因数大于sqrt(x),那么平方后就大于x,这显然不可能。

所以我们O(sqrt(x))判断一次

CODE

#include<cstdio>
#include<cmath>
using namespace std;
int n,m,x;
inline char tc(void)
{
static char fl[100000],*A=fl,*B=fl;
return A==B&&(B=(A=fl)+fread(fl,1,100000,stdin),A==B)?EOF:*A++;
}
inline void read(int &x)
{
x=0; char ch=tc();
while (ch<'0'||ch>'9') ch=tc();
while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=tc();
}
inline bool check(int x)
{
if (!(x^1)) return 0;
register int i; int bound=(int)sqrt(x);
for (i=2;i<=bound;++i)
if (!(x%i)) return 0;
return 1;
}
int main()
{
//freopen("CODE.in","r",stdin); freopen("CODE.out","w",stdout);
read(n); read(m);
while (m--)
{
read(x);
puts(check(x)?"Yes":"No");
}
return 0;
}

4.对于算法3的优化

首先我们看一个结论:

大于等于5的质数一定和6的倍数相邻。

证明等参考:dalao's blog

然后同3,我们只不过每次快进6个单位,然后常数就得到了难以言喻都优化(一跃成为此题最快的算法)

CODE

#include<cstdio>
#include<cmath>
using namespace std;
int n,m,x;
inline char tc(void)
{
static char fl[100000],*A=fl,*B=fl;
return A==B&&(B=(A=fl)+fread(fl,1,100000,stdin),A==B)?EOF:*A++;
}
inline void read(int &x)
{
x=0; char ch=tc();
while (ch<'0'||ch>'9') ch=tc();
while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=tc();
}
inline bool check(int x)
{
if (!(x^1)) return 0;
if (!(x^2)||!(x^3)) return 1;
if ((x%6)^1&&(x%6)^5) return 0;
register int i; int bound=(int)sqrt(x);
for (i=5;i<=bound;i+=6)
if (!(x%i)||!(x%(i+2))) return 0;
return 1;
}
int main()
{
//freopen("CODE.in","r",stdin); freopen("CODE.out","w",stdout);
read(n); read(m);
while (m--)
{
read(x);
puts(check(x)?"Yes":"No");
}
return 0;
}

5.Miller-Rabin算法

这是历史上判断素数最快的方法了吧(但在此题中被算法4吊打了)

首先,这个算法基于费马小定理和二次探测定理:

二次探测定理:如果p是奇素数,则 x2≡1(modp)的解为x = 1或x = p - 1(mod p)

所以我们可以把x变成r*2^t的形式,其中r是一个奇数

然后我们结合两种算法&&快速幂就可以稳定O(log x)进行单次判断了

但是这个算法是一个非完美算法,它每一次都25%的概率是错的,所以我们可以多选择几个数多弄几次

但是偶然在网上看到一段话:

对于大数的素性判断,目前Miller-Rabin算法应用最广泛。一般底数仍然是随机选取,但当待测数不太大时,选择测试底数就有一些技巧了。比如,如果被测数小于4 759 123 141,那么只需要测试三个底数2, 7和61就足够了。当然,你测试的越多,正确的范围肯定也越大。如果你每次都用前7个素数(2, 3, 5, 7, 11, 13和17)进行测试,所有不超过341 550 071 728 320的数都是正确的。如果选用2, 3, 7, 61和24251作为底数,那么10^16内唯一的强伪素数为46 856 248 255 981。这样的一些结论使得Miller-Rabin算法在OI中非常实用。通常认为,Miller-Rabin素性测试的正确率可以令人接受,随机选取k个底数进行测试算法的失误率大概为4^(-k)。

所以对于这一题n=10000000的范围就只需要选择2,7,61即可

CODE

#include<cstdio>
using namespace std;
typedef long long LL;
const int prime[3]={2,7,61};
int n,m,x;
inline char tc(void)
{
static char fl[100000],*A=fl,*B=fl;
return A==B&&(B=(A=fl)+fread(fl,1,100000,stdin),A==B)?EOF:*A++;
}
inline void read(int &x)
{
x=0; char ch=tc();
while (ch<'0'||ch>'9') ch=tc();
while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=tc();
}
inline int quick_pow(int x,int p,int mod)
{
int tot=1;
while (p)
{
if (p&1) tot=((LL)tot*x)%mod;
x=((LL)x*x)%mod; p>>=1;
}
return tot;
}
inline bool Miller_Rabin(int x)
{
if (!(x^2)) return 1;
if (x<2||!(x&1)) return 0;
int t=0,u=x-1;
while (!(u&1)) ++t,u>>=1;
for (register int i=0;i<3;++i)
{
if (!(x^prime[i])) return 1;
if (!(x%prime[i])) return 0;
int lst=quick_pow(prime[i],u,x);
for (register int j=1;j<=t;++j)
{
int now=((LL)lst*lst)%x;
if (!(now^1)&&lst^1&&lst^(x-1)) return 0; lst=now;
}
if (lst^1) return 0;
}
return 1;
}
int main()
{
//freopen("CODE.in","r",stdin); freopen("CODE.out","w",stdout);
read(n); read(m);
while (m--)
{
read(x);
puts(Miller_Rabin(x)?"Yes":"No");
}
return 0;
}

最后给出5个算法的运行结果(无O2)

  1. 埃氏筛

  2. 线性筛(欧拉筛)

  3. 基础素性测试

  4. 优化的算法3

  5. Miller-Rabin

有关素数判断的一些算法(总结&&对比)的更多相关文章

  1. POJ 1811 大素数判断

    数据范围很大,用米勒罗宾测试和Pollard_Rho法可以分解大数. 模板在代码中 O.O #include <iostream> #include <cstdio> #inc ...

  2. #C++初学记录(素数判断2)

    素数判断2 比较简单的算法,没有技术含量 A prime number is a natural number which has exactly two distinct natural numbe ...

  3. 数学:随机素数测试(Miller_Rabin算法)和求整数素因子(Pollard_rho算法)

    POJ1811 给一个大数,判断是否是素数,如果不是素数,打印出它的最小质因数 随机素数测试(Miller_Rabin算法) 求整数素因子(Pollard_rho算法) 科技题 #include< ...

  4. C语言 · 素数判断

     算法提高 素数判断   时间限制:1.0s   内存限制:512.0MB      编写一函数IsPrime,判断某个大于2的正整数是否为素数. 样例输入: 5样例输出:yes 样例输入: 9样例输 ...

  5. 实现100以内的素数输出(Python与C++对比)

    今天从链接http://www.2cto.com/kf/201302/187699.html中看到了Python实现100以内的素数输出的算法,颇受感触.尤其是被其中的Python的列表生成器的使用方 ...

  6. POJ3641 Pseudoprime numbers(快速幂+素数判断)

    POJ3641 Pseudoprime numbers p是Pseudoprime numbers的条件: p是合数,(p^a)%p=a;所以首先要进行素数判断,再快速幂. 此题是大白P122 Car ...

  7. JAVA语言的素数判断,随机数,函数调用

    近来刚学JAVA,就从JAVA写起吧,JAVA判别素数,其实方法和C/C++没什么区别,主要就是想谈一下,其中包括的3个点. (1)JAVA语言产生随机数,random函数,定义参数max的作用是给出 ...

  8. #C++初学记录(素数判断)

    练习题目二 素数判断 A prime number is a natural number which has exactly two distinct natural number divisors ...

  9. 2019-11-29-C#-字典-Dictionary-的-TryGetValue-与先判断-ContainsKey-然后-Get-的性能对比

    原文:2019-11-29-C#-字典-Dictionary-的-TryGetValue-与先判断-ContainsKey-然后-Get-的性能对比 title author date CreateT ...

随机推荐

  1. XSS(跨站脚本攻击)漏洞解决方案

    首先,简单介绍一下XSS定义: 一 . XSS介绍 XSS是跨站脚本攻击(Cross Site Scripting)的缩写.为了和层叠样式表CSS(Cascading Style Sheets)加以区 ...

  2. 洗礼灵魂,修炼python(61)--爬虫篇—【转载】requests模块

    requests 1.简介 Requests 是用Python语言编写的第三方库,所以你需要pip安装,安装过程就略过了.它基于urllib,采用 Apache2 Licensed 开源协议的 HTT ...

  3. rpm安装时出现循环依赖

    在安装git包时提示要安装perl-git,当安装perl-git时又提示要安装git包.报错如下: [root@racdb1 Packages]# rpm -ivh perl-Git-1.7.1-4 ...

  4. centos7执行umount提示:device is busy或者target is busy解决方法

    问题描述: 因为挂载错了,想取消挂载,但是umount报告如下错误: [root@zabbix /]# umount /dev/sdc1 umount: /data1: target is busy. ...

  5. Flask消息闪现

    目录 Flask消息闪现 简单的例子 闪现消息的类别 过滤闪现消息 Message Flashing 参考 Flask消息闪现 一个好的应用和用户界面都需要良好的反馈.如果用户得不到足够的反馈,那么应 ...

  6. C++中的istringstream

    istringstream用于执行C++风格的串流操作. 下面的示例是使用一个字符串初始化istringstream类,然后再使用>>操作符来依次输出字符串中的内容. temp_mon=& ...

  7. 阿里八八β阶段Scrum(1/5)

    今日进度 叶文滔: 修改了α阶段遗留的部分界面BUG,比如状态栏白底等 张岳: 修复用户模块信息修改返回失败的BUG 林炜鸿: 重构了添加事件的代码,增加可修改的特性 黄梅玲: 绘制了新的日程显示格界 ...

  8. Navicat 连接Oracle时提示oracle library is not loaded的问题解决

    笔者使用的Navicat Premium 12启动界面截屏: 请注意是64位的.笔者win7 64位系统. 连接Oracle时提示“oracle library is not loaded”. 解决方 ...

  9. 【转】URL编码(encodeURIComponent和decodeURIComponent)

    转自http://blog.jhonse.com/archives/2032.jhonse 最近在用CI框架的时候,发现一个问题,URL的GET方式链接时,如果用中文字符的话,就会出现问题,提示:链接 ...

  10. Python--Windows下安装虚拟环境

    为什么需要虚拟环境 在python开发中,我们可能会遇到一种情况:就是当前的项目依赖的是某一个版本,但是另一个项目依赖的是另一个版本,这样就会造成依赖冲突.在这种情况之下,我们就需要一个工具能够将这两 ...