[論文レビュー] Benford's law for the $3x+1$ function

Jeffrey C. Lagarias, K. Soundararajan|ArXiv.org|Sep 7, 2005

Benford’s Law and Fraud Detection参考文献 12被引用数 4

ひとこと要約

本稿は、$3x+1$ 関数のほとんどの初期値に対して、最初の $N$ 回の反復が任意の基数 $B \geq 2$ において、$B$ が 2 のべきである場合を含め、ほぼベンフォードの法則を満たすことを示している。$3x+1$ 動的システムの確率的モデルと差異理論を用いて、ほとんどの初期値に対して、$ olimits\{\log_B x_k\}$ がモジュロ 1 で一様分布することを証明し、差異が $O(N^{-1/36})$ の速度で減少することを示した。これにより、$3x+1$ の反復に対する強い定量的形でのベンフォードの法則が確立された。

ABSTRACT

We show that for most choices of an initial seed $x_0$, the sequence of the first $N$ iterates of $x_0$ under the $3x+1$ map approximately satisfies Benford's law.

研究の動機と目的

• 大規模な $N$ に対して、$3x+1$ 関数の初期反復が基数 $B$ においてベンフォードの法則をほぼ満たすかどうかを調査すること。
• 初期値 $x_0$ のうち、列 $\{\log_B T^{(k)}(x_0)\}$ がモジュロ 1 で一様分布するものの割合を定量的に特定すること。
• 従来の研究で奇数番目の反復のみを扱ったモデルにおいて、2 のべきである基数 $B$ が異常とされたのを解消し、すべての基数 $B \geq 2$ に一般化したベンフォードの法則の結果を拡張すること。
• $3x+1$ の反復から得られる対数列の差異に対して、厳密かつ定量的な上限を提示すること。

提案手法

• $3x+1$ 動的システムを、各ステップで独立に $T(n) = \frac{3n}{2}$ または $T(n) = \frac{n}{2}$ を確率 $1/2$ で選択する確率的過程でモデル化し、$\log_B x_k$ が $\theta_1 = \log_B(3/2)$ または $\theta_2 = \log_B(1/2)$ だけ増加するようにする。
• $y_k = \log_B x_k$ の列 $\{y_1, \dots, y_N\}$ の差異 $D(\{y_1, \dots, y_N\})$ を、モジュロ 1 における一様分布からの逸脱度として定義する。
• ディオファントス近似に関する結果（指数 $\alpha = 7.616$）を用いて、確率的過程 $\mathcal{P}(\theta_1, \theta_2, y_0 = 0)$ の期待値としての差異が有界であることを用いる。
• マルコフの不等式を用いて、差異が $N^{-1/36}$ を超える確率を評価し、例外的初期値集合の定量的上限を導出する。
• 持ち上げ補題を用いて、実際の $3x+1$ の反復と、確率的過程の切り詰め版との間に接続を確立し、ほとんどの初期値に対して差異が保存されることを示す。
• 周期性と区間被覆の議論を用いて、$X = 2^N$ から一般の $X > 2^N$ への結果の拡張を実現する。

実験結果

リサーチクエスチョン

• RQ1 ほとんどの初期値に対して、$3x+1$ 関数の最初の $N$ 回の反復は、基数 $B$ においてベンフォードの法則をほぼ満たすか？
• RQ2 $3x+1$ の反復から得られる $\{\log_B x_k\}$ の差異に対して、定量的な上限を確立できるか？
• RQ3 奇数番目の反復のみを扱ったモデルで以前に異常とされたが、$B$ が 2 のべきである場合を含め、すべての基数 $B \geq 2$ に対してこの結果が成り立つか？
• RQ4 差異は $N$ が増加するにつれてどのように減少するか？ また、ベンフォードの法則に背く「例外的」な初期値の割合はどの程度か？

主な発見

• $X \geq 2^N$ である $1 \leq x_0 \leq X$ のほとんどの初期値 $x_0$ に対して、列 $\{\log_B x_k\}_{k=1}^N$ の差異は $D \leq N^{-1/36} + 2^{2 - \frac{1}{100}N}$ を満たす。
• 差異が $N^{-1/36}$ を超える初期値 $x_0$ の例外的集合の濃度は、ある定数 $C$ を用いて、$C N^{-1/36} 2^N + 2^{1 + \frac{99}{100}N}$ 以下である。
• 差異の上限は、$B$ が 2 のべきである場合を含め、すべての基数 $B \geq 2$ に対して一様に成り立ち、従来の奇数反復モデルで生じていた異常を解消した。
• $3x+1$ の反復をモデル化する確率的過程の期待差異は $O(N^{-1/18})$ で有界であり、マーティンの不等式により、これが主結果を示唆する。
• 本結果は、実際の $3x+1$ の反復の差異が、ディオファントス近似によってよく理解されている関連する確率的過程の差異に非常に近いことを示すことで確立された。