[论文解读] Computing $\pi(x)$ Analytically

David J. Platt|arXiv (Cornell University)|Mar 26, 2012

Analytic Number Theory Research参考文献 10被引用 9

一句话总结

本文提出了一种严格的解析算法，使用黎曼的显式公式和区间算法，无条件地计算至 x = 10^24 的 π(x)，即素数计数函数。该方法通过利用误差函数平滑的梅林变换对，改进了先前的方法，确保通过验证的误差界和区间算法实现数值稳定性和正确性，首次实现了对 π(10^24) 的无条件计算。

ABSTRACT

We describe a rigorous implementation of the Lagarias and Odlyzko Analytic Method to evaluate the prime counting function and its use to compute unconditionally the number of primes less than $10^{24}$.

研究动机与目标

• 开发一种严格且无条件的方法，用于计算传统组合筛法无法企及的 π(x)。
• 克服高精度解析数论计算中数值不稳定性和浮点误差累积的限制。
• 使用区间算法和显式误差界，实现并验证拉加里亚斯-奥德利兹科解析方法的变体。
• 无条件地计算 π(10^24)，解决计算素数定理领域长期存在的开放问题。
• 为未来高精度 π(x) 评估提供可重现、可验证且数值鲁棒的框架。

提出的方法

• 采用佩罗公式，将 π*(x) 表示为 log ζ(s) 的复积分，乘以由互补误差函数导出的平滑核函数 φ(t)。
• 使用加瓦伊定义的梅林变换对：φ(t) = ½ erfc( log(t/x)/(√2 λ) ) 和 φ̂(s) = x^s / s ⋅ exp(λ²s²/2)，在收敛性和局部化之间取得平衡。
• 使用区间算法在整个计算过程中严格界定浮点误差，取代启发式误差估计。
• 应用对非平凡黎曼 zeta 函数零点求和截断误差的显式界，结合黎曼-冯·曼戈尔特函数 N(t) 和尾部估计。
• 使用引理 A.1 通过已知的渐近展开式和误差项界 N(t)，实现对零点求和截断的精确控制。
• 针对 [T₁, T₂] 内的零点和 [−T₂, T₂] 外的零点，分别引入误差界 E₁ 和 E₂，利用斯蒂尔杰斯积分和指数衰减估计。

实验结果

研究问题

• RQ1 拉加里亚斯-奥德利兹科解析方法能否通过使用区间算法和严格误差界，实现无条件计算大规模 π(x)？
• RQ2 在高精度解析数论计算中，如何控制数值不稳定性与浮点误差累积？
• RQ3 平滑参数 λ 与截断阈值 T₁ 和 T₂ 之间应如何权衡，以在保持计算可行性的同时最小化误差？
• RQ4 能否对非平凡 zeta 函数零点求和施加数学上严格的误差界，同时考虑已知和未验证的零点位置？
• RQ5 是否可能在不假设黎曼猜想的前提下，使用解析方法和经验证的数值技术，无条件计算 π(10^24)？

主要发现

• 本文成功无条件地计算出 π(10^24)，标志着在该规模上首次实现此类计算，且不依赖黎曼猜想。
• 实现过程中使用区间算法严格控制浮点误差，确保即使在硬件或编译器异常情况下仍保持正确性。
• 该方法的时间复杂度为 O(x^{1/2+ε})，显著优于组合筛法在大 x 下的 O(x^{2/3} log^{-2}x) 复杂度。
• 通过黎曼-冯·曼戈尔特函数 N(t) 导出零点求和截断的显式误差界，引理 A.1 提供了对零点分布的紧密控制。
• 被忽略零点的误差贡献由 E₁ 和 E₂ 界定，其中 E₂ 因核函数的指数衰减而可忽略不计，即使实部不为 ½ 亦然。
• 最终结果通过多种独立验证手段确认，包括与较小尺度下已知的 π(x) 值交叉验证，以及与素数定理的一致性。