斐波那契循环节的一种求法

原链接

---

题意

求$f(f(n))\mathrm{m}\mathrm{o}\mathrm{d}p,n\le 10^{100}$

---

循环节问题

设$g(n)$表示$\mathrm{m}\mathrm{o}\mathrm{d}n$意义下斐波那契数列的循环节长度。

设模数$n=\prod _{i=1}^m{p}_i^{k_i}$，则：
$$g(n)=lcm(g(p_1^{k_1}),g(p_2^{k_2}),...,g(p_m^{k_m}))$$
且
$$g(p^k)=g(p)\ast p^{k-1}$$
($g$均指最小循环节长度)

对于$g(p)$又有以下定理：

如果5是模$p$的二次剩余那么循环节的长度是$p-1$的因子否则长度为$2(p+1)$
对于小于等于5的素数特殊判断，$g(2)=3,g(3)=8,g(5)=20$。

(证明先咕了)
所以算出两个二次剩余之后矩乘即可。

---

代码

代码也咕了。