大素数检测: Miller_Rabin算法

最新推荐文章于 2021-11-11 20:53:37 发布

原创最新推荐文章于 2021-11-11 20:53:37 发布 · 898 阅读

3 ·

CC 4.0 BY-SA版权

文章标签：

#算法

本文介绍了一种解决哥德巴赫猜想的方法，即任一大于5的整数可以写成三个质数之和。文章通过埃拉托斯特尼筛法找出小范围内的质数，并使用Miller-Rabin算法验证较大数的素性。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

2021年ACM竞赛班训练(三)

A题: 验证哥德巴赫猜想

原题链接

思路：

1、原初哥德巴赫猜想的现代陈述为：任一大于5的整数都可以拆成三个质数之和。通过观察样例输入，不妨假设输入 $n$ 足够大时： $n = a + b + c$ 其中 $a, b$ 在小质数内可实现（ $\in (2, 600)$ ，可利用（埃拉托斯特尼）筛法得到该范围内素数表).

2、对于较大数 $c = n - a - b$ 的素性判断考虑Miller-Rabbin算法.

代码：

#include <bits/stdc++.h>
#define LL long long
using namespace std;

LL pl[601]; //用来装筛法得到的素数pl[i] = 0表示素数，pl[i] = 1表示被筛掉的合数。

//qmul用于计算两数的积，定义该函数目的是防溢出
LL qmul(LL x, LL y, LL mod) { 
    return (x * y - (long long)(x / (long double)mod * y + 1e-3) *mod + mod) % mod;
}
//快速幂算法
LL qpow(LL a, LL n, LL mod) {
    LL res = 1, tmp = a;
    for(LL i = n; i; i >>= 1){
        if(i&1) res = (res*tmp)%mod;
        tmp = (tmp*tmp) %mod;
    }
    return res;
}
//Miller_Rabin算法
LL prime[6] = {2, 3, 5, 233, 331}; 
bool Miller_Rabin(LL p) {
    if(p < 2) return 0;
    if(p != 2 && p % 2 == 0) return 0;
    LL s = p - 1;
    while(! (s & 1)) s >>= 1;
    for(int i = 0; i < 5; ++i) {
        if(p == prime[i]) return 1;
        LL t = s, m = qpow(prime[i], s, p);
        while(t != p - 1 && m != 1 && m != p - 1) {
            m = qmul(m, m, p);
            t <<= 1;
        }
        if(m != p - 1 && !(t & 1)) return 0;
    }
    return 1;
}

void solve(LL num) {
	//当a, b, c中a, b为小范围内的素数且c(=n-a-b)通过Miller_Rabin检测时即为答案输出
    for(int i = 2; i <= 600; i++) {
        for(int j = i; j <= 600; j++) {
            if(pl[i] == 0 && pl[j] == 0 && Miller_Rabin(num - i - j)) {
                printf("%d %d %lld\n", i, j, num - i -j);
                return;
            }
        }
    }
}

int main() {
	//埃拉托斯特尼筛法实现
    for(int i = 2; i <= sqrt(600) + 1; i++) {
        if(!pl[i]) for(int j = i*i; j <= 600; j+=i) {
            pl[j] = 1;
        }
    }
	
	//读取数据
    LL num;
    while(~scanf("%lld", &num)) {
        solve(num);
    }
    return 0;
}

Miller_Rabin算法：

已知:

（1）费马定理：
$a^{p-1} = 1 (mod\,p)$
其中 $p$ 为质数， $a$ 不是 $p$ 的倍数
（2）二次探测定理：
$x^2 \equiv 1 (mod\,p) \Rightarrow x = 1 \,or\, x = p-1$

简要证明：

设 $n$ 是奇素数，则 $n - 1$ 必为偶数。令 $n-1 = 2^qm$ .
对 $\forall a \in (1, n)$ , 利用费马定理知, $(a^{2^qm} = a ^{n-1}) \equiv 1 (mod\,n)$ 。
由二次探测定理可知： $a^{2^{q-1}m} \equiv 1 (mod\,n)$ 或 $a^{2^{q-1}m} \equiv n - 1(mod\,n).$
若 $a^{2^{q-1}m} \equiv 1 (mod\,n)$ 成立，则再次由二次探测定理可知:
$a^{2^{q-2}m} \equiv 1 (mod\,n)\,$ 或 $\,a^{2^{q-2}m} \equiv n-1(mod\,n)$
$\cdots$ 。反复利用刺耳探测定理直到：
$a^{m} \equiv 1 (mod\,n)\,$ 或 $a^{m} \equiv n-1 (mod\,n)$ 。
综上：如果 $n$ 是素数，则必有 $a^m = 1\,(mod\,n)$ ，或存在 $a^{2^rm} \equiv n-1(mod\,n)\quad r \in [0, q-1]$ 而使二次检测不能传递到 $m$ 。

算法应用

1、给定（input）奇数 $n$ , 为判断其素性，先测试 $a^{2^pm} \equiv 1(mod\,n)$ 是否成立。若不成立则 $n$ 必为合数，若成立需要进一步测试。
2、考虑对于Miller序列：
$a^m\,(mod\,n), a^{2m}\,(mod\,n),a^{4m}\,(mod\,n), \cdots, a^{2^qm}\,(mod\,n)$ ，若 $a^m\,(mod\,n)$ 成立，或存在 $a^{2^rm} \equiv n-1(mod\,n)\quad r \in [0, q-1]$ 则认为n通过Miller_Rabin测试。

注意：1、Miller_Rabin是素数的必要条件，即素数一定可以通过Miller_Rabin测试，不能通过Miller_Rabin的数必为合数。
2、Miller_Rabin不是素数的充分条件，即通过该测试的数可能不是素数。但该算法给出错误结果的概率 $\leq (\frac{1}{4})^k$ 其中 $k$ 为测试次数。在实际应用中，效果往往好得多！