Miller_Rabin算法详解

目录

基本引理：

1，费马定理：

2，二次探测定理：

作用：

证明：

代码实现：

---

目录

基本引理：

1，费马定理：

2，二次探测定理：

---

基本引理：

1，费马定理：

费马定理的证明链接

2，二次探测定理：

二次探测定理的证明链接

作用：

有效的检测大整数是否为素数。

证明：

由费马定理，可以排除大部分非素数的情况（满足费马定理是素数的必要条件），给出一个奇素数n，显然n-1为一个偶数，存在,显然（q，m为任意整数）是成立的，所以,显然.

所以根据二次探测定理，可以推断,……,.

所以我们只需要证明大整数在每次分解都满足上式就可以大概率的确定n是素数。每次判断的错误率给、大概为0.25，k次运算后的错误率大概是0.25^k,大概8次以上这个错误概率就低到忽略不计了

代码实现：

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<ctime>
using namespace std;
typedef long long ll;
const int N = 1e5 + 7;
const int times = 10;
ll fast_mod(ll a,ll b,ll mod)//计算2^q的过程
{
    ll res = 0;
    while(b){
        if(b & 1) res = res + a;
        a <<= 1;
        if(a >= mod) a -= mod;
        if(res >= mod) res -= mod;
        b >>= 1;
    }
    return res;
}
ll fast_pow_mod(ll a,ll b,ll mod)//快速幂算出a^m
{
    ll res = 1;
    while(b){
        if(b & 1) res = (res * a) % mod;
        a = (a * a) % mod;
        b >>= 1;
    }
    return res;
}
bool check(ll a,ll m,ll p,ll n)//对于每次随机的a进行测试
{
    ll temp = fast_pow_mod(a,m,n),ret = temp;
    for(int i = 0;i < p;++i){
        ret = fast_mod(temp,temp,n);
        if(ret == 1 && temp != n - 1 && temp != 1) return true;
        temp = ret;
    }
    return ret != 1;
}
bool Miller_Pabin(ll n)//Miller测试的主体结构
{
    if(n < 2) return false;
    if(n == 2) return true;
    if(n & 1 == 0) return false;//对于偶数的优化
    ll p = 0,x = n - 1;//p为Miller测试的q，x为Miller测试的m
    while(x & 1 == 0){
        x >>= 1;
        p++;
    }
    srand(time(NULL));
    for(int i = 0;i < times;++i){
        ll o = rand() % (n - 1) + 1;//o就是Miller测试的底数a
        if(check(o,x,p,n)) return false;
    }
    return true;
}

int main()
{
    ios::sync_with_stdio(false);
    int t;
    cin >> t;
    while(t--){
        long long n;
        cin >> n;
        cout << (Miller_Pabin(n) ? "Prime" : "Not a Prime") << endl;
    }
    return 0;
}