[论文解读] 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$下（包括2的幂）近似满足本福德定律。通过建立$3x+1$动力系统的随机模型并结合分歧理论，证明了对于绝大多数初值，对数序列$\{\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.

研究动机与目标

• 研究$3x+1$函数的初始迭代在大$N$时是否在基$B$下近似满足本福德定律。
• 量化满足序列$\{\log_B T^{(k)}(x_0)\}$模1均匀分布的初值$x_0$的比例。
• 将先前关于$3x+1$本福德定律的结果扩展至所有$B \geq 2$的基，包括此前在仅考虑奇数次迭代的模型中表现异常的2的幂。
• 为$3x+1$迭代导出的对数序列的分歧度提供严格且定量的上界。

提出的方法

• 使用一个随机过程对$3x+1$动力系统建模，其中每一步独立地以相等概率应用$T(n) = \frac{3n}{2}$或$T(n) = \frac{n}{2}$，对应于$\log_B x_k$增加$\theta_1 = \log_B(3/2)$或$\theta_2 = \log_B(1/2)$。
• 将序列$y_k = \log_B x_k$的分歧度$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\}$模1的分歧度建立定量上界？
• RQ3 该结果是否对所有$B \geq 2$的基（包括2的幂）均成立，即使在仅考虑奇数次迭代的模型中曾出现异常？
• RQ4 随着$N$增大，分歧度如何衰减？有多少比例的初值属于‘异常’集合，即违反本福德定律？

主要发现

• 对于满足$1 \leq x_0 \leq X$且$X \geq 2^N$的绝大多数初值$x_0$，序列$\{\log_B x_k\}_{k=1}^N$的分歧度满足$D \leq N^{-1/36} + 2^{2 - \frac{1}{100}N}$。
• 分歧度超过$N^{-1/36}$的初值异常集合的基数至多为$C N^{-1/36} 2^N + 2^{1 + \frac{99}{100}N}$，其中$C$为某常数。
• 分歧度上界对所有$B \geq 2$一致成立，包括$B$为2的幂的情况，解决了此前仅考虑奇数次迭代模型中的异常问题。
• 建模$3x+1$迭代的随机过程的期望分歧度有界于$O(N^{-1/18})$，结合马尔可夫不等式可推出主结论。
• 该结果通过证明实际$3x+1$迭代的分歧度与一个相关随机过程的分歧度极为接近而建立，而后者可通过丢番图逼近理论得到充分理解。