Miller_Rabin素数测试算法模板对比

昨天在USACO做了一道判断素数的题,就想着学习一下Miller_Rabin素数测试算法,在网上找到两种模版,第一种十分简洁,运行速度也很快,但是会判错极少的几个非素数;第二种比较麻烦,运行速度很慢,所以我便想找到第一种模版不能判断的非素数特判一下,结果用了一天,电脑只找到10^8以下的,10^9内还有2个没找到,但正确的模版运行速度太慢,我的电脑又太渣,耗不起时间了,姑且先这样,等以后有深入理解有更好的方法再更新一下。

第一种:源自吉林大学ACM模版

刚开始用的是随机数测试,我想到以前了解过只要用2,7,61便可能判断unsigned内的数,便做了修改,交了USACO,由于特判了5,7,11,竟然水过了,然后就愉快地去完了。。。

后来做另一道素数题,终于暴露了它的缺点,但是感觉它判断的很快,而且只有极少的不能判断,便先找出它的特例特判一下,经过一天的无聊枚举终于找到了。

文本对比还是自己写来的快,网上下的用了结果内存爆了而且运行慢,自己写的10s出结果(50mb的文本)。

目前只针对2,3,5,找到的特例是:4097,1048577,16777217,25326001,(960946321是一个10^9内的特例,这个算法枚举一次10^9,我的渣电脑都要8min左右,别说更严密的那个了)

bool witness(long long a,long long n) {
    long long d=1,x;
    int i=ceil(log(n-1.0)/log(2.0))-1;//ceil()是向上取整
    for(;i>=0;i--) {
        x=d;
        d=(d*d)%n;
        if(d==1&&x!=1&&x!=n-1)
            return true;
        if(((n-1)&(1<<i))>0)
            d=(d*a)%n;
    }
    return d!=1;
}

bool miller_rabin(long long n) {
    int s[]={2,3,5};//只用2,3,5和下面的特例判断就能正确判断10^8以内的所有素数,枚举一遍用时40s左右
    if(n==2||n==3||n==5||n==7)
        return true;
    if(n==1||(n&1)==0||0==n%3||0==n%5||n==4097||n==1048577||n==16777217||n==25326001)
        return false;
    for(int i=0;i<3;i++)
        if(witness(s[i], n))
            return false;
    return true;
}

第二种模版:网上大多数用的模版(模版修改自:http://blog.csdn.net/z690933166/article/details/9860937)

就算只用2,7,61判断,也奇慢无比,枚举一次10^8要跑17min左右,但unsigned int内没有特例

typedef unsigned long long LL;

LL modular_multi(LL x,LL y,LL mo) {
	LL t;
	x%=mo;
	for(t=0;y;x=(x<<1)%mo,y>>=1)
		if (y&1)
			t=(t+x)%mo;
	return t;
}

LL modular_exp(LL num,LL t,LL mo) {
	LL ret=1,temp=num%mo;
	for(;t;t>>=1,temp=modular_multi(temp,temp,mo))
		if (t&1)
			ret=modular_multi(ret,temp,mo);
	return ret;
}

bool miller_rabin(LL n) {
	if (n==2||n==7||n==61)
        return true;
	if (n==1||(n&1)==0)
        return false;
	int t=0,num[3]={2,7,61};//2,7,61对unsigned int内的所有数够用了,最小不能判断的数为4 759 123 141;用2,3,7,61在 10^16 内唯一不能判断的数是 46 856 248 225 981
	LL a,x,y,u=n-1;
	while((u&1)==0)
        t++,u>>=1;
	for(int i=0;i<3;i++) {
		a=num[i];
		x=modular_exp(a,u,n);
		for(int j=0;j<t;j++) {
			y=modular_multi(x,x,n);
			if (y==1&&x!=1&&x!=n-1)
				return false;
            //其中用到定理,如果对模n存在1的非平凡平方根,则n是合数。
            //如果一个数x满足方程x^2≡1 (mod n),但x不等于对模n来说1的两个‘平凡’平方根:1或-1,则x是对模n来说1的非平凡平方根
			x=y;
		}
		if (x!=1)//根据费马小定理,若n是素数,有a^(n-1)≡1(mod n).因此n不可能是素数
			return false;
	}
	return true;
}

版权声明:本文为博主原创文章,未经博主允许不得转载。

时间: 2024-11-05 21:58:51

Miller_Rabin素数测试算法模板对比的相关文章

Miller_Rabin 素数测试算法

HDU How many prime numbers Give you a lot of positive integers, just to find out how many prime numbers there are. 根据费马小定理 要求 P 是质数 虽然不是充要条件 但实际上可以根据这个 来测试 素数 注意 a不能是p 的倍数 code: // #include<bits/stdc++.h> using namespace std; #define ll long long lo

Miller-Rabin素数测试算法

由费马小定理可以知道,若p是素数且a是整数,则满足a^p==a(mod p).若存在正整数a不满足a^p==a(mod p),那么n是合数. 定义:令a是一个正整数,若p是合数且满足a^p==a(mod p),则p称为以a为基的伪素数. Miller-Rabin素数测试算法原理: 假如p是素数,且(a,p)==1,(a为任意小于p的正整数),那么a^p-1==1(mod p).如果a^p-1==1(mod p), 则可认为n是素数,取多个底进行试验,次数越多,n为素数概率越大.(我的个人理解多次

优化后的二次测试Miller_Rabin素性测试算法

ll random(ll n) { return (ll)((double)rand()/RAND_MAX*n + 0.5); } ll pow_mod(ll a,ll p,ll n) { if(p == 0) return 1; ll ans = pow_mod(a,p/2,n); ans = ans*ans%n; if(p%2) ans = ans*a%n; return ans; } bool Witness(ll a,ll n) { ll m = n-1; int j = 0; whil

数论 - Miller_Rabin素数测试 + pollard_rho算法分解质因数 ---- poj 1811 : Prime Test

Prime Test Time Limit: 6000MS   Memory Limit: 65536K Total Submissions: 29046   Accepted: 7342 Case Time Limit: 4000MS Description Given a big integer number, you are required to find out whether it's a prime number. Input The first line contains the

Miller_Rabin素数测试

1 #include<iostream> 2 #include<cmath> 3 #include<cstdio> 4 #include<cstring> 5 #include<algorithm> 6 using namespace std; 7 8 long long mul(long long a,long long n,long long mo){ 9 long long ans=0; 10 while (n){ 11 if (n&

Miller_Rabin(米勒拉宾)素数测试算法

首先需要知道两个定理: 1: 费马小定理: 假如p是素数,且gcd(a,p)=1,那么 a(p-1)≡1(mod p). 2:二次探测定理:如果p是素数,x是小于p的正整数,且,那么要么x=1,要么x=p-1. 证明:这是显然的,因为相当于p能整除,也即p能整除(x+1)(x-1). 由于p是素数,那么只可能是x-1能被p整除(此时x=1) 或 x+1能被p整除(此时x=p-1). 接着 如果a^(n-1) ≡ 1 (mod n)成立,Miller-Rabin算法不是立即找另一个a进行测试,而是

素数测试算法(基于Miller-Rabin的MC算法) // Fermat素数测试法

在以往判断一个数n是不是素数时,我们都是采用i从2到sqrt(n)能否整除n.如果能整除,则n是合数;否则是素数.但是该算法的时间复杂度为O(sqrt(n)),当n较大时,时间性能很差,特别是在网络安全和密码学上一般都是需要很大的素数.而从目前来看,确定性算法判断素数的性能都不好,所以可以用MC概率算法来解决,其中Miller Rabin算法就是其中的很经典的解决方法.下面首先介绍下相关的数学理论. 数学原理 Fermat小定理:若n是素数,则对所有1≤a≤n-1的整数a,有a^(n-1)mod

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

POJ1811 给一个大数,判断是否是素数,如果不是素数,打印出它的最小质因数 随机素数测试(Miller_Rabin算法) 求整数素因子(Pollard_rho算法) 科技题 1 #include<cstdlib> 2 #include<cstdio> 3 const int maxn=10005; 4 const int S=20; 5 int tot; 6 long long n; 7 long long factor[maxn]; 8 long long muti_mod(

HDU2138 随机素数测试 Miller-Rabin算法

题目描述 Give you a lot of positive integers, just to find out how many prime numbers there are.. In each case, there is an integer N representing the number of integers to find. Each integer won’t exceed 32-bit signed integer, and each of them won’t be