0.引入

那年,机房里来了个新教练, 口胡鼻祖lhy

第一节课,带我们体验了暴力的神奇,

第二节课,带我们体验了随机数的玄妙,

……

那节课,便是我第一次接触到Miller Rabbin算法,

直到现在,终于搞懂了一些。


该算法(名字过长,不想打了)主要是解决快速判断一个极大的数是否是质数的问题。

我们知道,能保证正确的最快的算法,就是的复杂度,不能再小了,对于一个很大的long long,复杂度达到O(1e9)

但是该算法却能在的时间复杂度内判断,(如果用了光速乘,还可以变为

那究竟是为什么呢?

我们先特判了小质数后,再进入下面的算法

1.费马定理

我们知道,对于一个质数p,和比它小的数a,有

所以我们可以先判断 p 是否符合费马定理,不符合就肯定不是质数。

但是符合也不一定是质数,先不说随便在小于一个合数 p 的数中选一些数可能恰好都符合费马定理,就算把1~p-1都穷举完了,也有合数是可以满足上面的式子的,那就是卡迈克尔数。

但是这一步判断至少确定了枚举到的 a 是有逆元的吧……

2.二次探测

若有

其中 p 是质数,y < p,则

 

这就是二次探测的原理

然后大家知道,这篇文章不是在证定理,而是在讲算法,所以注重算法过程而不是算法推导,所以笔者就直接讲过程了!

由于上面费马定理已经判断了

而若p是质数,p-1就一定是偶数了(2已经特判掉了)可以分解出很多2出来,

那么先算出 ,然后每次  把它平方,只要第一次碰到 a==1 的情况,就看前一次的 a 是多少,如果不是p-1那么该数就是合数,否则就继续算法过程,总会碰到 a==1 的情况,因为最后会乘到 

3.随机数

随便选一个数 a 来做二次探测,让合数躲过的几率只有10%左右,

那么随机选10个、20个、30个数做二次探测,让合数躲过的几率就小得多了,可以小到同出现相同指纹的概率这么小,(但前提是足够随机)

所以我们就放心做十次二次探测,每次选一个1~p-1中的随机 a 跑上面的过程,把这个10算作一个log的话,复杂度也只有了。

4.例题

FZU1649 质数还是合数

一道该算法的板题

CODE

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
#include<ctime>
#include<queue>
#include<stack>
#include<map>
#include<vector>
#include<algorithm>
#define MAXN 65545
#define MAXM 35
#define ENDL putchar('\n')
#define LL long long
#define DB double
#define lowbit(x) ((-x)&(x))
//#define int LL
using namespace std;
inline LL read() {
LL f = 1,x = 0;char s = getchar();
while(s < '0' || s > '9') {if(s == '-')f = -1;s = getchar();}
while(s >= '0' && s <= '9') {x = x * 10 + (s - '0');s = getchar();}
return x * f;
}
const LL jzm = 1000000007;
int n,m,i,j,s,o,k;
LL safemul(LL a,LL b,LL zxy) {
LL res = 0;if(a < b) swap(a,b);
while(b>0) {
if(b & 1ll) (res += a) %= zxy;
(a <<= 1ll) %= zxy;
b >>= 1;
}
return res;
}
LL qkpow(LL a,LL b,LL zxy) {
LL res = 1;
while(b>0) {
if(b & 1ll) res = safemul(res,a,zxy);
a = safemul(a,a,zxy);
b >>= 1;
}
return res % zxy;
}
LL Rand() {return rand()*rand()+rand();}
bool Miller_rabbin(LL p) {
if(p == 2 || p == 3 || p == 5) return 1;
if(p <= 1) return 0;
LL nm = p-1;
nm /= lowbit(nm);
for(int id = 1;id <= 10;id ++) {
LL a = Rand() % p;
if(qkpow(a,p-1,p) != 1ll) return 0;
a = qkpow(a,nm,p);
LL pre = p-1;
while(a != 1) {
pre = a;
a = safemul(a,a,p);
}
if(pre != p-1) return 0;
}
return 1;
}
int main() {
srand(time(0));
LL N;
while(scanf("%lld",&N) == 1) {
printf(Miller_rabbin(N) ? "It is a prime number.\n":"It is not a prime number.\n");
}
return 0;
}

Miller Rabbin 算法—费马定理+二次探测+随机数 (讲解+例题:FZU1649 Prime number or not)的更多相关文章

  1. 【数论基础】素数判定和Miller Rabin算法

    判断正整数p是否是素数 方法一 朴素的判定   

  2. Miller Rabin算法详解

    何为Miller Rabin算法 首先看一下度娘的解释(如果你懒得读直接跳过就可以反正也没啥乱用:joy:) Miller-Rabin算法是目前主流的基于概率的素数测试算法,在构建密码安全体系中占有重 ...

  3. Miller-Robin与二次探测

    素数在数论中经常被用到.也是数论的基础之一. 人们一直在讨论的问题是,怎样快速找到素数?或者判断一个数是素数? 1.根号n枚举 原始暴力方法. 2.埃氏筛 每个合数会被筛质因子次数次.复杂度O(Nlo ...

  4. Miller Rabin 算法简介

    0.1 一些闲话 最近一次更新是在2019年11月12日.之前的文章有很多问题:当我把我的代码交到LOJ上,发现只有60多分.我调了一个晚上,尝试用{2, 3, 5, 7, 11, 13, 17, 1 ...

  5. Miller Rabin算法学习笔记

    定义: Miller Rabin算法是一个随机化素数测试算法,作用是判断一个数是否是素数,且只要你脸不黑以及常数不要巨大一般来讲都比\(O(\sqrt n)\)的朴素做法更快. 定理: Miller ...

  6. 算法(二)之遗传算法(SGA)

    算法(二)之遗传算法(SGA) 遗传算法(Genetic Algorithm)又叫基因进化算法或进化算法,是模拟达尔文的遗传选择和自然淘汰的生物进化过程的计算模型,属于启发式搜索算法一种. 下面通过下 ...

  7. (Miller Rabin算法)判断一个数是否为素数

    1.约定 x%y为x取模y,即x除以y所得的余数,当x<y时,x%y=x,所有取模的运算对象都为整数. x^y表示x的y次方.乘方运算的优先级高于乘除和取模,加减的优先级最低. 见到x^y/z这 ...

  8. java 解决Hash(散列)冲突的四种方法--开放定址法(线性探测,二次探测,伪随机探测)、链地址法、再哈希、建立公共溢出区

    java 解决Hash(散列)冲突的四种方法--开放定址法(线性探测,二次探测,伪随机探测).链地址法.再哈希.建立公共溢出区 标签: hashmaphashmap冲突解决冲突的方法冲突 2016-0 ...

  9. 70 数组的Kmin算法和二叉搜索树的Kmin算法对比

    [本文链接] http://www.cnblogs.com/hellogiser/p/kmin-of-array-vs-kmin-of-bst.html [分析] 数组的Kmin算法和二叉搜索树的Km ...

随机推荐

  1. JavaScript之parseInt()方法

    parseInt(string, radix):用于解析一个字符串并返回指定基数的十进制整数或者NaN string参数为被解析的值,如果该值不是一个字符串,则会隐式的使用toString()方法转化 ...

  2. vue虚拟dom和diff算法

    vue的虚拟dom和diff算法 1.虚拟dom 虚拟dom,我的理解就是通过js对象的方式来具体化每一个节点,把dom树上面的每个节点都变为对象里的一个元素,元素的子元素变为子节点,节点上面的cla ...

  3. Quartus II 13.0 sp1的官方下载页面

    今天为了下个ModelSim跑到网上去找下载资源,清一色的百度网盘,下载速度60k/s,简直有病,于是跑到Intel官网上把连接挖出来了,供各位直接下载 实测使用IDM多线程下载速度可以轻松上到数MB ...

  4. sublime_text 3安装Emmet时出现PyV8警告

    使用Emmet是需要在PyV8依赖下才可以的.1. 下面是下载网址:PyV8下载地址 下载自己系统版本的压缩包,然后解压,自己创建一个名为PyV8文件夹.将解压后的文件放入该文件夹里. 打开首选项里的 ...

  5. 物联网?快来看 Arduino 上云啦

    作者:HelloGitHub-Anthony 这里是 HelloGitHub 推出的讲解开源硬件开发平台 Arduino 的系列教程. 第一篇:Arduino 介绍和开发环境搭建 第二篇:制作温湿度显 ...

  6. VisonPro · 视觉定位工具包示例

    一.概述 视觉定位工具包一般包含: 1.相机取像: 2.图像九点标定: 3.Mark点粗定位: 4.建立粗定位坐标系: 5.Mark点精定位 6.输出Mark点坐标,角度等信息. 二.分类 1.单特征 ...

  7. oracle备份数据库数据及导入数据库

    1.oracle数据库备份和导入 bat 脚本 scott oracle数据库用户名称 123456 数据库scott用户下的密码 192.168.124.8 本电脑IP orcl 为oracle库 ...

  8. CesiumJS 2022^ 源码解读[7] - 3DTiles 的请求、加载处理流程解析

    目录 1. 3DTiles 数据集的类型 2. 创建瓦片树 2.1. 请求入口文件 2.2. 创建树结构 2.3. 瓦片缓存机制带来的能力 3. 瓦片树的遍历更新 3.1. 三个大步骤 3.2. 遍历 ...

  9. k8s之有状态服务部署基石(基础知识)

    PV&PVC&HeadlessService 4.1.什么是无状态/有状态服务? 无状态服务: 1.没有实时的数据需要存储 (即使有,也是静态数据) 2.服务集群网络中,拿掉一个服务后 ...

  10. S32K148-CAN收发

    最近在搞一个转换板开发,大概意思把CAN信号转换成SPI信号,方案有两种:1)通过硬件电路直接把信号的bit位一位一位移给两个集成芯片:2)通过MCU接收CAN信号,再把信号变量转换成SPI信号发送给 ...