梅森素数是形如2^p−1的素数且p本身必须为素数;判定需先验证n是否等于2^p−1(即n+1为2的幂),再检查p是否为素数,最后验证n是否为素数。

美利坚强素数(Mersenne prime)不是标准数学术语,实际应为“梅森素数”(Mersenne prime),指形如 2^p - 1 的素数,其中 p 本身也必须是素数。C++ 中无法直接“判断一个任意整数是否为梅森素数”,因为梅森素数有严格形式约束——它必须是 2^p - 1 形式且该值为素数。所以第一步永远是:先验证输入是否符合梅森数结构。
如何快速判断 n 是否为梅森数(即是否等于 2^p - 1)
对给定整数 n,需检查是否存在整数 p > 1,使得 n == (1 。注意:<code>n 必须为正整数,且 n + 1 必须是 2 的幂。
-
n + 1必须大于 1 且为 2 的整数次幂 → 可用(n + 1) & n == 0判断(前提是n > 0) - 若成立,再用
std::bit_width((unsigned long long)(n + 1)) - 1得到p(C++20),或手动右移计数 - 若
n == 0或n == 1,直接排除(2^1 - 1 = 1不是素数;2^2 - 1 = 3才是第一个梅森素数) - 常见误判:把
31(=2^5 - 1)当成梅森素数是对的,但把15(=2^4 - 1)当成梅森素数就错了——因为4不是素数,15也不是素数
为什么 p 必须是素数?
这是梅森素数定义的核心条件。若 p 是合数(如 p = ab,a,b > 1),则 2^p - 1 = (2^a)^b - 1 可被 2^a - 1 整除,必然为合数。因此即使 n = 2^p - 1 是素数,若 p 不是素数,n 就不被视为梅森素数。
- 例如:
p = 4→n = 15,虽满足形式但p非素,且15非素 → 排除 -
p = 11→n = 2047 = 23 × 89,虽p是素数,但n不是素数 → 不是梅森素数 - 所以完整判定链是:
n→ 检查是否为2^p - 1→ 提取p→ 检查p是否为素数 → 检查n是否为素数
小整数范围可用试除法,大数需注意溢出和性能
对 n 做素性判定时,n 可能很大(如 2^31 - 1 = 2147483647)。C++ 内置类型容易溢出,int 最大仅约 2e9,long long 最大约 9e18 —— 而第 10 个梅森素数已超 30 位十进制数。
立即学习“C++免费学习笔记(深入)”;
- 对
n ,可用优化试除:只试到 <code>sqrt(n),且跳过偶数、模 6±1 数 - 对更大的
n,必须用 Miller-Rabin 概率素性测试(确定性版本对uint64_t范围可行,需预设底数如{2, 325, 9375, 28178, 450775, 9780504, 1795265022}) - 不要用
std::sqrt(n)直接比较循环上限——浮点误差可能导致越界;改用i 或整数平方根二分 - 若输入是
int类型,2^31 - 1已超出int表示范围(有符号 int 通常上限为2^31 - 1,但左移1 是未定义行为)→ 务必用 <code>unsigned或uint64_t
真正麻烦的不是写判定逻辑,而是:你得先接受——绝大多数输入根本连梅森数的边都沾不上。先做 n + 1 是否为 2 的幂,再验 p 是否为素,最后验 n 是否为素。三步缺一不可,且每一步都有隐式类型陷阱和边界 case。


















