2025-11-21T02:43:15.649030

An effective analytic recurrence for prime numbers

Cloitre
The Golomb--Keller formula expresses the next prime $p_{n+1}$ as a recurrence relation in terms of the first $n$ primes $p_1, \ldots, p_n$ using the Riemann zeta function and an Euler product, but requires taking a limit as $s \to \infty$, rendering it non-constructive. We transform this asymptotic formula into an effective recurrence by proving that a finite parameter $s \leq p_n$ suffices when combined with the ceiling function, establishing a constructive method valid for all $n \geq 1$. The minimal integer parameter $s_n$ (OEIS A389650) reveals deep connections to prime constellations. We prove $\liminf_{n\to\infty} σ_n = 0$ unconditionally, where $σ_n = s_n/p_n$. The limit superior $C = \limsup σ_n$ satisfies $\log ψ\lesssim C \leq 0.4332$, where $ψ\approx 1.46557$ is the supergolden ratio. The lower bound is conditional on the twin prime conjecture; the upper bound is unconditional. The constant $C$ relates to the densest admissible prime constellation, connecting to the Hardy--Littlewood conjectures. The method extends to Dirichlet L-functions, yielding other effective formulas for calculating $p_{n+1}$ but also for predicting residues of $p_{n+1}$ modulo any integer with reduced precision requirements.
academic

素数に対する有効な解析的漸化式

基本情報

  • 論文ID: 2508.02690
  • タイトル: An effective analytic recurrence for prime numbers
  • 著者: Benoit Cloitre
  • 分類: math.NT (数論)、math.HO (歴史と概観)
  • 発表日時: 2025年10月12日 (arXiv v2)
  • 論文リンク: https://arxiv.org/abs/2508.02690v2

要約

Golomb-Keller公式はリーマンζ関数とオイラー積を通じて、次の素数pn+1p_{n+1}を最初のnn個の素数p1,…,pnp_1, \ldots, p_nの漸化式として表現しますが、極限s→∞s \to \inftyを取る必要があるため、構成的に使用することができません。本論文は、有限パラメータs≤pns \leq p_nと天井関数の組み合わせで十分であることを証明することにより、この漸近公式を有効な漸化式に変換し、すべてのn≥1n \geq 1に対して有効な構成的方法を確立しました。

最小整数パラメータsns_n (OEIS A389650)は、素数星座との深い関連性を明らかにします。著者は無条件にlim inf⁡n→∞σn=0\liminf_{n\to\infty} \sigma_n = 0を証明しました。ここでσn=sn/pn\sigma_n = s_n/p_nです。上極限C=lim sup⁡σnC = \limsup \sigma_nはlog⁡ψ≲C≤0.4332\log \psi \lesssim C \leq 0.4332を満たします。ここでψ≈1.46557\psi \approx 1.46557は超黄金比です。下界はツインプライム予想に依存し、上界は無条件です。定数CCは最も密集した許容可能な素数星座に関連し、Hardy-Littlewood予想と結びついています。

この方法はDirichlet L関数に拡張され、pn+1p_{n+1}を計算するための他の有効な公式が得られ、精度要件を低減してpn+1p_{n+1}を任意の整数で割った余りを予測することができます。

研究背景と動機

問題背景

素数の明示的公式を見つけることは数論における古典的問題です。直接的な非漸化公式(Willans公式、Mills公式など)は存在しますが、計算上実行不可能です。本論文は漸化式、すなわちp1,…,pnp_1, \ldots, p_nからpn+1p_{n+1}を表現することに焦点を当てています。

歴史的発展

  • Gandhi 4は最初に原始数とメビウス関数を使用してこのような漸化式を提供しました
  • Vanden Eynden 19は証明を簡略化しました
  • Jakimczuk 9は方法を一般化しました
  • Golomb 5は解析的数論を使用して独立に公式を発見し、後にKeller 10によって再発見されました

既存方法の限界

古典的なGolomb-Keller公式は以下の通りです: pn+1=lim⁡s→∞[(∏k=1n(1−1pks))ζ(s)−1]−1/sp_{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→∞s \to \inftyを取る必要があり、実際の計算では使用できないことです。

研究動機

本論文は逆のアプローチを採用しています。完全なζ関数級数を作業精度まで計算しますが、有限のssを使用します。このように指数を切り詰め、級数ではなく、公式を構成的にします。

核心的貢献

  1. 構成的漸化公式: すべてのn≥1n \geq 1に対して、最小整数sns_nが存在して以下が成り立つことを証明しました: pn+1=⌈(−1+ζ(sn)∏j=1n(1−1pjsn))−1/sn⌉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仮説を使用してsn≤2pns_n \leq 2p_nを証明(定理10)
    • Nagura定理を使用してsn≤pns_n \leq p_nを証明(定理12)
  3. 漸近挙動の分析:
    • 無条件にlim inf⁡n→∞σn=0\liminf_{n\to\infty} \sigma_n = 0を証明(命題13)
    • C:=lim sup⁡n→∞σnC := \limsup_{n\to\infty} \sigma_nの界を確立: 0.3823≲C≤0.43320.3823 \lesssim C \leq 0.4332
  4. 素数星座との関連性: 下界log⁡ψ≈0.3823\log \psi \approx 0.3823(ツインプライム予想に依存)を発見。ここでψ\psiは超黄金比です
  5. Dirichlet L関数への拡張: pn+1(mod4)p_{n+1} \pmod 4などの余数特性を予測することが可能になります
  6. 数値データ: n=1n = 1から200200までのsns_n値を提供(OEIS A389650)

方法の詳細

タスク定義

最初のnn個の素数p1,p2,…,pnp_1, p_2, \ldots, p_nが与えられたとき、構成的に次の素数pn+1p_{n+1}を計算します。

コア方法アーキテクチャ

1. Dirichlet級数メカニズム

主要関数を定義します: Dn(s)=∑k≥1gcd⁡(k,Pn)=1k−sD_n(s) = \sum_{\substack{k \geq 1 \\ \gcd(k,P_n)=1}} k^{-s}

ここでPn=∏j=1npjP_n = \prod_{j=1}^n p_jは第nnの原始数です。

2. オイラー積表現

補題3: ℜ(s)>1\Re(s) > 1に対して、 Dn(s)=ζ(s)∏j=1n(1−pj−s)D_n(s) = \zeta(s) \prod_{j=1}^n (1-p_j^{-s})

3. 収束特性の分析

h(s)=(Dn(s)−1)−1/sh(s) = (D_n(s)-1)^{-1/s}を定義し、以下を証明します:

  • すべてのs>1s > 1に対してh(s)<pn+1h(s) < p_{n+1}
  • lim⁡s→∞h(s)=pn+1\lim_{s\to\infty} h(s) = p_{n+1}
  • h(s)h(s)は(1,∞)(1,\infty)上で厳密に増加

4. 臨界パラメータの決定

命題6: 各n≥1n \geq 1に対して、h(sn∗)=pn+1−1h(s_n^*) = p_{n+1} - 1を満たす唯一のsn∗>1s_n^* > 1が存在します。

sn=⌊sn∗⌋+1s_n = \lfloor s_n^* \rfloor + 1を最小整数パラメータとして定義します。

技術的革新点

  1. 精度と切り詰めのバランス: Keller級数和を切り詰める方法と異なり、本論文は完全なζ(s)を保持しながら有限のssを使用します
  2. 天井関数のトリック: ⌈h(s)⌉=pn+1\lceil h(s) \rceil = p_{n+1}当且つつs>sn∗s > s_n^*の場合のみという性質を巧妙に使用
  3. 積分界限技術: 精密な積分比較を使用して尾項誤差を制御

実験設定

数値計算ツール

  • PARI/GPおよびPythonのmpmath ライブラリを使用した高精度計算
  • n≤200n \leq 200に対して100桁の精度が必要
  • n≈500n \approx 500に対して約2500桁の精度が必要(消去効果の増強による)

検証方法

理論的界限sn≤pns_n \leq p_nがすべてのn=1,…,200n = 1, \ldots, 200に対して成立することを直接計算により検証します。

作業例

p7=17p_7 = 17を計算(n=6n = 6から開始):

  • s=2p6=26s = 2p_6 = 26を使用: h(26)≈16.941817904h(26) \approx 16.941817904、p7=⌈h(26)⌉=17p_7 = \lceil h(26) \rceil = 17を得る
  • 実際の最小値s6=8s_6 = 8: h(8)≈16.5189076h(8) \approx 16.5189076

実験結果

主要な数値結果

界限の検証

  • 定理10: すべてのn≥1n \geq 1に対してsn≤2pns_n \leq 2p_n
  • 定理12: すべてのn≥1n \geq 1に対してsn≤pns_n \leq p_n(Nagura定理による)

漸近挙動

  • 命題13: lim inf⁡n→∞σn=0\liminf_{n\to\infty} \sigma_n = 0(無条件)
  • 定理14: C=lim sup⁡n→∞σn≤0.4332C = \limsup_{n\to\infty} \sigma_n \leq 0.4332(無条件)
  • 定理15: ツインプライム予想の下で、C>log⁡ψ≈0.38225C > \log \psi \approx 0.38225

経験的分布分析

n=1n = 1から200200までのデータ分析は以下を示します:

  • 5つの外れ値を除去した後、σn\sigma_nの分布はベータ分布に類似
  • 平均≈0.291、中央値≈0.277、標準偏差≈0.087
  • 適合パラメータ: Beta(α≈7.64, β≈18.62)
  • 理論的最頻値≈0.274、経験的中央値と一致

固定係数公式

定理18: c>c0≈0.5956c > c_0 \approx 0.5956に対して、N0(c)N_0(c)が存在して、すべてのn≥N0(c)n \geq N_0(c)に対してs=cpns = cp_nを使用する公式がpn+1p_{n+1}を正しく計算します。

関連研究

歴史的発展の流れ

  1. Gandhi (1971): 原始数とメビウス関数を使用した最初の漸化式
  2. Golomb (1976): 解析的数論方法の導入
  3. Keller (2007): 独立した再発見と異なる導出の提供
  4. 本論文 (2025): 公式を初めて構成的にする

他の方法との比較

  • Willans公式: 直接的だが計算不可能
  • Mills公式: 定数に基づくが素数の事前知識が必要
  • ふるい法: 実用的だが漸化構造を提供しない
  • 本方法: 理論的に興味深いが計算複雑度が高い

結論と考察

主要な結論

  1. 構成的突破: Golomb-Keller公式を初めて漸近から有効計算可能に変換
  2. 深層的関連性: 素数漸化パラメータと素数星座、間隙分布の内在的関連性を明らかに
  3. 理論的界限: パラメータsns_nの精密な界限と漸近挙動を確立
  4. 超黄金比: 素数理論における新しい応用を発見

限界

  1. 計算複雑度: O(npn3log⁡pn)O(np_n^3 \log p_n)ビット演算。多項式的ですが実際には実行不可能
  2. 精度要件: nnの増大に伴い指数級の作業精度増加が必要
  3. 依存性: 部分的な結果は未証明予想(ツインプライム予想など)に依存

今後の方向

  1. 計算最適化: 精度要件を低減する方法の探索
  2. 理論的完成: 未証明予想への依存性を排除
  3. 推広応用: 他の数論関数への拡張
  4. 数値探索: より大きな範囲のsns_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

本論文は理論数論の分野において重要な意義を持ちます。実際の計算応用は限定的ですが、その理論的洞察と方法的革新は素数理論研究に新しい方向を切り開いています。特に発見された超黄金比と素数星座との関連性は、数論における予期しない深層構造を示しています。