[論文レビュー] A natural prime-generating recurrence

Eric Rowland|ArXiv.org|Oct 17, 2007

Advanced Combinatorial Mathematics参考文献 8被引用数 16

ひとこと要約

この論文は、初期値 $ a(1) = 7 $ を持つ自然に生じる素数生成再帰関係 $ a(n) = a(n-1) + \gcd(n, a(n-1)) $ を提示し、差分 $ a(n) - a(n-1) $ が常に 1 または素数であることを証明する。主な貢献は、GCD の数論的性質と、$ \gcd(n+k, a(n-1)+k) \neq 1 $ を満たす最小の $ k $ を特定する構造的補題を用いた、素数性の性質に関する初等的証明であり、非自明な素数出力を効率的に計算可能である。

ABSTRACT

For the sequence defined by a(n) = a(n-1) + gcd(n, a(n-1)) with a(1) = 7 we prove that a(n) - a(n-1) takes on only 1s and primes, making this recurrence a rare "naturally occurring" generator of primes. Toward a generalization of this result to an arbitrary initial condition, we also study the limiting behavior of a(n)/n and a transience property of the evolution.

研究の動機と目的

• 非設計的・自然に生じる再帰関係が、信頼性を持って素数を生成できるかどうかを調査すること。
• 初期値 $ a(1) = 7 $ を持つ再帰関係 $ a(n) = a(n-1) + \gcd(n, a(n-1)) $ が、差分として 1 または素数のみを生成することを証明すること。
• 任意の初期条件に対する再帰関係の挙動を一般化し、その長期的挙動と収束性を分析すること。
• GCD の挙動の構造に基づいた計算の最適化手法を開発し、中間ステップをすべて計算せずに素数出力を効率的に抽出すること。

提案手法

• 初期値 $ a(1) = 7 $ を持つ再帰関係 $ a(n) = a(n-1) + \gcd(n, a(n-1)) $ を定義し、差分 $ \gcd(n, a(n-1)) $ の列を分析する。
• 任意の整数線形結合 $ rn + sm $ を割り切るという基本的性質 $ \gcd(n, m) $ を用いて、非自明な GCD の構造を分析する。
• 初期値 $ \Delta = a(n-1) - n $ の素因数に基づき、$ \gcd(n+k, a(n-1)+k) \neq 1 $ を満たす最小の $ k \geq 1 $ を特定する補題を導入する。
• この最小の $ k $ を $ \min\{ \textnormal{mod}_1(-n, p) \mid p \text{ は素数}, p \mid \Delta \} $ として表現する。
• この補題を応用して、差分列における 1 の連続を飛ばし、次に現れる非自明な GCD（1 または素数）を直接計算する。
• 変換 $ n' = n/2 $、$ a'(n') = a(n) - n/2 $ を用いて、$ r=2 $ の進化を $ r=3 $ で模倣し、計算を高速化するとともに構造的比較を可能にする。

実験結果

リサーチクエスチョン

• RQ1 初期値 $ a(1) = 7 $ を持つ再帰関係 $ a(n) = a(n-1) + \gcd(n, a(n-1)) $ は、差分として 1 または素数のみを生成するか？
• RQ2 この再帰関係の挙動は任意の初期条件に一般化可能か？また、$ a(n)/n $ の極限挙動は何か？
• RQ3 1 の連続の列の構造は何か？また、それらを効率的に飛ばして、素数出力のみを抽出できるか？
• RQ4 すべての素数（2 を除く）は、この再帰関係の差分列に現れるか？
• RQ5 異なる初期条件ではどのように進化するか？また、それらは最終的に同値な軌道に収束するか？

主な発見

• 初期値 $ a(1) = 7 $ を持つ再帰関係 $ a(n) = a(n-1) + \gcd(n, a(n-1)) $ は、GCD の数論的性質を用いた証明により、差分が 1 または素数に限定されることを示した。
• 差分列における 1 の連続を飛ばすための最小の $ k \geq 1 $ は、$ \min\{ \textnormal{mod}_1(-n, p) \mid p \text{ は素数}, p \mid (a(n-1) - n) \} $ として与えられ、1 の飛躍が効率的に行える。
• 初期値 $ a(1) = 7 $ から 10,000 回の最適化を適用した後、まだ生成されていない最小の奇素数は 587 であり、2 を除くすべての素数が最終的に現れる可能性が高いことを示唆している。
• r=3 の場合、6 で割った余りが 5 である素数は、常に大きなギャップの後に生成され、次に現れる非自明な GCD は常に 3 である。
• 異なる初期条件からの進化は頻繁に同じ軌道に収束し、初期値 $ 4 \leq a(1) \leq 2^{13} $ に対して $ n < 2^{23} $ の範囲で、203 個の異なる同値類が得られた。
• $ r=2 $ の進化は、変換 $ n' = n/2 $、$ a'(n') = a(n) - n/2 $ を用いて $ r=3 $ の進化で模倣可能であり、GCD の構造を保ちつつ計算を高速化できる。