码迷,mamicode.com
首页 > 其他好文 > 详细

【数论】Miller_Rabin

时间:2016-11-27 07:08:41      阅读:147      评论:0      收藏:0      [点我收藏+]

标签:判断   步骤   pre   class   long   for   turn   提高   iostream   

Miller_Rabin素数测试

    Miller_Rabin判断单个素数的方法运用了费马小定理,可以说非常之快了。

    Miller_Rabin曾经被称作“黑科技”,但是根据费马小定理其实完全可以自己写出来大半。

 

    其算法的运行过程如下:

    (1)对于奇数M,使得N=(2^r)*M+1

    (2)选取随机数使得A<N

    (3)对于任意i(i<r),若(A^(2^i)) Mod N=N - 1,则N为素数

    (4)或者,若(A^M) Mod N=1,则N通过随机数A的测试

    若对素数N进行T次测试,那么失误率为1/4^T,我们可以进一步提高其效率,如省去步骤3

   

代码如下:

#include<iostream>
#include<cstring>
#include<cstdio>
#include<time.h>
#include<cstdlib>
#include<cmath>
using namespace std;
long long exp(long long a,long long m,long long n){//快速幂
	if(m==0) return 1;
	if(m==1) return (a%n);
	long long w=exp(a,m/2,n);
	w=w*w%n;
	if(m&1) w=w*a%n;
	return w%n;
}
bool Witness(long long a,long long n)
{
    long long m=n-1;//满足原先条件
    int j=0;
    while(!(m&1)){
        j++;
        m>>=1;
    }
    long long x=exp(a,m,n);
    if(x==1||x==n-1) return false;
    while(j--){
        x=x*x%n;
        if(x==n-1) return false;
    }
    return true;
}
bool Miller_Rabin(long long n){
	if(n==2) return true;
	if(n&1==0) return false;
	for(int i=1;i<=10;i++){
		long long a=rand()%(n-2)+2;//一定为a<N
		if(Witness(a,n)) return false;
	}
	return true;
}
bool prime(long long N){
	long long k=sqrt(N);
	for(int i=2;i<=k;i++) if(N%i==0) return false;
	return true;
}
int main(){
    srand(time(NULL));
    for(long long i=3;i<=10000000;i++)
    if(Miller_Rabin(i)!=prime(i)) cout<<i<<endl;
}

  

【数论】Miller_Rabin

标签:判断   步骤   pre   class   long   for   turn   提高   iostream   

原文地址:http://www.cnblogs.com/wxjor/p/6105554.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!