[論文レビュー] 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)。収束性と局在化のバランスを取る。
• 浮動小数点誤差を厳密に境界付けるために区間演算を用い、ヒューリスティックな誤差推定を置き換える。
• 非自明なゼータ関数の零点の和の切り捨て誤差に対して明示的な境界を適用し、リーマン＝フォン・マンゴルト関数 N(t) とテイル推定を用いる。
• レムマ A.1 を用いて、既知の漸近展開と誤差項を用いて N(t) を境界付け、零点和の切り捨てに対する精密な制御を可能にする。
• 零点が [T₁, T₂] 内にある場合と [−T₂, T₂] 外にある場合のそれぞれに E₁ と E₂ の別々の誤差境界を導入し、ステイエルジス積分と指数的減衰推定を用いる。

実験結果

リサーチクエスチョン

• RQ1 ラガリアス＝オドリツコの解析的手法を、区間演算と厳密な誤差境界を用いて実装することで、スケールに応じた条件なしの π(x) の計算が可能になるか？
• RQ2 高精度な解析的整数論計算において、数値的不安定性と浮動小数点誤差の蓄積をどのように制御できるか？
• RQ3 スムージングパラメータ λ と切り捨て閾値 T₁, T₂ の最適なバランスは何か？誤差を最小限に抑えつつ、計算の実行可能性を維持するには？
• RQ4 既知の零点と未検証の零点の両方を考慮した、非自明なゼータ関数の零点の和の切り捨てに対して、数学的に厳密な誤差境界を適用できるか？
• 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) の既知の値とのクロスバリデーション、および素数定理との整合性を満たしている。