[论文解读] A natural prime-generating recurrence

Eric Rowland|ArXiv.org|Oct 17, 2007

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

一句话总结

本文提出了一种自然发生的素数生成递推关系，定义为 $ a(n) = a(n-1) + \gcd(n, a(n-1)) $，初始值 $ a(1) = 7 $，并证明了差值 $ a(n) - a(n-1) $ 始终为 1 或素数。其主要贡献在于利用最大公约数的数论性质以及一个结构引理，给出了一个初等证明，该引理识别出满足 $ \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(n) = a(n-1) + \gcd(n, a(n-1)) $ 在 $ a(1) = 7 $ 条件下，其差值仅产生 1 和素数。
• 将该递推关系的行为推广至任意初始条件，并分析其长期演化与收敛性。
• 基于最大公约数行为的结构，开发一种计算捷径，以高效提取素数输出，而无需计算所有中间步骤。

提出的方法

• 定义递推关系 $ a(n) = a(n-1) + \gcd(n, a(n-1)) $，初始值 $ a(1) = 7 $，并分析差值序列 $ \gcd(n, a(n-1)) $。
• 利用基本性质：$ \gcd(n, m) $ 整除任意整数线性组合 $ rn + sm $，以分析非平凡最大公约数的结构。
• 引入一个引理，识别出满足 $ \gcd(n+k, a(n-1)+k) \neq 1 $ 的最小 $ k \geq 1 $，其依据是 $ \Delta = a(n-1) - n $ 的素因子。
• 将最小 $ k $ 表达为 $ \min\{ \textnormal{mod}_1(-n, p) \mid p \text{ 为素数}, p \mid \Delta \} $，其中 $ p $ 是 $ \Delta $ 的素因子。
• 应用该引理跳过差值序列中的一系列 1，直接计算下一个非平凡最大公约数，其值为 1 或素数。
• 通过变换 $ n' = n/2 $，$ a'(n') = a(n) - n/2 $，利用 $ r=3 $ 的演化模拟 $ r=2 $ 的演化，实现更快计算并支持结构比较。

实验结果

研究问题

• RQ1 递推关系 $ a(n) = a(n-1) + \gcd(n, a(n-1)) $ 在 $ a(1) = 7 $ 条件下，是否仅生成 1 和素数作为差值？
• RQ2 该递推关系的行为能否推广至任意初始条件？$ a(n)/n $ 的极限行为是什么？
• RQ3 1 的序列具有何种结构？能否高效跳过这些 1，以仅提取素数输出？
• RQ4 是否所有素数（可能除外 2）都会出现在该递推关系的差值序列中？
• RQ5 不同初始条件的演化路径如何？它们是否会最终收敛到等价的轨迹？

主要发现

• 递推关系 $ a(n) = a(n-1) + \gcd(n, a(n-1)) $ 在 $ a(1) = 7 $ 条件下，其差值仅由 1 或素数组成，该结论通过最大公约数的数论性质得到证明。
• 满足 $ \gcd(n+k, a(n-1)+k) \neq 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 的素数都会在经历大段间隔后被生成，且下一个非平凡最大公约数恒为 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 $ 的演化实现模拟，保持最大公约数结构，从而实现更快计算。