An effective analytic recurrence for prime numbers

基本信息

• 论文ID: 2508.02690
• 标题: An effective analytic recurrence for prime numbers
• 作者: Benoit Cloitre
• 分类: math.NT (Number Theory), math.HO (History and Overview)
• 发表时间: 2025年10月12日 (arXiv v2)
• 论文链接: https://arxiv.org/abs/2508.02690v2

摘要

Golomb-Keller公式通过黎曼ζ函数和欧拉乘积，将下一个素数$p_{n+1}$表示为前$n$个素数$p_1, \ldots, p_n$的递推关系，但需要取极限$s \to \infty$，使其无法构造性地使用。本文通过证明有限参数$s \leq p_n$结合向上取整函数即可充分，将这一渐近公式转化为有效递推，建立了对所有$n \geq 1$有效的构造性方法。

最小整数参数$s_n$(OEIS A389650)揭示了与素数星座的深层联系。作者无条件证明了$\liminf_{n\to\infty} \sigma_n = 0$，其中$\sigma_n = s_n/p_n$。上极限$C = \limsup \sigma_n$满足$\log \psi \lesssim C \leq 0.4332$，其中$\psi \approx 1.46557$是超黄金比例。下界依赖于孪生素数猜想；上界是无条件的。常数$C$与最密集的可容许素数星座相关，连接到Hardy-Littlewood猜想。

该方法扩展到Dirichlet L函数，产生了计算$p_{n+1}$的其他有效公式，并能以降低的精度要求预测$p_{n+1}$模任意整数的余数。

研究背景与动机

问题背景

寻找素数的显式公式是数论中的经典问题。虽然存在直接的非递归公式（如Willans公式、Mills公式），但它们在计算上不可行。本文专注于递推关系，即用$p_1, \ldots, p_n$表示$p_{n+1}$。

历史发展

• Gandhi 4 首先使用原数和Möbius函数提供了此类递推
• Vanden Eynden 19 简化了证明
• Jakimczuk 9 推广了方法
• Golomb 5 使用解析数论独立发现了公式，后被Keller 10 重新发现

现有方法的局限性

经典的Golomb-Keller公式为： $p_{n+1} = \lim_{s\to\infty} \left[\left(\prod_{k=1}^n \left(1-\frac{1}{p_k^s}\right)\right) \zeta(s) - 1\right]^{-1/s}$

该公式的主要问题是需要取极限$s \to \infty$，使其在实际计算中无法使用。

研究动机

本文采用相反的方法：保留完整的ζ函数级数计算到工作精度，但使用有限的$s$。这样截断指数而非级数，使公式变为构造性的。

核心贡献

1. 构造性递推公式：证明了对所有$n \geq 1$，存在最小整数$s_n$使得： $p_{n+1} = \left\lceil \left(-1 + \zeta(s_n) \prod_{j=1}^n \left(1-\frac{1}{p_j^{s_n}}\right)\right)^{-1/s_n} \right\rceil$
2. 有效界限：
   • 使用Bertrand假设证明$s_n \leq 2p_n$（定理10）
   • 使用Nagura定理证明$s_n \leq p_n$（定理12）
3. 渐近行为分析：
   • 无条件证明$\liminf_{n\to\infty} \sigma_n = 0$（命题13）
   • 建立$C := \limsup_{n\to\infty} \sigma_n$的界限：$0.3823 \lesssim C \leq 0.4332$
4. 与素数星座的联系：发现下界$\log \psi \approx 0.3823$（依赖孪生素数猜想），其中$\psi$是超黄金比例
5. 扩展到Dirichlet L函数：允许预测$p_{n+1} \pmod 4$等余数性质
6. 数值数据：提供了$n = 1$到$200$的$s_n$值（OEIS A389650）

方法详解

任务定义

给定前$n$个素数$p_1, p_2, \ldots, p_n$，构造性地计算下一个素数$p_{n+1}$。

核心方法架构

1. Dirichlet级数机制

定义关键函数： $D_n(s) = \sum_{\substack{k \geq 1 \\ \gcd(k,P_n)=1}} k^{-s}$

其中$P_n = \prod_{j=1}^n p_j$是第$n$个原数。

2. 欧拉乘积表示

引理3：对$\Re(s) > 1$， $D_n(s) = \zeta(s) \prod_{j=1}^n (1-p_j^{-s})$

3. 收敛性质分析

定义$h(s) = (D_n(s)-1)^{-1/s}$，证明：

• $h(s) < p_{n+1}$对所有$s > 1$成立
• $\lim_{s\to\infty} h(s) = p_{n+1}$
• $h(s)$在$(1,\infty)$上严格递增

4. 临界参数确定

命题6：对每个$n \geq 1$，存在唯一的$s_n^* > 1$使得$h(s_n^*) = p_{n+1} - 1$。

定义$s_n = \lfloor s_n^* \rfloor + 1$为最小整数参数。

技术创新点

1. 精度与截断的权衡：不同于Keller截断级数求和的方法，本文保留完整的$\zeta(s)$但使用有限的$s$
2. 向上取整技巧：巧妙使用$\lceil h(s) \rceil = p_{n+1}$当且仅当$s > s_n^*$的性质
3. 积分界限技术：使用精细的积分比较来控制尾项误差

实验设置

数值计算工具

• 使用PARI/GP和Python的mpmath库进行高精度计算
• 对$n \leq 200$需要100位精度
• 对$n \approx 500$需要约2500位精度（由于消除效应增强）

验证方法

通过直接计算验证理论界限$s_n \leq p_n$对所有$n = 1, \ldots, 200$成立。

工作示例

计算$p_7 = 17$（从$n = 6$开始）：

• 使用$s = 2p_6 = 26$：$h(26) \approx 16.941817904$，得$p_7 = \lceil h(26) \rceil = 17$
• 实际最小值$s_6 = 8$：$h(8) \approx 16.5189076$

实验结果

主要数值结果

界限验证

• 定理10：$s_n \leq 2p_n$对所有$n \geq 1$成立
• 定理12：$s_n \leq p_n$对所有$n \geq 1$成立（通过Nagura定理）

渐近行为

• 命题13：$\liminf_{n\to\infty} \sigma_n = 0$（无条件）
• 定理14：$C = \limsup_{n\to\infty} \sigma_n \leq 0.4332$（无条件）
• 定理15：在孪生素数猜想下，$C > \log \psi \approx 0.38225$

经验分布分析

对$n = 1$到$200$的数据分析显示：

• 移除5个异常值后，$\sigma_n$的分布类似Beta分布
• 均值≈0.291，中位数≈0.277，标准差≈0.087
• 拟合参数：Beta(α≈7.64, β≈18.62)
• 理论众数≈0.274，与经验中位数一致

固定系数公式

定理18：对任何$c > c_0 \approx 0.5956$，存在$N_0(c)$使得对所有$n \geq N_0(c)$，使用$s = cp_n$的公式正确计算$p_{n+1}$。

相关工作

历史发展脉络

1. Gandhi (1971)：首个使用原数和Möbius函数的递推
2. Golomb (1976)：引入解析数论方法
3. Keller (2007)：独立重新发现并提供不同推导
4. 本文 (2025)：首次使公式变为构造性

与其他方法的比较

• Willans公式：直接但计算不可行
• Mills公式：基于常数但需要预知素数
• 筛法：实用但不提供递推结构
• 本方法：理论上有趣但计算复杂度高

结论与讨论

主要结论

1. 构造性突破：首次将Golomb-Keller公式从渐近转化为有效可计算
2. 深层联系：揭示了素数递推参数与素数星座、间隙分布的内在联系
3. 理论界限：建立了参数$s_n$的精确界限和渐近行为
4. 超黄金比例：发现其在素数理论中的新应用

局限性

1. 计算复杂度：$O(np_n^3 \log p_n)$位运算，虽然多项式但实际不可行
2. 精度要求：随着$n$增大需要指数级增长的工作精度
3. 依赖性：部分结果依赖于未证明的猜想（如孪生素数猜想）

未来方向

1. 计算优化：寻找降低精度要求的方法
2. 理论完善：移除对未证明猜想的依赖
3. 推广应用：扩展到其他数论函数
4. 数值探索：计算更大范围的$s_n$值以验证猜想

深度评价

优点

1. 理论创新：成功解决了一个长期存在的构造性问题
2. 方法严谨：使用了深入的解析数论技巧
3. 联系深刻：建立了递推公式与素数星座理论的桥梁
4. 数据丰富：提供了详实的数值验证和统计分析

不足

1. 实用性有限：虽然理论上构造性，但计算成本过高
2. 精度挑战：高精度要求限制了方法的可扩展性
3. 部分结果条件性：某些重要结果依赖未证明猜想

影响力

1. 理论贡献：为素数递推理论提供新视角
2. 方法论价值：展示了如何将渐近公式转化为有效算法
3. 跨领域联系：连接解析数论、计算数论和素数几何

适用场景

1. 理论研究：素数分布、间隙理论、星座研究
2. 数值实验：小规模验证和模式探索
3. 教学展示：解析数论方法的经典应用

参考文献

关键参考文献包括：

• 5 S. W. Golomb, Formulas for the next prime, Pacific J. Math. 63 (1976), 401–404
• 10 J. B. Keller, A recursion equation for prime numbers, arXiv:0711.3940, 2007
• 4 J. M. Gandhi, Formulae for the nth prime, Proc. Washington State Univ. Conf. Number Theory (1971), 96–107
• 12 J. Nagura, On the interval containing at least one prime number, Proc. Japan Acad. 28 (1952), 177–181

这篇论文在理论数论领域具有重要意义，虽然实际计算应用有限，但其理论洞察和方法创新为素数理论研究开辟了新的方向。特别是发现的与超黄金比例和素数星座的联系，展现了数论中意想不到的深层结构。