Research & findings

Investigations, surveys, third-party tool comparisons, and concrete decisions distilled from the conversation.

多 agent 4 阶段头脑风暴：研究路线图

2026-05-23 · session 69442929

10 题 × 4 阶段（10+5+2+1）的发散→综合架构，最终产出 5 个跨题大主题、3 层时间目标、5 个高优先级可攻命题。

  • 5 个跨题大主题：Frobenius 轨道结构化 / 指数-素性 Dirichlet 级数 D(s) / 条件反证硬度分级 / 大素因子枢纽目标 / 算术随机性语言建设
  • 短期（5-10y）目标：大素因子下界提升 / Bateman-Horn 奇异级数正性 / GRH 下本原素因子矛盾
  • 中期（10-20y）：D(s) 解析延拓 / 指数序列有效筛法 / Lucas-Lehmer Galois 字典
  • 长期（20+y）：算术随机性公理化 / 指数型 Schinzel H 弱化版本
  • 5 个高优先级命题：本原素因子密度下界（无条件）/ M_p 大素因子比例下界（GRH）/ Bateman-Horn 严格正性（无条件，可立即动笔）/ 有限假设导致 Chebotarev 矛盾（GRH）/ Mersenne 素性序列算术不可压缩性（position paper）
  • 局限性诚实告示：文献遗漏 / 伪深刻 / 技术跳跃 / 独立性假设 / 策略覆盖不全
Sources: /tmp/mersenne_brainstorm/phase1_q01..q10.md, /tmp/mersenne_brainstorm/phase2_a1..a5.md, /tmp/mersenne_brainstorm/phase3_a1.md, phase3_a2.md, /tmp/mersenne_brainstorm/phase4_section13.html

Mersenne 素数无穷性问题：基础与 Wagstaff 启发式

2026-05-23 · session 69442929

整理 Mersenne 素数无穷性猜想的命题、必要条件、Wagstaff 启发式的推导链以及 EFF 奖金现状。

  • M_p = 2^p - 1 形式；p 必须为素数（n 合数 ⟹ M_n 必合，反之不真：M_11 = 23×89）
  • Lenstra-Pomerance-Wagstaff: #{Mersenne primes with exponent ≤ N} ~ (e^γ/ln 2)·ln N，常数 ≈ 2.5695
  • ln ln N 类极慢函数：Mersenne 永远稀疏但无界
  • EFF 奖金：10^7 位（已领 2009）；10^8 位 \$150k（未领，预计 2030 年代中期）；10^9 位 \$250k
Sources: Wagstaff 1983, GIMPS, EFF Cooperative Computing Awards

如何验证 M_136279841（4100 万位）的素性

2026-05-23 · session 69442929

Lucas-Lehmer 检验是确定性多项式时间算法，依靠 FFT 把朴素的 2^{p/2} 试除降到 p²log p。代数数论（ω=2+√3 在 Q(√3) 中的阶）解释为何这种特殊形式数反而易验。

  • Lucas-Lehmer: s_0=4, s_{k+1}=s_k²-2 mod M_p；M_p 素 ⟺ s_{p-2} ≡ 0
  • 工程：FFT 乘法 O(p log p)；模 M_p 等价位旋转（无除法）；总成本约 4×10^17 位运算
  • 正确性保障：Gerbicz 错误检测 + 不同硬件/软件双独立复跑 + 公布末位 64 比特
  • M_136279841 由 Luke Durant 用 Nvidia A100 云集群约 30 GPU·天验证
Sources: Lucas 1878, Lehmer 1930, Schönhage-Strassen, Gerbicz 2017, GIMPS / gpuOwL

现有方法是暴力还是解析？

2026-05-23 · session 69442929

把'验证单个 M_p'、'预测密度'、'证明无穷性'三件事分清。前者是代数算法（非暴力），中者是启发式（非证明），后者完全无解析入口。

  • Lucas-Lehmer 不是暴力，加速来自代数数论
  • 解析方法全部失效：Euclid 不能递推；Dirichlet/L 函数不适用指数稀疏序列；筛法要求密度 ≥ 1/log^C x 而 Mersenne 密度仅 ln ln x / x
  • 重述存在（Φ_p(2) ⟺ Schinzel H、动力系统 T(x)=x²-2）但都同等困难
  • AKS / ECPP 通用素性测试常数太大，对 4×10^7 位数完全够不到
Sources: Friedlander-Iwaniec 1997, Maynard-Tao 2013, Schinzel H hypothesis

GIMPS 怎么挑选要测试的 p ——筛选漏斗

2026-05-23 · session 69442929

实际工作流是 6 阶段筛选：素数枚举 → 试除 → P-1 → ECM → PRP → LL → 双独立复跑。每阶段越来越贵但也越来越准；进入 LL 的只剩约 10%。

  • 试除（TF）利用因子结构 q = 2kp+1 且 q ≡ ±1 (mod 8)，秒级单价砍掉 ~65%
  • Pollard P-1 撬出 q-1 光滑因子（典型 B_1 ~ 10^6、B_2 ~ 3×10^7）
  • PRP（2018 起取代 LL 为默认主测试）：3^{M_p-1} ≡ 1 mod M_p，支持 Gerbicz 中间校验
  • PrimeNet 服务器分发任务，按硬件能力分发不同阶段任务
  • 当前主战线 p ∈ [10^8, 1.2×10^8]，p ≤ 10^8 已全测两遍
Sources: GIMPS / Prime95, gpuOwL, mlucas, PrimeNet

52 个已知 Mersenne 素数：完整列表 + 增长分布图

2026-05-23 · session 69442929

整理了从古代发现的 M_2=3 到 2024 年发现的 M_136279841 的全部 52 个，配 Chart.js 渲染的 log10(p) vs n 散点 + Wagstaff 直线拟合 + 累计发现年表。

  • 前 48 个序号已确认；#49–#52 为暂定序号（中间可能尚有未发现）
  • 1952 年 Robinson 用 SWAC 一次发现 5 个（#13–17）；1996 年 GIMPS 启动后曲线匀速上升
  • log10(p_n) vs n 实际数据几乎完美贴合 Wagstaff 直线斜率 0.169
  • #29 vs #30、#46 vs #47 出现序号-发现年逆置
  • #47（UCLA Smith 2008）赢 EFF \$100k；#52（Durant 2024）首个非个人 PC 发现
Sources: GIMPS history, Mersenne Wiki

谁在做这件事 · 真实动机

2026-05-23 · session 69442929

学术数学界几乎不参与；GIMPS 是唯一的有组织搜索项目。真实动机不是'证明无穷'，而是硬件压力测试 / 分布式计算示范 / 算法工程练兵场。

  • GIMPS 1996 创立（早 SETI@home 3 年），是史上首个 BOINC 风格志愿项目
  • Prime95 是 Intel/AMD/Nvidia 内部 CPU/GPU 极限稳定性测试标准
  • EFF 设奖目的是推广合作计算理念 ——Mersenne 是载体不是目标
  • 主要发现者：Curtis Cooper、Edson Smith、Luke Durant、Slowinski/Gage、Robinson
  • Maynard、Tao、Granville、Soundararajan 等顶级解析数论人都不在做这个问题
  • 中国算力贡献 8-12%（次于美/德），但还无中国发现者
Sources: mersenne.org, EFF Cooperative Computing Awards, PrimeNet stats

近 12 年学界文献全景（arXiv 2014-2026）

2026-05-23 · session 69442929

通过 arXiv API 拉取 47 篇直接涉及 Mersenne 的论文，分 6 类整理。0 篇声称证明无穷性。函数域类比和算法改进是主要正经方向。

  • 0 篇声称证明无穷性
  • 2 篇函数域类比有实质内容（Lucas 2025、Nguyen 2014）
  • 2 篇等价重述：Shlossberg 2022（拓扑群极小性）、Prunescu 2025（统一计数封闭式）
  • 约 8 篇改进 Lucas-Lehmer 或周边算法
  • 约 10 篇研究 M_p 在其他对象（完美多项式、关联方案）中的角色
  • 5-10 篇 math.GM 类未审稿稿件，多为民科作品
  • 已下载 17 篇 PDF 到 ./reference_mersenne/
Sources: arXiv 2512.08060 (Lucas 2025), arXiv 2204.08302 (Shlossberg 2022), arXiv 2512.01680 (Prunescu 2025), arXiv 2603.08994 (Dominguez 2026), arXiv 1407.7206 (Nguyen 2014)