Pisano Period

参考自一篇全英文的paperThe 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\)的情况。

posted @ 2019-08-20 20:10  Shiina_Mashiro  阅读(2698)  评论(3编辑  收藏