Pisano Period

参考自一篇全英文的issueThe Period of the Fibonacci Sequence Modulo j

\(\text{Definition}\)

\(\text{Fibonacci}\)数列

\(\forall n\in\mathbb{N_+},F_{n+1}=F_n+F_{n-1}\quad(F_0=0,F_1=1)\)

令\(\phi=\frac{1+\sqrt5}{2},\overline{\phi}=\frac{1-\sqrt5}{2}\)，易知有\(\phi^2=\phi+1\)与\(\overline{\phi}^2=\overline{\phi}+1\)。
实际上这就是\(\text{Fibonacci}\)数列的特征方程\(x^2=x+1\)的两实数根。
并且由此我们可以推出\(\text{Fibonacci}\)数列的通项公式\(F_n=\frac{1}{\sqrt5}(\phi^n-\overline{\phi}^n)\)。

\(\text{Pisano Period}\)

\(F_n\bmod P\)的循环节\(m\)是使得\(F_m\equiv0\pmod P\wedge F_{m+1}\equiv1\pmod P\)的最小正整数\(m\)。

易知\(m\mid k\Leftrightarrow F_k\equiv0\pmod P,F_{k+1}\equiv1\pmod P\)。

\(\text{Part.1}\)

我们需要先证明一些引理。

引理\(1\)

设\(p\in\mathbb P,n\in\mathbb{N_+}\)，\(a\equiv1\pmod p\Rightarrow a^{p^n}\equiv1\pmod{p^{n+1}}\)。

证：
我们使用数学归纳法。
不妨设\(a=rp+1\quad(r\in\mathbb Z)\)。
运用二项式定理，\(a^p\equiv(rp+1)^p\equiv1+\sum\limits_{i=1}^p\left(_i^p\right)(rp)^i\equiv1\pmod{p^2}\)。
因此当\(n=1\)时引理成立。
假设\(n=m\)时定理成立，即有\(a^{p^m}\equiv1\pmod{p^{m+1}}\)。
不妨设\(a^{p^m}=1+sp^{m+1}(s\in\mathbb Z)\)。
同样运用二项式定理，\(a^{p^{m+1}}=(a^{p^m})^p=(1+sp^{m+1})^p=1+\sum\limits_{i=1}^p\left(_i^p\right)(sp^{m+1})^i\equiv1\pmod{p^{m+2}}\)。
因此当\(n=m\)成立时，\(n=m+1\)也成立。
命题得证。

推论\(1\)

设\(p\in\mathbb P,k\in\mathbb{N_+}\)，\(m\)为\(F_n\bmod p\)的循环节，有\(\phi^{mp^{k-1}}\equiv\overline{\phi}^{mp^{k-1}}\equiv1\pmod{p^k}\)

证：
由\(F_m\equiv\frac{\phi^m-\overline\phi^m}{\sqrt5}\equiv0\pmod p\)可知\(\phi^m\equiv\overline{\phi}^m\pmod p\)。
我们有\(F_m=F_{m+1}-F_1=\frac{\phi(\phi^m-1)-\overline{\phi}(\overline{\phi}^m-1)}{\sqrt5}\)。
用\(\phi^m\)替换\(\overline{\phi}^m\)得到\(F_m\equiv\frac{(\phi^m-1)(\phi-\overline{\phi})}{\sqrt5}\equiv(\phi^m-1)F_1\equiv\phi^m-1\equiv0\pmod p\)。
\(\therefore\phi^m\equiv\overline{\phi}^m\equiv1\pmod p\)
运用引理\(1\)，\(\phi^{mp^{k-1}}\equiv\overline{\phi}^{mp^{k-1}}\equiv1\pmod{p^k}\)。
命题得证。

引理\(2\)

设\(p\in\mathbb P\wedge p\ne5\)，\((\frac5p)=1\Leftrightarrow p\equiv\pm1\pmod5,(\frac5p)=-1\Leftrightarrow p\equiv\pm2\pmod5\)

证：
运用Gauss二次互反律，\((\frac5p)(\frac p5)=(-1)^{(\frac{5-1}2)(\frac{p-1}2)}=((-1)^{\frac{p-1}2})^2=1\)。
运用Euler准则，\((\frac5p)=(\frac p5)\equiv p^2\pmod5\)。
利用枚举法发现该命题在任意情况下成立，故命题得证。

\(\text{Part.2}\)

定理\(1\)

设\(P=\prod\limits_{i=1}^{s}p_i^{k_i}\)，\(m_i\)为\(F_n\bmod p_i^{k_i}\)的循环节，\(M\)为\(F_n\bmod P\)的循环节，则\(M=\operatorname{lcm}(m_1,\cdots,m_s)\)。

证：
\(\because F_M\equiv0\pmod P\Leftrightarrow \forall i\in[1,s],F_M\equiv0\pmod{p_i^{k_i}}\)
\(\therefore\forall i\le s,m_i\mid M\)
\(\therefore M=\operatorname{lcm}(m_1,\cdots,m_s)\)
命题得证。

定理\(2\)

设\(p\in\mathbb P\)，\(m\)是\(F_n\bmod p\)的循环节，\(M\)是\(F_n\bmod p^k\)的循环节，则\(M\mid mp^{k-1}\)。

证：
由推论\(1\)我们有\(\phi^{mp^{k-1}}\equiv\overline{\phi}^{mp^{k-1}}\equiv1\pmod{p^k}\)。
\(\therefore F_{mp^{k-1}}\equiv\frac{\phi^{mp^{k-1}}-\overline\phi^{mp^{k-1}}}{\sqrt5}\equiv\pmod{p^k}\)
\(\therefore F_{mp^{k-1}+1}\equiv\frac{\phi^{mp^{k-1}+1}-\overline\phi^{mp^{k-1}+1}}{\sqrt5}\equiv\frac{\phi-\overline{\phi}}{\sqrt5}\equiv1\pmod{p^k}\)
\(\therefore M\mid mp^{k-1}\)
命题得证。

事实上有猜想\(M=mp^{k-1}\)一定成立，目前尚未找到反例。

定理\(3\)

设\(p\in\mathbb P\)，\(m\)为\(F_n\bmod p\)的循环节，若\(p\equiv\pm1\pmod5\)，则\(m\mid p-1\)

证：
由引理2知\((\frac 5p)=1\)，所以\(\sqrt5\)在\(\mathbb{F_p}\)中。
运用Euler定理，\(\phi^{p-1}\equiv\overline{\phi}^{p-1}\equiv 1\pmod p\)。
\(\therefore F_{p-1}\equiv\frac{\phi^{p-1}-\overline\phi^{p-1}}{\sqrt5}\equiv0\pmod p,F_p\equiv\frac{\phi^p-\overline\phi^p}{\sqrt5}\equiv\frac{\phi-\overline{\phi}}{\sqrt5}\equiv1\pmod p\)
\(\therefore m\mid p-1\)
命题得证。

定理\(4\)

设\(p\in\mathbb P\)，\(m\)为\(F_n\bmod p\)的循环节，若\(p\equiv\pm2\pmod5\)，则\(m\mid 2p+2\wedge 2\nmid\frac{2p+2}m\)

证：
由引理2知\((\frac5p)=-1\)，所以我们定义扩域\(\mathbb{F_p}(\sqrt5)=\{a+b\sqrt5|a,b\in\mathbb{F_p}\}\)。
利用二项式定理，我们有\((a+b)^p\equiv a^p+b^p\pmod p\)。
\(\therefore\phi^p\equiv(\frac12+\frac{\sqrt5}{2})^p\equiv(\frac12)^p+(\frac{\sqrt5}2)^p\equiv(\frac12)^p(1+\sqrt5^p)\equiv\frac12(1+5^\frac{p-1}2\sqrt5)\equiv\frac12(1-\sqrt5)\equiv\overline{\phi}\pmod p\)
同理，可以得到\(\overline{\phi}^p\equiv\phi\pmod p\)。
因此我们有
\(F_{p}\equiv\frac{\phi^p-\overline\phi^p}{\sqrt5}\equiv\frac{\overline\phi-\phi}{\sqrt5}\equiv p-1\pmod p\)
\(F_{p+1}\equiv\frac{\phi^{p+1}-\overline\phi^{p+1}}{\sqrt5}\equiv\frac{\overline\phi\phi-\phi\overline\phi}{\sqrt5}\equiv0\pmod p\)
\(F_{p+2}\equiv F_p+F_{p+1}\equiv p-1\pmod p\)
\(\therefore m\nmid p+1\)
同理有
\(F_{2p+1}\equiv\frac{\phi^{2p+1}-\overline\phi^{2p+1}}{\sqrt5}\equiv\frac{\overline\phi^2\phi-\phi^2\overline\phi}{\sqrt5}\equiv\frac{\phi\overline\phi(\overline\phi-\phi)}{\sqrt5}\equiv1\pmod p\)
\(F_{2p+2}\equiv\frac{\phi^{2p+2}-\overline\phi^{2p+2}}{\sqrt5}\equiv\frac{\overline\phi^2\phi^2-\phi^2\overline\phi^2}{\sqrt5}\equiv0\pmod p\)
\(F_{2p+3}\equiv F_{2p+1}+F_{2p+2}\equiv1\pmod p\)
\(\therefore m\mid 2p+2\)
命题得证。

定理5

若\(P=2\times5^k(k\in\mathbb{N_+})\)，则\(F_n\bmod P\)的循环节为\(6P\)，否则\(F_n\bmod P\)的循环节\(\le4P\)。

\(\text{Extra}\)

我们发现\(\text{Part.2}\)给出了Pisano Period的某个正整数倍数。
如果仅仅是求\(F_n\bmod m\)的话这就足够了。
如果需要求最小正周期的话，在运用定理\(3,4\)的时候枚举约数检查即可。

在计算\(F_n\bmod p\)的循环节时，需要对\(p=2\)和\(p=5\)进行特判：
\(m_i=\begin{cases}3&p_i=2\\20&p_i=5\end{cases}\)
最后还要记得特判\(P=1\)的情况。