[论文解读] The Fibonacci sequence modulo $p^2$ -- An investigation by computer for $p < 10^{14}$

Andreas-Stephan Elsenhans, Jörg Jahnel|arXiv (Cornell University)|Jun 4, 2010

Advanced Mathematical Theories and Applications参考文献 5被引用 10

一句话总结

本文通过大规模计算方法，研究了小于 $10^{14}$ 的素数 $p$ 的斐波那契序列模 $p^2$ 的周期长度。证明了对所有此类素数均有 $\kappa(p^2) \neq \kappa(p)$，从而确认标准公式 $\kappa(p^n) = \kappa(p) \cdot p^{n-1}$ 在 $p < 10^{14}$ 范围内普遍成立，解决了斐波那契模运算中长期存在的一个开放问题。

ABSTRACT

We show that for primes $p < 10^{14}$ the period length $\kappa (p^2)$ of the Fibonacci sequence modulo $p^2$ is never equal to its period length modulo $p$. The investigation involves an extensive search by computer. As an application, we establish the general formula $\kappa (p^n) = \kappa (p) \cdot p^{n-1}$ for all primes less than $10^{14}$.

研究动机与目标

• 本文旨在解决一个长期开放的问题：公式 $\kappa(p^n) = \kappa(p) \cdot p^{n-1}$ 是否对所有素数 $p$ 和 $n \geq 1$ 成立。
• 研究模 $p^2$ 的周期长度是否可能等于模 $p$ 的周期长度，这将与标准公式相矛盾。
• 研究重点在于验证斐波那契序列在素数幂模下的行为，特别是 $F_{\kappa(p)}$ 在模 $p^2$ 下的 $p$-进赋值。
• 本文提供了广泛素数范围内的计算证据和理论依据，以支持 $\kappa(p^n)$ 的通用公式（即 Wall 数）。

提出的方法

• 作者使用计算机搜索，对所有小于 $10^{14}$ 的素数 $p$ 计算了 $\kappa(p)$ 和 $\kappa(p^2)$，并利用优化的模斐波那契算术算法。
• 采用 Montgomery 表示法以加速大规模计算中的模算术运算。
• 关键的理论工具是基本引理 2.2.4，其内容为：若 $\kappa(p) \mid l$ 且 $p^e \mid F_l$，则 $\kappa(p^e) \mid l$，该结论使周期长度的高效测试成为可能。
• 算法通过检查 $F_{\kappa(p)} \equiv 0 \pmod{p^2}$ 来判断 $p$ 是否为异常素数（即 $\kappa(p^2) \ne \kappa(p)$）。
• 使用筛法排除合数候选，并通过斐波那契合数性测试 $p \mid F_{p-1}$ 验证素性。
• 对商 $R(p) = F_{p-1} \mod p^2 / p$（或 $F_{2p+2} \mod p^2 / p$）进行统计分析，以评估分布模式并检测异常现象。

实验结果

研究问题

• RQ1 是否存在某个素数 $p < 10^{14}$ 使得 $\kappa(p^2) = \kappa(p)$？
• RQ2 公式 $\kappa(p^n) = \kappa(p) \cdot p^{n-1}$ 是否对所有素数 $p < 10^{14}$ 成立？
• RQ3 对于满足 $p \equiv \pm1 \pmod{5}$ 的素数，商 $R(p) = (F_{p-1} \mod p^2)/p$ 的分布如何？
• RQ4 在满足 $p \equiv \pm1 \pmod{5}$ 的素数中，$Q(p) = 1$ 在商 $Q(p) = (p-1)/\kappa(p)$ 中出现的频率如何？
• RQ5 是否存在异常素数 $p$，使得 $p^2 \mid F_{\kappa(p)}$，即 $\kappa(p^2) \ne \kappa(p)$？

主要发现

• 对所有小于 $10^{14}$ 的素数 $p$，均有 $\kappa(p^2) \ne \kappa(p)$，即在此范围内不存在异常素数。
• 对所有小于 $10^{14}$ 的素数 $p$，通用公式 $\kappa(p^n) = \kappa(p) \cdot p^{n-1}$ 成立，确认其在极广范围内的有效性。
• 在小于 $2 \cdot 10^7$ 且满足 $p \equiv \pm2 \pmod{5}$ 的素数中，商 $Q(p) = 1$ 的频率约为 78.8%，表明 $Q(p) = 1$ 是常见现象。
• 在区间 $[0, 10^{12}]$ 内，有 30 个素数使得商 $R'(p)$ 为个位数，与预期值约 26.85 非常接近。
• 小商数 $R'(p)$ 的分布与理论预期高度一致：共有 12162 个商小于 10,000，216 个小于 100。
• 在区间 $[8 \cdot 10^{13}, 8.5765 \cdot 10^{13}]$ 内的 10,792 个素数中未发现任何合数，且所有素数均通过斐波那契测试及独立验证确认为素数。