素数判定算法 MILLER RABIN

入门级筛素数--试除法，复杂度O（n^2）

bool rmprime( long long n )
{
    for(long long i = 2; i <= sqrt(n) ; i++)    if(n%i == 0) return false;
    return true;
}

 

学了一段时间算法以后，应该会了解到筛法求素数，复杂度略高于O（n）

#include
#include
using namespace std;
const int MAXN = 500000;
bool isp[MAXN];
int p[MAXN];
int cnt = 0;//保存素数个数
void getp()
{
    for(int i = 1; i < MAXN; i++)    isp[i] = 1 ; //默认全部素数
    for(int i = 2; i < MAXN; i++){
        if(!isp[i]) continue; //如果不是素数就continue掉
        p[cnt++] = i;
        for(int j = i * 2 ; j < MX; j += i)        isp[j] = 0;//记录，往后滚处理
    }
}

int main(){
    getp();
    for(int i = 0; i < cnt; i++)    printf("%d ",p[i]);
        return 0;      
}

 当然，不难发现，如果MaX值过大的话，不只空间会炸，而且从头开始扫很玄学，是不怎么可取的。

于是引入MILLER RABIN算法。可以单独判断一个大数是否素数，但是不保证正确。我们只能通过多次执行算法让这个错误的概率很小。//不过幸运的是，除非你是非酋王，这个算法一般都是不会出问题的。

算法的理论基础：

1. Fermat定理：若n是奇素数，a是任意正整数(1≤ a≤ n−1)，则 a^(n-1) ≡ 1 mod n。

2. 推演自Fermat定理（具体过程我没看懂，Orz）, 如果n是一个奇素数，将n−1表示成2^s*r的形式，r是奇数，a与n是互素的任何随机整数，那么a^r ≡ 1 mod n或者对某个j (0 ≤ j≤ s−1, j∈Z) 等式a^(2jr) ≡ −1 mod n 成立。

 

聪明的读者已经发现上述的定理是一个数是素数的必要条件而非充分条件，所以这就是这个算法的不精确的原有，对于一些特定的检验算子a，存在一些合数也满足上述定理。所以我们要多找几个a反复检验，这样能让错误率大大降低。

 

那么我们按照上述的定理2，首先重复n次实验。对于每一次实验，随机取检验算子a，带入定理2进行检验，看看在算子a下，n能否满足

a^r ≡ 1 mod n或者对某个j (0 ≤ j≤ s−1, j∈Z) 等式a^(2jr) ≡ −1 mod n       **

如果任意一次实验不满足，则判定不是素数，如果都满足，可近似可以认为是素数（错误率极小）。

 

#include  
#include  
#include   
#include   
#include   
#include 

 

 

转载于:https://www.cnblogs.com/shuri/p/7638588.html