小象扑满用Python实现Miller-Rabin素性测试：竞赛选手的高效素数判定指南
在算法竞赛中，素数判定是一个经典问题。当面对需要快速判断大数是否为素数的题目时，传统的试除法往往力不从心。Miller-Rabin素性测试作为一种概率性算法，以其高效性和实用性成为竞赛选手的利器。本文将深入解析该算法的Python实现，并提供可直接用于竞赛的优化代码和参数表。
1. 为什么需要Miller-Rabin算法
素数判定在密码学、数论和算法竞赛中有着广泛应用。试除法虽然简单直观，但对于大数（如10^18量级）效率极低。Miller-Rabin算法则能在O(k log³n)时间内完成判定，其中k是测试次数。
传统方法的局限性：
- 试除法时间复杂度：O(√n)
- 确定性AKS算法虽准确但实现复杂
- Fermat测试无法识别Carmichael数
提示：Carmichael数是指能通过所有基的Fermat测试的合数，如561、1105等
Miller-Rabin通过引入二次探测，有效解决了这些问题。下面是一个简单对比：
| 方法 | 时间复杂度 | 确定性 | 适用场景 |
|---|---|---|---|
| 试除法 | O(√n) | 是 | 小规模数 |
| Fermat测试 | O(k logn) | 否 | 快速筛选 |
| Miller-Rabin | O(k log²n) | 概率性 | 竞赛常用 |
| AKS | O(log⁶n) | 是 | 理论证明 |
需要模型API调用？ 免费领10W Token，多模型网关一键接入 Claude、DeepSeek 等主流模型。
2. Miller-Rabin算法原理详解
Miller-Rabin算法的核心在于利用模算术和二次探测理论。其数学基础是：
定理：若p是奇素数，将p-1分解为2^s * d，则对于任意a(1 < a < p)，以下至少一个成立：
- a^d ≡ 1 mod p
- 存在某个r(0 ≤ r < s)使得a^(2^r * d) ≡ -1 mod p
算法步骤如下：
-
处理特殊情况：
- n < 2：非素数
- n == 2：素数
- n为偶数：非素数
-
分解n-1为2^s * d的形式
-
对于选定的基a：
- 计算x = a^d mod n
- 若x =