022Miller-Rabin 概率素性检验 - 概率与数论的完美联姻

Miller-Rabin 概率素性检验 - 概率与数论的完美联姻

米勒-拉宾：锻造密码学的牢不可破之钥

核心发明者：

加里·米勒（Gary L. Miller）：卡内基梅隆大学计算机科学教授，1976年在博士论文中提出了基于广义黎曼猜想（GRH）的确定性素性检验算法
迈克尔·拉宾（Michael O. Rabin，1931-）：以色列-美国计算机科学家，图灵奖得主（1976年，与Dana Scott共同获奖），希伯来大学和哈佛大学教授，1977年将米勒的确定性算法改造为概率版本，使其无需依赖未经证明的猜想

背景：拉宾是理论计算机科学的奠基人之一，他对不确定性图灵机、概率算法有根本性贡献。他对素性检验问题的兴趣源于密码学的需求：RSA密码体制（1977年）需要生成大素数，而当时没有足够快的确定性算法。

时间线：

时代背景：

地点：

算法：Miller-Rabin 概率素性检验（Probabilistic Primality Test）

数学基础：费马小定理的强化版------强概率素数条件

对奇素数 p，将 p-1 写成 2^s × d（d 为奇数）。对任意 a（1 < a < p-1），必有：

复制代码

a^d ≡ 1 (mod p)，或
a^(2^r × d) ≡ -1 (mod p)，对某个 0 ≤ r < s

若 n 是合数，随机选取的 a 通过上述检验的概率至多为 1/4。

核心突破：

Carmichael数问题 ：

费马素性检验（a^(n-1) ≡ 1 mod n）对Carmichael数（如561=3×11×17）失效，因为所有 a 都能通过。Miller-Rabin通过检验"强"条件解决了这个问题。

密码学驱动：RSA、DSA、椭圆曲线密码等公钥体制都需要生成大素数（通常512位到4096位）。朴素试除法对这种量级完全不可行。

效率需求：生成一个1024位的RSA素数，需要测试约710个候选数（素数定理估计）。每次检验必须快速。

可靠性需求：即使是概率算法，也需要误判率可量化、可控制。

时间复杂度：O(k × log²(n) × M(log n))，其中 M(log n) 是大数乘法时间

确定性使用：对小范围的 n，使用固定的基集合可保证100%正确：

历史影响：

需求名称：实现Miller-Rabin概率素性检验，支持概率版本与确定性版本

对应文件 : miller_rabin.c

编译运行:

bash 复制代码

gcc -std=c99 -Wall -o /tmp/miller_rabin_test miller_rabin.c
./miller_rabin_test