解题报告 (八) 矩阵快速幂【进阶】

一、矩阵快速幂 算法讲解

夜深人静写算法（二十）矩阵快速幂

二、矩阵快速幂 解题报告

HDU 3658 How many words

• 链接：HDU 3658 How many words
• 题意：一个字符串只能包含大小写字母，必须满足两条规则：
  1）任意两个相邻字符的 $ASCII$ 码差值必须小于等于 $32$；
  2）至少一对相邻字符的 $ASCII$ 码 差值正好等于 $32$；
• 给定长度 $m(2<=m<=10^9)$，问总共有多少不同字符；
• 难度：★★★☆☆
• 题解：用 X 表示任意相邻字符差值小于等于 32 的数量；Y表示任意相邻字符差值小于32的数量，Z表示任意相邻字符差值不等于32的数量；因为条件 1，我们必须取集合 X ；因为条件 2，我们必须取集合 Z 的非，所以画出如下图：
• X是小于 或 等于，Z是小于或大于，所以他们的交 Y 就是小于了，最后的答案就是 X - Y；
  
• 建立邻接矩阵，矩阵的行代表前一个字符，列代表当前字符，如果行列字符的 ASCII 码 差值小于等于 32，则为1，否则为0；
• 建立二分矩阵计算 X， 然后用同样方法二分计算 Y，相减就是答案；

HDU 6030 Happy Necklace

链接：HDU 6030 Happy Necklace
题意：一个字符串由蓝色字符 (A表示) 和红色字符 (a表示) 组成，要求任意素数长度的子串中蓝色字符不能少于红色字符，给定字符串长度 $n<=10^{18}$，求方案数 $\%(10^9+7)$；
难度：★★★☆☆
题解：

• 1）考虑如果一个字符串以a结尾，那么后面只能跟A，因为2是素数，不能出现aa，所以有状态转移如下：

• 2）如果一个字符串以 aA 结尾，那么后面只能跟 A，因为 3是素数，不能出现 aAa，所以有状态转移如下：

• 3）如果一个字符串以 aAA 结尾，那么后面可以跟 a，也可以跟 A，原因是以 aAA 结尾的字符串要么就是 aAA 要么一定也是以AaAA 结尾的（因为不可能出现 aaAA），所以有状态转移；

• 4）最后来看 aAAA 和 AAA 的状态其实是一样的，所以可以进行状态合并，并且可以得出 AAA 的状态转移如下：

• 5）最后将所有状态组合在一个图上，得到一个状态转移图如下：

• 6）用数字对几个状态进行编号，得到状态转移矩阵如下：
  $$A=\left[\begin{array}{cccc}0& 1& 0& 0\\ 0& 0& 1& 0\\ 1& 0& 0& 1\\ 1& 0& 0& 1\end{array}\right]$$
• 长度为 $n$ 的方案数为 $A^{n+2}$ 第一行的和；

HDU 1588 Gauss Fibonacci

• 链接：HDU 1588 Gauss Fibonacci

• 题意：

• 1）$f(0)=0;f(1)=1;f(n)=f(n-1)+f(n-2)(n>=2)$

• 2）$g(i)=k\ast i+b$

• 3）给定 $k,b,n,M(0<=k,b,n<10^9)$，求 ${\sum }_i^{n-1}f(g(i))\%M$

• 难度：★★★☆☆

• 题解：一维递推类问题，构造系数矩阵如下：

• 我们知道斐波那契数列的系数矩阵如下：
  $$\begin{array}{rl}\left[\begin{array}{c}f(n)\\ f(n-1)\end{array}\right]& =\left[\begin{array}{cc}1& 1\\ 1& 0\end{array}\right]\left[\begin{array}{c}f(n-1)\\ f(n-2)\end{array}\right]\\ & ={\left[\begin{array}{cc}1& 1\\ 1& 0\end{array}\right]}^{n-1}\left[\begin{array}{c}f(1)\\ f(0)\end{array}\right]\\ & =A^{n-1}\left[\begin{array}{c}1\\ 0\end{array}\right]\end{array}$$

• 所以求 $f(g(i))=f(k\ast i+b)$，就是求 $A^{k\ast i+b-1}$，那么求和 ${\sum }_{i=0}^{n-1}f(g(i))$ 就是求
  $$\begin{array}{rl}T=& A^{b-1}+A^{k+b-1}+A^{2k+b-1}+...+A^{(n-1)k+b-1}\\ & =A^{b-1}\ast (A^0+A^k+A^{2k}+...+A^{(n-1)k})\end{array}$$
  令 $S(nk)=A^k+A^{2k}+...+A^{nk}$，则可以构造求和矩阵如下：
  $$\begin{array}{rl}\left[\begin{array}{c}S(nk)\\ I\end{array}\right]& =\left[\begin{array}{cc}{A}^k& A^k\\ O& I\end{array}\right]\left[\begin{array}{c}S((n-1)k)\\ I\end{array}\right]\\ & ={\left[\begin{array}{cc}{A}^k& A^k\\ O& I\end{array}\right]}^{n-1}\left[\begin{array}{c}S(k)\\ I\end{array}\right]\\ & =S^{n-1}\left[\begin{array}{c}{A}^k\\ I\end{array}\right]\end{array}$$

• 所以，依次求出 $A^k$、$S$、$S^{n-2}$、$S((n-1)k)$、$A^{b-1}$、$T$、$f(g(i))$；

• 注意，这题的特殊判断比较多，比如 $n=0、1、2；b=0；$ 等情况都需要进行特殊处理；

HDU 3936 FIB Query

• 链接：HDU 3936 FIB Query

• 题意：

• 1）$f(0)=0;f(1)=1;f(n)=f(n-1)+f(n-2)(n>=2)$

• 2）$g(i)=4\ast i-1$

• 3）给定 $L,R(1<=L,R<10^{12})$，求 ${\sum }_L^{R}f(g(i))\%1000000007$

• 难度：★★★☆☆

• 题解：首先，部分和可以转化成求前缀和，如下：
  $$\sum _L^{R}f(g(i))\%1000000007=(\sum _1^{R}f(g(i))-\sum _1^{L-1}f(g(i)))\%1000000007$$

• 问题转化成求 ${\sum }_{i=1}^{n}f(g(i))\%1000000007$；

• 一维递推类问题，构造系数矩阵如下：

• 我们知道斐波那契数列的系数矩阵如下：
  $$\begin{array}{rl}\left[\begin{array}{c}f(n)\\ f(n-1)\end{array}\right]& =\left[\begin{array}{cc}1& 1\\ 1& 0\end{array}\right]\left[\begin{array}{c}f(n-1)\\ f(n-2)\end{array}\right]\\ & ={\left[\begin{array}{cc}1& 1\\ 1& 0\end{array}\right]}^{n-1}\left[\begin{array}{c}f(1)\\ f(0)\end{array}\right]\\ & =A^{n-1}\left[\begin{array}{c}1\\ 0\end{array}\right]\end{array}$$

• 所以求 $f(g(i))=f(4\ast i-1)$，就是求 $A^{4\ast i-2}$，那么求和 ${\sum }_{i=1}^{n}f(g(i))$ 就是求
  $$\begin{array}{rl}T& =A^{4\ast 1-2}+A^{4\ast 2-2}+A^{4\ast 3-2}+...+A^{4\ast n-2}\\ & =A^2\ast (A^{4\ast 0}+A^{4\ast 1}+A^{4\ast 2}+...+A^{4\ast (n-1)})\end{array}$$
  令 $S(nk)=A^k+A^{2k}+...+A^{nk}$，则可以构造求和矩阵如下：
  $$\begin{array}{rl}\left[\begin{array}{c}S(nk)\\ I\end{array}\right]& =\left[\begin{array}{cc}{A}^k& A^k\\ O& I\end{array}\right]\left[\begin{array}{c}S((n-1)k)\\ I\end{array}\right]\\ & ={\left[\begin{array}{cc}{A}^k& A^k\\ O& I\end{array}\right]}^{n-1}\left[\begin{array}{c}S(k)\\ I\end{array}\right]\\ & =S^{n-1}\left[\begin{array}{c}{A}^k\\ I\end{array}\right]\end{array}$$

• 所以，依次求出 $A^4$、$S$、$S^{n-2}$、$S((n-1)k)$、$A^2$、$T$、$f(g(i))$；

• 注意，这题的特殊判断比较多，比如 $n=0、1、2；b=0；$ 等情况都需要进行特殊处理；

HDU 2254 奥运

• 链接：HDU 2254 奥运
• 题意：给定一个有向图，最多30个点，求v1到v2在第t1到t2天的方案数 (0<=t1, t2<10000)；
• 难度：★★★☆☆
• 题解：令 $f(n,v)$ 为从任意点出发，经过 n 步以后到达 v 的方案数，那么有状态转移如下：
  $$f(n,v)=\sum ^{u->v}f(n-1,u)$$
• 举个例子，有向图如下：

• 这个图的邻接矩阵如下（$A[u][v]表示u到v有多少条边$）：
  $$A=\left[\begin{array}{ccc}0& 1& 1\\ 1& 0& 0\\ 0& 0& 0\end{array}\right]$$
• 根据状态转移构造出如下矩阵递推关系：
  $$\left[\begin{array}{c}f(n,0)\\ f(n,1)\\ f(n,2)\end{array}\right]=\left[\begin{array}{ccc}0& 1& 0\\ 1& 0& 0\\ 1& 0& 0\end{array}\right]\left[\begin{array}{c}f(n-1,0)\\ f(n-1,1)\\ f(n-1,2)\end{array}\right]$$
• 观察发现 系数矩阵 和 正好是 邻接矩阵 的 转置，于是根据邻接矩阵建立反图；令系数矩阵为A，顶点数为 m，从 $v_1点到v_2$ 点经过 $n$ 天的方案数，可以构造出如下等式：
  $$\left[\begin{array}{c}f(n,0)\\ f(n,1)\\ ...\\ f(n,v_2)\\ ...\\ f(n,m-1)\end{array}\right]=A^{n-1}\left[\begin{array}{c}f(1,0)\\ f(1,1)\\ ...\\ f(1,v_1)\\ ...\\ f(1,m-1)\end{array}\right]$$
• 其中 $$f(1,x)=\{\begin{array}{ll}0& x<>v_1\\ 1& x=v_1\end{array}$$
• 这题求的是 $t_1$ 到 $t_2$ 天的方案数总和，所以需要求 $$A^{t_1}+A^{t_1+1}+...+A^{t_2}$$
• 令 $S(n)=A^1+A^2+...+A^n$，我们要求的就是 $S(t_2)-S(t_1-1)$；于是问题转换成求矩阵 n 次幂的前缀和；
• 构造一个阶数是 A 的两倍的 矩阵 $A^{\prime }$，且满足如下：
  $$\begin{array}{rl}\left[\begin{array}{c}S(n)\\ I\end{array}\right]& =\left[\begin{array}{cc}A& A\\ O& I\end{array}\right]\left[\begin{array}{c}S(n-1)\\ I\end{array}\right]\\ & ={\left[\begin{array}{cc}A& A\\ O& I\end{array}\right]}^n\left[\begin{array}{c}O\\ I\end{array}\right]\\ & ={A^{\prime }}^n\left[\begin{array}{c}O\\ I\end{array}\right]\end{array}$$
• 于是问题转化成求 矩阵 $A^{\prime }$ 的 n 次幂；

HDU 5607 graph

• 链接：HDU 5607 graph
• 题意：给定一张有向图，从任意点出发的概率都是均等的，问从起点 u 到 所有点 经过 K 步的概率；因为输出的概率是个小数，所以如果答案是 $\frac{X}{Y}$，则需要输出 $X\times Y^{10^9+5}\mathrm{m}\mathrm{o}\mathrm{d}10^9+7$.
• 难度：★★★☆☆
• 题解：直接举例更加直观；对于如下的有向图：

• 建立矩阵如下：
  $$\begin{array}{r}A=\left[\begin{array}{ccc}0& \frac{1}{3}& 0\\ \frac{1}{2}& \frac{1}{3}& 1\\ \frac{1}{2}& \frac{1}{3}& 0\end{array}\right]\end{array}$$
• 其中$A[i][j]$表示的是 $j$ 到 $i$ 的概率，所以如果 $j$ 到 $i$ 有边，那么令 $j$ 的出度为 $D$，$A[i][j]=\frac{1}{D}$；且 $A$ 这个矩阵的每一列的和恒为1；
• 由于这个矩阵存在小数，可以用 $D$ 对 $10^9+7$ 求逆元化成整数，得到新的矩阵如下：
  $$\begin{array}{r}A=\left[\begin{array}{ccc}0& 333333336& 0\\ 500000004& 333333336& 1\\ 500000004& 333333336& 0\end{array}\right]\end{array}$$
• 最后 $X$ 代表 $u$ 到某个点的概率值，$Y$ 代表 $u$ 到所有点的概率值，利用矩阵二分求解后按照题意中的进行二分幂即可；

HDU 3292 No more tricks, Mr Nanguo

• 链接：HDU 3292 No more tricks, Mr Nanguo
• 题意：等于等式 $X^2=N\times Y^2+1（X，Y>0）$，给定 $N(N<=29)、K$，求第 $K$ 大的 $X$ 的值；
• 难度：★★★☆☆
• 题解：首先，当 $N$ 为完全平方数的时候，则无解；否则，利用枚举得到最小的解 $X_1,Y_1$，则有如下等式成立：
  $$X_k+\sqrt{N}{Y}_k=(X_1+\sqrt{N}{Y}_1{)}^k$$，转化成递推的形式如下：
  $$\begin{array}{rl}{X}_k+\sqrt{N}{Y}_k& =(X_{k-1}+\sqrt{N}{Y}_{k-1})(X_1+\sqrt{N}{Y}_1)\\ & =(X_1{X}_{k-1}+N{Y}_1{Y}_{k-1})+\sqrt{N}(X_1{Y}_{k-1}+Y_1{X}_{k-1})\end{array}$$
• 得到递推矩阵如下：
  $$\begin{array}{rl}\left[\begin{array}{c}{X}_k\\ Y_k\end{array}\right]& =\left[\begin{array}{cc}{X}_1& N{Y}_1\\ Y_1& X_1\end{array}\right]\left[\begin{array}{c}{X}_{k-1}\\ Y_{k-1}\end{array}\right]\\ & =A^{k-1}\left[\begin{array}{c}{X}_1\\ Y_1\end{array}\right]\end{array}$$
• 矩阵二分求解 A 即可；

HDU 2256 Problem of Precision

• 链接：HDU 2256 Problem of Precision
• 题意：求 $\lfloor (\sqrt{2}+\sqrt{3}{)}^n\rfloor \mathrm{m}\mathrm{o}\mathrm{d}1024$
• 难度：★★★☆☆
• 题解：首先令$f(n)={(\sqrt{2}+\sqrt{3})}^{2n}=(5+2\sqrt{6}{)}^n$，无论 n 多大，一定可以表示成 $f(n)=A_n+B_n\sqrt{6}$，则有如下递推式：
  $$\begin{array}{rl}f(n)& =A_n+B_n\sqrt{6}\\ & =(A_{n-1}+B_{n-1}\sqrt{6})(5+2\sqrt{6})\\ & =(5A_{n-1}+12B_{n-1})+(2A_{n-1}+5B_{n-1})\sqrt{6}\end{array}$$
• 表示成矩阵的形式如下：
  $$\left[\begin{array}{c}{A}_n\\ B_n\end{array}\right]=\left[\begin{array}{cc}5& 12\\ 2& 5\end{array}\right]\left[\begin{array}{c}{A}_{n-1}\\ B_{n-1}\end{array}\right]={\left[\begin{array}{cc}5& 12\\ 2& 5\end{array}\right]}^{n-1}\left[\begin{array}{c}5\\ 2\end{array}\right]$$
  又由于 $(5+2\sqrt{6})(5-2\sqrt{6})=1$，得到 $(5+2\sqrt{6}{)}^n(5-2\sqrt{6}{)}^n=1$，则 $(A_n+B_n\sqrt{6})(A_n-B_n\sqrt{6})=1$，得到：
  $$6B_n^2=A_n^2-1$$，从而推导出$$B_n\sqrt{6}=\sqrt{A_n^2-1}$$，于是 $B_n\sqrt{6}$ 的整数部分就是 $A_n-1$，所以我们要求的答案就是 $(2A_n-1)\%1024$，矩阵求解 $A_n$；

HDU 4549 M斐波那契数列

链接：HDU 4549 M斐波那契数列
题意：$$F(0)=a,F(1)=b,F(n)=F(n-1)\ast F(n-2)(n>1)$$

• 给定 $a、b、n(n<=10^9)$，求 $F(n)\mathrm{m}\mathrm{o}\mathrm{d}1000000007$

• 难度：★★★☆☆

• 题解：扩展欧拉定理 + 矩阵二分

• 根据递推公式容易得到 $F(n)=a^{f(n-2)}{b}^{f(n-1)}$，其中 $f(0)=f(1)=1$，且满足 $f(n)=f(n-1)+f(n-2)$，即斐波那契数列；

• 令 $G(n,k)=k^{f(n)}\mathrm{m}\mathrm{o}\mathrm{d}1000000007$，那么我们只需要求出 $G(n,k)$ 就可以求得原式的解；

• 扩展欧拉定理
  $$a^b\mathrm{m}\mathrm{o}\mathrm{d}c=\{\begin{array}{ll}{a}^{b\mathrm{m}\mathrm{o}\mathrm{d}\phi (c)}\mathrm{m}\mathrm{o}\mathrm{d}c& gcd(a,c)=1\\ a^b\mathrm{m}\mathrm{o}\mathrm{d}c& gcd(a,c)\ne 1，b<\phi (c)\\ a^{b\mathrm{m}\mathrm{o}\mathrm{d}\phi (c)+\phi (c)}\mathrm{m}\mathrm{o}\mathrm{d}c& gcd(a,c)\ne 1，b>=\phi (c)\end{array}$$

• 当 $f(n)<\phi (1000000007)$ 说明 n 比较小，直接求递推求出来即可；

• 当 $f(n)>=\phi (1000000007)$ 时，根据扩展欧拉定理，得到：
  $$\begin{array}{rl}G(n,k)& =k^{f(n)}\mathrm{m}\mathrm{o}\mathrm{d}1000000007\\ & =k^{f(n)\mathrm{m}\mathrm{o}\mathrm{d}\phi (1000000007)+\phi (1000000007)}\mathrm{m}\mathrm{o}\mathrm{d}1000000007\end{array}$$

• $f(n)$ 是最基本的斐波那契数列，直接用矩阵二分求解即可；

HDU 5895 Mathematician QSC

• 链接：HDU 5895 Mathematician QSC
• 题意：已知公式：$$f(n)=\{\begin{array}{ll}0& (n==0)\\ 1& (n==1)\\ 2\ast f(n-1)+f(n-2)& (n>=2)\end{array}$$

$$g(n)=\sum _{i=0}^{n}f(i{)}^2$$
给出 $n、y、x、s$ 求：
$$x^{g(n\ast y)}\%(s+1)$$

• 难度：★★★☆☆
• 题解：扩展欧拉定理 + 矩阵二分
• f(n) 的矩阵比较容易构造，如下：
  $$\left[\begin{array}{c}f(n)\\ f(n-1)\end{array}\right]=\left[\begin{array}{cc}2& 1\\ 1& 0\end{array}\right]\left[\begin{array}{c}f(n-1)\\ f(n-2)\end{array}\right]$$
  那么如何构造 $f(n{)}^2$ 呢？我们尝试将原式带入，得到：
  $$f(n{)}^2=4f(n-1{)}^2+4f(n-1)f(n-2)+f(n-2{)}^2$$，所以递推式里面需要引入当前项和前一项的乘积；矩阵构造的时候，列向量应该是:
  $$\left[\begin{array}{c}f(n{)}^2\\ f(n)f(n-1)\\ f(n-1{)}^2\end{array}\right]$$
  进行递推矩阵的构造，得到如下矩阵：
  $$\left[\begin{array}{c}f(n{)}^2\\ f(n)f(n-1)\\ f(n-1{)}^2\end{array}\right]=\left[\begin{array}{ccc}4& 4& 1\\ 2& 1& 0\\ 1& 0& 0\end{array}\right]\left[\begin{array}{c}f(n-1{)}^2\\ f(n-1)f(n-2)\\ f(n-2{)}^2\end{array}\right]$$
  然而，这里又要求是求和，所以需要引用一个求和项，即
  $$\left[\begin{array}{c}sum(n)\\ f(n{)}^2\\ f(n)f(n-1)\\ f(n-1{)}^2\end{array}\right]$$
  整理后，得到最终的递推矩阵为：
  $$\left[\begin{array}{c}sum(n)\\ f(n{)}^2\\ f(n)f(n-1)\\ f(n-1{)}^2\end{array}\right]=\left[\begin{array}{cccc}1& 4& 4& 1\\ 0& 4& 4& 1\\ 0& 2& 1& 0\\ 0& 1& 0& 0\end{array}\right]\left[\begin{array}{c}sum(n-1)\\ f(n-1{)}^2\\ f(n-1)f(n-2)\\ f(n-2{)}^2\end{array}\right]$$
  即：
  $$\left[\begin{array}{c}sum(n)\\ f(n{)}^2\\ f(n)f(n-1)\\ f(n-1{)}^2\end{array}\right]={\left[\begin{array}{cccc}1& 4& 4& 1\\ 0& 4& 4& 1\\ 0& 2& 1& 0\\ 0& 1& 0& 0\end{array}\right]}^{n-1}\left[\begin{array}{c}1\\ 1\\ 0\\ 0\end{array}\right]$$
  根据扩展欧拉定理，有
  $$a^b\mathrm{m}\mathrm{o}\mathrm{d}c=\{\begin{array}{ll}{a}^{b\mathrm{m}\mathrm{o}\mathrm{d}\phi (c)}\mathrm{m}\mathrm{o}\mathrm{d}c& gcd(a,c)=1\\ a^b\mathrm{m}\mathrm{o}\mathrm{d}c& gcd(a,c)\ne 1，b<\phi (c)\\ a^{b\mathrm{m}\mathrm{o}\mathrm{d}\phi (c)+\phi (c)}\mathrm{m}\mathrm{o}\mathrm{d}c& gcd(a,c)\ne 1，b>=\phi (c)\end{array}$$
  然后直接套公式，矩阵二分求解；

HDU 2238 机器人的舞蹈II

• 链接：HDU 2238 机器人的舞蹈II

• 题意：一天四个相同的机器人在一张跳舞毯上跳舞，这是一张特殊的跳舞毯，他由4个正方形毯子组成一个大的正方形毯子，一开始四个机器人分别站在4块毯子上，舞蹈的每一步机器人可以往临近（两个毯子拥有同一条边视为临近）的一个毯子移动或停留在原来的毯子（同一块毯子可以有多个机器人停留），这个时候机器人的制造者你想知道经过n步的移动有多少种方式可以让每个毯子上都有机器人停留。
  

• 难度：★★★☆☆

• 题解：矩阵 + DFS

• 比较直观的理解就是每个机器人在每一步都有四个毯子选择，机器人一共有四个，所以任意情况总共有 $4^4$ 种状态，用一个五维数组表示状态，第一维表示经过到哪一步，后面四维分别代表每个机器人在那个步数下处于哪个毯子上；
  $$dp[step][s_0][s_1][s_2][s_3](其中0<s_0,s_1,s_2,s_3<4)$$，列出状态转移方程，如下：
  $$dp[step][s_0][s_1][s_2][s_3]=\sum dp[step-1][p_0][p_1][p_2][p_3]（p_i和s_i相邻）$$
  将后面四维状态编码到 0 - 255 的整数中，就可以列出一个 $256\times 256$ 的矩阵， 但是这样计算的话，时间复杂度是$O(log(n)N^3)$，所以需要做状态压缩；由于每个机器人都是一样的，所以这里其实有很多状态是可以合并的，考虑：
  1）将 每个机器人在哪个毯子上 的状态表示改成 每个毯子上有几个机器人；
  2）只有当四个数字加起来等于 4 的才是合法状态，状态数压缩了一大半，最后只有 35 个状态；
  3）对所有状态进行编码，使之能够连续，方便把稀疏矩阵转为稠密矩阵；

• 对每个状态枚举它能够到达的下一个状态，枚举方式采用深度优先搜索，比如当前状态为 $$\left[\begin{array}{cc}2& 1\\ 1& 0\end{array}\right]$$

• 枚举每个机器人能够去的位置（三个），然后得到下一个状态，这里需要注意的是，因为每个机器人都是一样的，所以对于2个机器人的毯子，其实就是将它分拆成三个非负整数的方案，即：
  $$\begin{array}{rl}2& =0+0+2\\ & =0+1+1\\ & =0+2+0\\ & =1+0+1\\ & =1+1+0\\ & =2+0+0\end{array}$$

• 通过深搜枚举出状态后构造状态转移矩阵，然后进行矩阵二分，再乘上一个 列向量B 就是答案， 列向量B 只有一个数字为 1（四个位置都是1的那个状态），其他均为0；

HDU 2855 Fibonacci Check-up

• 链接：HDU 2855 Fibonacci Check-up
• 题意：$$S(n)=\sum _{k=0}^n{C}_n^{k}F(k)$$其中 $F(k)$ 为斐波那契数列: $F(0)=0,F(1)=1$，给定 n 和 m，求 $S(n)\%m$ 的值；
• 难度：★★★★☆（推荐）
• 题解：二项式定理 + 矩阵二分
• 对于任意整数 $x$ 满足二项式定理如下：
  $$\begin{array}{r}(1+x{)}^n=C_n^0+C_n^1{x}^1+C_n^2{x}^2+...C_n^n{x}^n\end{array}$$
• 将 x 看作是一个矩阵，这个定理依然成立，即：
  $$\begin{array}{r}(I+A{)}^n=C_n^0+C_n^1{A}^1+C_n^2{A}^2+...C_n^n{A}^n\end{array}$$
• 其中 $I$ 为单位矩阵，我们知道斐波那契数列的系数矩阵的为 $\begin{array}{r}\left[\begin{array}{cc}1& 1\\ 1& 0\end{array}\right]\end{array}$，并且满足 $${\begin{array}{r}\left[\begin{array}{cc}1& 1\\ 1& 0\end{array}\right]\end{array}}^n=\begin{array}{r}\left[\begin{array}{cc}F(n)& F(n-1)\\ F(n-1)& F(n-2)\end{array}\right]\end{array}$$
• 令 $A=\begin{array}{r}\left[\begin{array}{cc}1& 1\\ 1& 0\end{array}\right]\end{array}，I=\begin{array}{r}\left[\begin{array}{cc}1& 0\\ 0& 1\end{array}\right]\end{array}$，套入二项式定理，得到：
  $$\begin{array}{rl}(I+A{)}^n& =C_n^0+C_n^1{A}^1+C_n^2{A}^2+...C_n^n{A}^n\\ & ={\left[\begin{array}{cc}2& 1\\ 1& 1\end{array}\right]}^n\\ & =B^n\end{array}$$
• 这时我们发现，$S(n)$ 就是 矩阵 $B^n$ 的 左上角那个元素， 矩阵二分求解矩阵 B 即可；

HDU 3509 Buge’s Fibonacci Number Problem

• 链接：HDU 3509 Buge’s Fibonacci Number Problem
• 题意：$$F(1)=f_1,F(2)=f_2,F(n)=aF(n-1)+bF(n-2)$$
• 给定 $f_1、f_2、a、b、n，m，k(k<=50)$，求 ${\sum }_{i=1}^{n}F(i{)}^k\mathrm{m}\mathrm{o}\mathrm{d}m$
• 难度：★★★★☆
• 题解：递推公式k次幂的
• 令 $G(n,k)=F(n{)}^k$，则根据递推公式，对 $G(n,k)$ 进行二项式展开，得到如下公式：
  $$\begin{array}{rl}G(n,k)& =F(n{)}^k\\ & =[aF(n-1)+bF(n-2){]}^k\\ & =C_k^0[aF(n-1){]}^k+C_k^1[aF(n-1){]}^{k-1}[bF(n-2)]+...+C_k^k[aF(n-1){]}^0[bF(n-2){]}^k\\ & =C_k^0{a}^{k}G(n-1,k)+C_k^1{a}^{k-1}bG(n-1,k-1)G(n-2,1)+...+C_k^k{a}^0{b}^{k}G(n-1,0)G(n-2,k)\end{array}$$
• 根据递推公式，列出所有涉及的列向量的值：

$$\begin{array}{r}\left[\begin{array}{c}G(n,k)\\ G(n,k-1)G(n-1,1)\\ ...\\ G(n,1)G(n-1,k-1)\\ G(n,0)G(n-1,k)\end{array}\right]\end{array}$$

• 再根据递推公式，将整个系数矩阵构造出来：

$$\begin{array}{rl}\left[\begin{array}{c}G(n,k)\\ G(n,k-1)G(n-1,1)\\ ...\\ G(n,1)G(n-1,k-1)\\ G(n,0)G(n-1,k)\end{array}\right]& =\left[\begin{array}{ccccc}{C}_k^0{a}^k& C_k^1{a}^{k-1}b& ...& C_k^{k-1}{a}^1{b}^{k-1}& C_k^k{a}^0{b}^k\\ C_{k-1}^0{a}^{k-1}& C_{k-1}^1{a}^{k-2}{b}^1& ...& C_{k-1}^{k-1}{a}^{k-1}{b}^0& 0\\ ...& ...& ...& ...& ...\\ C_1^0{a}^1{b}^0& C_1^1{a}^0{b}^1& ...& 0& 0\\ C_0^0{a}^0{b}^0& 0& ...& 0& 0\end{array}\right]\left[\begin{array}{c}G(n-1,k)\\ G(n-1,k-1)G(n-2,1)\\ ...\\ G(n-1,1)G(n-2,k-1)\\ G(n-1,0)G(n-2,k)\end{array}\right]\\ & =A^{n-2}\left[\begin{array}{c}G(2,k)\\ G(2,k-1)G(1,1)\\ ...\\ G(2,1)G(1,k-1)\\ G(2,0)G(1,k)\end{array}\right]\\ & =A^{n-2}\left[\begin{array}{c}{f}_2^k\\ f_2^{k-1}{f}_1\\ ...\\ f_2^1{f}_1^{k-1}\\ f_2^0{f}_1^k\end{array}\right]\end{array}$$

• 这里的矩阵 $A$ ，用 $A[i][j]$ 表示 第 $i$ 行 第 $j$ 列的元素，那么当 $i+j>k$ 的时候，元素为 0；否则，它的表示式如下：
  $$C_{k-i}^j{a}^{k-i-j}{b}^j$$
• 令 $S(n,k)={\sum }_{i=1}^{n}F(i{)}^k={\sum }_{i=1}^{n}G(i,k)$，套入原矩阵，得到求和矩阵如下：

$$\begin{array}{rl}\left[\begin{array}{c}S(n,k)\\ G(n,k)\\ G(n,k-1)G(n-1,1)\\ ...\\ G(n,1)G(n-1,k-1)\\ G(n,0)G(n-1,k)\end{array}\right]& ={\left[\begin{array}{cc}1& B\\ 0& A\end{array}\right]}^{n-2}\left[\begin{array}{c}{f}_1^k+f_2^k\\ f_2^k\\ f_2^{k-1}{f}_1\\ ...\\ f_2^1{f}_1^{k-1}\\ f_2^0{f}_1^k\end{array}\right]\end{array}$$

• 其中 B 为 一个行向量，并且只有一行，且为 矩阵 A 的第一行；二分求解即可；

HDU 3483 A Very Simple Problem

• 链接：HDU 3483 A Very Simple Problem
• 题意：给定如下求和式：
• 
• 其中 $1\le N,M\le 2\times 10^9,1\le x\le 50$，求最后的值；
• 难度：★★★★☆
• 题解：令 $f(N,x,a)=N^x{a}^N$，将 $N$ 拆分成 $N-1$ 和 $1$，然后对 $N^x$ 进行二项展开，得到：
  $$\begin{array}{rl}f(N,x,a)& =N^x{a}^N\\ & =[(N-1)+1{]}^x{a}^N\\ & =[C_x^0(N-1{)}^x+C_x^1(N-1{)}^{x-1}+...+C_x^{x-1}(N-1{)}^1+C_x^x(N-1{)}^0]a^N\\ & =a^N{C}_x^0(N-1{)}^x+a^N{C}_x^1(N-1{)}^{x-1}+...+a^N{C}_x^{x-1}(N-1{)}^1+a^N{C}_x^x(N-1{)}^0\\ & =a{C}_x^0(N-1{)}^x{a}^{N-1}+a{C}_x^1(N-1{)}^{x-1}{a}^{N-1}+...+a{C}_x^{x-1}(N-1{)}^1{a}^{N-1}+a{C}_x^x(N-1{)}^0{a}^{N-1}\\ & =a{C}_x^{0}f(N-1,x,a)+a{C}_x^{1}f(N-1,x-1,a)+...+a{C}_x^{x-1}f(N-1,1,a)+a{C}_x^{x}f(N-1,0,a)\end{array}$$
• 得到递推式后，可以构造出列向量和系数矩阵关系如下，其中 $C_n^m$ 为组合数：
  $$\begin{array}{rl}\left[\begin{array}{c}f(N,x,a)\\ f(N,x-1,a)\\ ...\\ f(N,1,a)\\ f(N,0,a)\end{array}\right]& =\left[\begin{array}{ccccc}a{C}_x^0& a{C}_x^1& ...& a{C}_x^{x-1}& a{C}_x^x\\ 0& a{C}_{x-1}^0& ...& a{C}_{x-1}^{x-2}& a{C}_{x-1}^{x-1}\\ ...& ...& ...& ...& ...\\ 0& 0& ...& a{C}_1^0& a{C}_1^1\\ 0& 0& ...& 0& a{C}_0^0\end{array}\right]\left[\begin{array}{c}f(N-1,x,a)\\ f(N-1,x-1,a)\\ ...\\ f(N-1,1,a)\\ f(N-1,0,a)\end{array}\right]\\ & =A^{N-1}\left[\begin{array}{c}f(1,x,a)\\ f(1,x-1,a)\\ ...\\ f(1,1,a)\\ f(1,0,a)\end{array}\right]\\ & =A^{N-1}\left[\begin{array}{c}a\\ a\\ ...\\ a\\ a\end{array}\right]\end{array}$$
• 现在是要求和，即 $s(N,x,a)={\sum }_{k=1}^{N}f(k,x,a)$，这里的 $a$ 为常数 $x$；
• 将列向量和矩阵各扩大一阶，得到如下式子：
  $$\begin{array}{rl}\left[\begin{array}{c}s(N,x,a)\\ f(N,x,a)\\ f(N,x-1,a)\\ ...\\ f(N,1,a)\\ f(N,0,a)\end{array}\right]& =\left[\begin{array}{cc}1& B\\ 0& A\end{array}\right]\left[\begin{array}{c}s(N-1,x,a)\\ f(N-1,x,a)\\ f(N-1,x-1,a)\\ ...\\ f(N-1,1,a)\\ f(N-1,0,a)\end{array}\right]\\ & ={\left[\begin{array}{cc}1& B\\ 0& A\end{array}\right]}^{N-1}\left[\begin{array}{c}a\\ a\\ a\\ ...\\ a\\ a\end{array}\right]\end{array}$$
• 其中矩阵 B 为 $1\times (x+1)$ 的行向量，且为矩阵 A 的第一行，构造出矩阵后进行二分求解即可；

HDU 4291 A Short problem

• 链接：HDU 4291 A Short problem
• 题意：已知 $g(0)=0，g(1)=1，g(n)=3g(n-1)+g(n-2)$，求 $g(g(g(n)))\mathrm{m}\mathrm{o}\mathrm{d}(10^9+7)$。
• 难度：★★★★☆
• 题解：循环节 + 矩阵二分
• 根据递推公式，构造矩阵如下：
  $$\left[\begin{array}{c}g(n)\\ g(n-1)\end{array}\right]={\left[\begin{array}{cc}3& 1\\ 1& 0\end{array}\right]}^{n-1}\left[\begin{array}{c}1\\ 0\end{array}\right]$$
• 将 $g(n)、g(g(n))$ 分别代入 $n$，得到递推矩阵如下：
  $$\left[\begin{array}{c}g(g(n))\\ g(g(n)-1)\end{array}\right]={\left[\begin{array}{cc}3& 1\\ 1& 0\end{array}\right]}^{g(n)-1}\left[\begin{array}{c}1\\ 0\end{array}\right]$$
  $$\left[\begin{array}{c}g(g(g(n)))\\ g(g(g(n))-1)\end{array}\right]={\left[\begin{array}{cc}3& 1\\ 1& 0\end{array}\right]}^{g(g(n))-1}\left[\begin{array}{c}1\\ 0\end{array}\right]$$
• 通过枚举矩阵次幂，得到
  $${\left[\begin{array}{cc}3& 1\\ 1& 0\end{array}\right]}^{222222224}\mathrm{m}\mathrm{o}\mathrm{d}(10^9+7)=\left[\begin{array}{cc}1& 0\\ 0& 1\end{array}\right]$$
• 于是，我们发现，原等式转化成求：
  $$\left[\begin{array}{c}g(g(g(n)))\\ g(g(g(n))-1)\end{array}\right]={\left[\begin{array}{cc}3& 1\\ 1& 0\end{array}\right]}^{g(g(n))\mathrm{m}\mathrm{o}\mathrm{d}222222224-1}\left[\begin{array}{c}1\\ 0\end{array}\right]$$
• 那么，问题转变为求 $g(g(n))\mathrm{m}\mathrm{o}\mathrm{d}222222224$，比原题降了一维；
• 通过同样方法进行枚举，得到：
  $${\left[\begin{array}{cc}3& 1\\ 1& 0\end{array}\right]}^{183120}\mathrm{m}\mathrm{o}\mathrm{d}(222222224)=\left[\begin{array}{cc}1& 0\\ 0& 1\end{array}\right]$$
• 于是有：
  $$\left[\begin{array}{c}g(g(n))\\ g(g(n)-1)\end{array}\right]={\left[\begin{array}{cc}3& 1\\ 1& 0\end{array}\right]}^{g(n)\mathrm{m}\mathrm{o}\mathrm{d}183120-1}\left[\begin{array}{c}1\\ 0\end{array}\right]$$
• 那么，问题转变为求 $g(n)\mathrm{m}\mathrm{o}\mathrm{d}183120$，回到了我们熟悉的矩阵二分；
• 这里需要注意的是，因为求矩阵次幂的时候指数上有 $-1$，所以次幂上的数取模后如果为0，需要变成模数本身，以防出现求矩阵的 $-1$ 次幂的情况；

HDU 5451 Best Solver

• 链接：HDU 5451 Best Solver
• 题意：$$y={(5+2\sqrt{6})}^{2^x+1}$$
  给定 $0<=x<2^{32}$，以及一个素数 $M(M\le 46337)$， 求 $[y]\%M$ 的值（ $[y]$ 表示 $y$ 的整数部分）；
• 难度：★★★☆☆
• 题解：$y$ 无论当 $x$ 取什么值，都一定可以表示成 $A_n+B_n\sqrt{6}$ 的形式，所以我们可以得到 $A_n和B_n$ 的递推关系如下：
  $$\left[\begin{array}{c}{A}_n\\ B_n\end{array}\right]=\left[\begin{array}{cc}5& 12\\ 2& 5\end{array}\right]\left[\begin{array}{c}{A}_{n-1}\\ B_{n-1}\end{array}\right]={\left[\begin{array}{cc}5& 12\\ 2& 5\end{array}\right]}^{n-1}\left[\begin{array}{c}5\\ 2\end{array}\right]$$
  又由于 $(5+2\sqrt{6})(5-2\sqrt{6})=1$，得到 $(5+2\sqrt{6}{)}^n(5-2\sqrt{6}{)}^n=1$，则 $(A_n+B_n\sqrt{6})(A_n-B_n\sqrt{6})=1$，得到：
  $$6B_n^2=A_n^2-1$$，从而推导出$$B_n\sqrt{6}=\sqrt{A_n^2-1}$$，于是 $B_n\sqrt{6}$ 的整数部分就是 $A_n-1$，所以给定 $x$ 的值，我们要求的就是 $(2A_{2^x+1}-1)\%M$，只要先求 $A_{2^x+1}\%M$ 的值即可；
  $$\left[\begin{array}{c}{A}_{2^x+1}\\ B_{2^x+1}\end{array}\right]={\left[\begin{array}{cc}5& 12\\ 2& 5\end{array}\right]}^{2^x}\left[\begin{array}{c}5\\ 2\end{array}\right]$$
  令矩阵 $$g(n)={\left[\begin{array}{cc}5& 12\\ 2& 5\end{array}\right]}^{2^n}$$
• 当 M = 2 时，答案永远是1，直接输出不用纠结；
• 当 M > 2 时，因为是素数，如果能够找到一个矩阵的K次幂，满足 ${\left[\begin{array}{cc}5& 12\\ 2& 5\end{array}\right]}^K=\left[\begin{array}{cc}1& 0\\ 0& 1\end{array}\right]$，那么必定可以满足：
• $$g(n)={\left[\begin{array}{cc}5& 12\\ 2& 5\end{array}\right]}^{2^n}={\left[\begin{array}{cc}5& 12\\ 2& 5\end{array}\right]}^{2^n\mathrm{m}\mathrm{o}\mathrm{d}K}$$
• 所以求解顺序就是： 利用循环节求 $K$、利用矩阵二分求 $g(n)$、最后计算 $(2A_{2^x+1}+M-1)\%M$；

HDU 1977 Odd Loving Bakers

• 链接：HDU 1977 Odd Loving Bakers
• 题意：有 $n(n<=100)$ 个人，每个人都有一个标记值，第 $i$ 个人初始标记值为 $X_i$ ，给出一张关系网，$a->b$ 有一条边代表 $a$ 能够给 $b$ 提供一点 标记值（但是需要满足提供人的标记值为奇数时才能提供），问经过 $t(t<=10^9)$次传递后，标记值为奇数的人的个数；
• 难度：★★★★☆
• 题解：令第 $i$ 个人在第 $t$ 次传递后的 标记值 数量为 $f(t,i)$，则它的标记值来源为对他有边并且上次标记值是奇数的玩家数量加上自己第 $t-1$ 次的标记值，有状态转移方程如下：
  $$f(t,i)=f(t-1,i)+\sum _{j=0}^{n-1}e(j,i)$$
  $$e(j,i)=\{\begin{array}{ll}1& A[j][i]=1,f(t-1,j)\mathrm{m}\mathrm{o}\mathrm{d}2=1\\ 0& \end{array}$$
• 第 $t-1$ 次 和 第 $t$ 次 矩阵递推关系如下：
  $$\begin{array}{rl}\left[\begin{array}{c}f(t,0)\\ f(t,1)\\ ...\\ f(t,n-2)\\ f(t,n-1)\end{array}\right]& =E\left[\begin{array}{c}f(t-1,0)\\ f(t-1,1)\\ ...\\ f(t-1,n-2)\\ f(t-1,n-1)\end{array}\right]+A\left[\begin{array}{c}f(t-1,0)\%2\\ f(t-1,1)\%2\\ ...\\ f(t-1,n-2)\%2\\ f(t-1,n-1)\%2\end{array}\right]\end{array}$$
• 其中 $E$ 代表单位矩阵, $A[i][j]$ 代表 $j$ 到 $i$ 有一条边；
• 由于最后需要求的就是 $f(t,i)$ 的奇偶性，所以如果 $f(t-1,i)$ 是偶数，对 $f(t,i)$ 相当于没有贡献，则将矩阵等式化简为：
  $$\begin{array}{rl}\left[\begin{array}{c}f(t,0)\%2\\ f(t,1)\%2\\ ...\\ f(t,n-2)\%2\\ f(t,n-1)\%2\end{array}\right]& =(E+A)\left[\begin{array}{c}f(t-1,0)\%2\\ f(t-1,1)\%2\\ ...\\ f(t-1,n-2)\%2\\ f(t-1,n-1)\%2\end{array}\right]\end{array}$$
• 化简后得到如下矩阵递推式，利用矩阵二分进行求解即可；
  $$\begin{array}{rl}\left[\begin{array}{c}f(t,0)\\ f(t,1)\\ ...\\ f(t,n-2)\\ f(t,n-1)\end{array}\right]\%2& =(E+A)\left[\begin{array}{c}f(t-1,0)\\ f(t-1,1)\\ ...\\ f(t-1,n-2)\\ f(t-1,n-1)\end{array}\right]\%2\\ & =(E+A{)}^{t-1}\left[\begin{array}{c}f(0,0)\\ f(0,1)\\ ...\\ f(0,n-2)\\ f(0,n-1)\end{array}\right]\%2\end{array}$$

HDU 5986 Fibonacci

HDU 5986 Fibonacci

• 题意：$$\begin{gathered} f(0)=0,f(1)=1 \\ f(n)=f(n-1)+f(n-2)(n\ge 2) \\ x(n)=f(x(n-1)) \end{gathered}$$ 给定 $x(0)和p$, 求最小的 n，满足 $x(n)\equiv x(0)(\mathrm{m}\mathrm{o}\mathrm{d}p)$ $(0<=x(0)<=10^9,p<=200000)$
• 难度：★★★★★
• 题解：循环节预处理 + 矩阵二分；
• 关于斐波那契数取余，观察如下表，对于任何一个模数，必定能够找到一个循环节，使得后面的数字产生循环，假设有 $f(i)\mathrm{m}\mathrm{o}\mathrm{d}p$ 的循环节为 $K_p$，则必然有 ：
  $$f(i)\mathrm{m}\mathrm{o}\mathrm{d}p=f(i\mathrm{m}\mathrm{o}\mathrm{d}{K}_p)\mathrm{m}\mathrm{o}\mathrm{d}p$$

i	0	1	2	3	4	5	6	7	8	9	10	11	12	13	14	15	16	17	18
f(i)	0	1	1	2	3	5	8	13	21	34	55	89	144	233	377	610	987	1597	2584
f(i)%1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
f(i)%2	0	1	1	0	1	1	0	1	1	0	1	1	0	1	1	0	1	1	0
f(i)%3	0	1	1	2	0	2	2	1	0	1	1	2	0	2	2	1	0	1	1
f(i)%4	0	1	1	2	3	1	0	1	1	2	3	1	0	1	1	2	3	1	0

• 如上图所示，$f(i)\mathrm{m}\mathrm{o}\mathrm{d}4$ 的循环节为6，则 $$f(15)\mathrm{m}\mathrm{o}\mathrm{d}4=f(15\mathrm{m}\mathrm{o}\mathrm{d}6)\mathrm{m}\mathrm{o}\mathrm{d}4=f(3)\mathrm{m}\mathrm{o}\mathrm{d}4=2$$

• 那么，来看如果在 $n$ 已知的情况下，对于$x(n)\mathrm{m}\mathrm{o}\mathrm{d}p$，满足：$$x(n)\mathrm{m}\mathrm{o}\mathrm{d}p=f(x(n-1))\mathrm{m}\mathrm{o}\mathrm{d}p$$，那么这里只要知道 $f(x)$ 相对于模 $p$ 的循环节为 $K_p$，则可以得出：$$x(n)\mathrm{m}\mathrm{o}\mathrm{d}p=f(x(n-1)\mathrm{m}\mathrm{o}\mathrm{d}{K}_p)\mathrm{m}\mathrm{o}\mathrm{d}p$$，而因为 $p$ 不大，所以这里的循环节可以通过枚举得出；

• 于是问题转化成求 $x(n-1)\mathrm{m}\mathrm{o}\mathrm{d}{K}_p$ 的问题，同样继续递归求解，得到：$$x(n-1)\mathrm{m}\mathrm{o}\mathrm{d}{K}_p=f(x(n-2)\mathrm{m}\mathrm{o}\mathrm{d}{K}_{K_p})\mathrm{m}\mathrm{o}\mathrm{d}{K}_p$$，通过这种方式不断递归，直到最后求到 $x(0)$，再回溯回来得到 $x(n)\mathrm{m}\mathrm{o}\mathrm{d}p$ 的值；

• 那么可以通过枚举 n，计算 $x(n)\mathrm{m}\mathrm{o}\mathrm{d}p$ ，如果等于 $x(0)\mathrm{m}\mathrm{o}\mathrm{d}p$ 则输出最小的 n；

• 然后再来考虑最后一个问题，什么情况下无解？

• 为了把问题具象化，这里用几个数组来代表中间的計算过程，$P_i$ 代表循环节链，$F_i$ 代表斐波那契数模上循环节的值，则有如下公式： $$F_i=f(F_{i+1})\mathrm{m}\mathrm{o}\mathrm{d}{P}_i$$
  

• 如上图所示，其中 $n$ 是我们需要求的，是个未知数，$P_0$ 已知，$P_i$ 可以通过枚举预处理出来，并且一定能在有限步里达到循环，满足如下规律：
  $$P_i=\{\begin{array}{ll}p& i=0\\ GetCycle(P_{i-1})& 0<=i<=m\\ P_{m-1}& i>m\end{array}$$

• 其中 GetCycle 函数实现如下，利用记忆化避免重复计算：

int GetCycle(int p) {
	if (p == 1) {
		Cycle[p] = 1;
	}
	if (Cycle[p]) {
		return Cycle[p];
	}
	int fprepre = 0;
	int fpre = 1;
	int fnow;
	for (int i = 1;; ++i) {
		fnow = (fpre + fprepre) % p;
		if (fnow == 1 && fpre == 0) {
			Cycle[p] = i;
			return i;
		}
		fprepre = fpre;
		fpre = fnow;
	}
}

• 由于 $m$ 一般比较小，所以 $F_0到F_{m-1}$ 计算可以采用枚举法；$F_m到F_n$ 计算过程中，模数都是 $P_m$，所以必然会产生循环，所以当产生循环还是没有得到结果的时候，就是无解了；

HDU 6755 Fibonacci Sum

• 链接：HDU 6755 Fibonacci Sum
• 题意：给定$N,C,K(N,C<=10^{18},K<=10^5)$，对于斐波那契数列 $f$，求
  $$S(N)=\sum _{i=0}^{N}f(iC{)}^K\mathrm{m}\mathrm{o}\mathrm{d}1000000009$$
• 难度：★★★★★
• 题解：二次剩余 + 费马小定理 + 斐波那契数列通项公式 + 二分快速幂
• 斐波那契数列的通项公式如下：
  $$f(n)=\frac{1}{\sqrt{5}}[(\frac{1+\sqrt{5}}{2}{)}^n-(\frac{1-\sqrt{5}}{2}{)}^n]$$
• 令 $x=\frac{1+\sqrt{5}}{2}$, $y=\frac{1-\sqrt{5}}{2}$, $z=\frac{1}{\sqrt{5}}$，将原公式化简为：
  $$f(n)=z[x^n-y^n]$$，通过代入 $iC$ 和 $K$，得到：
  f ( i C ) K = { z [ x i C − y i C ] } K f(iC)^K = \{ z [ x^{iC} - y^{iC}] \} ^ K f(iC)K={z[xiC−yiC]}K
• 利用二项式定理进行展开，得到：
  f ( i C ) K = { z [ x i C − y i C ] } K = z K [ x i C − y i C ] K = z K [ C K 0 x i C K ( − 1 ) 0 + C K 1 x i C ( K − 1 ) ( − 1 ) 1 y i C + . . . + C K K ( − 1 ) K y i C K ] = z K ∑ j = 0 K C K j x i C ( K − j ) ( − 1 ) j y i C j \begin{aligned} f(iC)^K &= \{ z [ x^{iC} - y^{iC}] \} ^ K \\ &= z^K [ x^{iC} - y^{iC}]^K \\ &= z^K [ C_K^{0}x^{iCK}(-1)^0+ C_K^{1}x^{iC(K-1)}(-1)^1y^{iC} + ... + C_K^{K}(-1)^Ky^{iCK}] \\ &= z^K \sum_{j=0}^{K} C_K^{j}x^{iC(K-j)}(-1)^jy^{iCj} \end{aligned} f(iC)K​={z[xiC−yiC]}K=zK[xiC−yiC]K=zK[CK0​xiCK(−1)0+CK1​xiC(K−1)(−1)1yiC+...+CKK​(−1)KyiCK]=zKj=0∑K​CKj​xiC(K−j)(−1)jyiCj​
• 从而得到求和公式 $S(N)$ 满足：
  $$\begin{array}{rl}S(N)& =\sum _{i=0}^{N}f(iC{)}^K\mathrm{m}\mathrm{o}\mathrm{d}1000000009\\ & =z^K\sum _{i=0}^N[\sum _{j=0}^K{C}_K^j{x}^{iC(K-j)}(-1{)}^j{y}^{iCj}]\mathrm{m}\mathrm{o}\mathrm{d}1000000009\end{array}$$
• 这里出现了二维求和，基于 $N(N<=10^{18})$ 和 $K(N<=10^5)$ 的取值范围，可以将两个求和循环进行一个交换，枚举 $j(0<=j<=K)$，计算 ${\sum }_{i=0}^N{C}_K^j{x}^{iC(K-j)}(-1{)}^j{y}^{iCj}]$；
  $$z^K\sum _{i=0}^N[\sum _{j=0}^K{C}_K^j{x}^{iC(K-j)}(-1{)}^j{y}^{iCj}]=z^K\sum _{j=0}^K[\sum _{i=0}^N{C}_K^j{x}^{iC(K-j)}(-1{)}^j{y}^{iCj}]$$
• 枚举 $j$ 之后，$j$ 就成了常数，那么我们可以令 $G_j(i)=C_K^j{x}^{iC(K-j)}(-1{)}^j{y}^{iCj}$，来看下：
  $$\begin{array}{rl}\frac{G_j(i)}{G_j(i-1)}& =\frac{C_K^j{x}^{iC(K-j)}(-1{)}^j{y}^{iCj}}{C_K^j{x}^{(i-1)C(K-j)}(-1{)}^j{y}^{(i-1)Cj}}\\ & =x^{C(K-j)}{y}^{Cj}\\ & =q_j\end{array}$$
• 我们发现 $G_j(i)$ 是一个等比数列，那么通过等比数列通项公式，得到 $G_j(i)$ 的表达式如下：
  $$\begin{array}{rl}{G}_j(N)& =G_j(0)q_j^N\\ & =(-1{)}^j{C}_K^j{q}_j^N\end{array}$$
• 根据等比数列求和公式，得到：
  $$\begin{array}{rl}S(N)& =\sum _{i=0}^{N}f(iC{)}^K\mathrm{m}\mathrm{o}\mathrm{d}1000000009\\ & =z^K\sum _{j=0}^K\sum _{i=0}^N{G}_j(i)\mathrm{m}\mathrm{o}\mathrm{d}1000000009\\ & =z^K\sum _{j=0}^K[(-1{)}^j{C}_K^j\frac{q_j^{N+1}-1}{q_j-1}]\mathrm{m}\mathrm{o}\mathrm{d}1000000009\end{array}$$
• 对于二次剩余，我们可以通过枚举求出最小的 $X$（我的机子大概 5 秒硬算出来） ：
  $$\begin{gathered} X^2\equiv 5(\mathrm{m}\mathrm{o}\mathrm{d}1000000009) \\ X=383008016 \end{gathered}$$
• 那么就能够得到：
  $$383008016\equiv \sqrt{5}(\mathrm{m}\mathrm{o}\mathrm{d}1000000009)$$
• 由于 $x=\frac{1+\sqrt{5}}{2}$, $y=\frac{1-\sqrt{5}}{2}$, $z=\frac{1}{\sqrt{5}}$，所以现在需要求这三者模上 $1000000009$ 的值，利用线性同余分别求出 2 和 $\sqrt{5}$ 在模 $1000000009$ 域上的逆元，得到：
  $$\begin{gathered} 276601605\equiv {\sqrt{5}}^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}1000000009) \\ 500000005\equiv 2^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}1000000009) \end{gathered}$$
• 于是可以根据这个，计算所有 $x,y,z$ 在模 1000000009 域上的值，如下：
  $$\begin{gathered} x=\frac{1+\sqrt{5}}{2}=(1+383008016)\ast 500000005(\mathrm{m}\mathrm{o}\mathrm{d}1000000009)=691504013 \\ y=\frac{1-\sqrt{5}}{2}=(1-383008016)\ast 500000005(\mathrm{m}\mathrm{o}\mathrm{d}1000000009)=308495997 \\ z=\frac{1}{\sqrt{5}}={\sqrt{5}}^{-1}(\mathrm{m}\mathrm{o}\mathrm{d}1000000009)=276601605 \end{gathered}$$
• 由于 $q_j=x^{C(K-j)}{y}^{Cj}$，利用二分快速幂求解得出 $q_j$，再求 $q_j-1$的逆元，从而求得 $S(N)$ 的值。
• 这里需要求组合数 $C_K^j\mathrm{m}\mathrm{o}\mathrm{d}1000000009$，可以利用递推公式如下：
  $$C_K^j=\{\begin{array}{ll}1& j=0\\ C_K^{j-1}\frac{K-j+1}{j}& 0<j<=K\end{array}$$
• 递推公式 O(1)，需要求 $j$ 的逆元；
• 这题时间卡的比较紧，需要注意一些优化点， 如下：
• 1）求逆元采用费马小定理：$a^{-1}\equiv a^{p-2}\mathrm{m}\mathrm{o}\mathrm{d}p$；
• 2）预处理 $[1,10^5]$ 的 逆元 存储在数组中；
• 3）预处理 5 的二次剩余的 $i(i<=10^5)$ 次幂存储在数组中；
• 4）利用欧拉定理进行降幂：$a^b\equiv a^{b\mathrm{m}\mathrm{o}\mathrm{d}p-1}\mathrm{m}\mathrm{o}\mathrm{d}p$；
• 5）$q_j=x^{C(K-j)}{y}^{Cj}$ 的可以通过递推把复杂度从 $O(lo{g}_n)$ 将为 $O(1)$：$q_j=q_{j-1}\frac{y^C}{x^C}$；
• 6）$+1$和$-1$用加减法代替乘法；
• 7）用 G++ 提交；

HDU 4959 Poor Akagi

• 链接：HDU 4959 Poor Akagi
• 题意：给定递推数列 $$L(n)=\{\begin{array}{ll}2& n=0\\ 1& n=1\\ L(n-1)+L(n-2)& n>1\end{array}$$，给出 $n(n<=10^{18})，k<=10^5$，求 ${\sum }_{i=0}^{n}L(i{)}^{k}mod10^9+7$
• 难度：★★★★★
• 题解：费马小定理 + 卢卡斯数通项公式 + 矩阵二分快速幂
• 这类问题一般有以下几种做法：
• 1）矩阵：线性递推式 加上取模，很容易联想到的是矩阵二分幂，构造的矩阵一般是 $(k+1)\times (k+1)$ 这个量级的，但是这里 $k$ 很大，所以无法用矩阵来做。
• 2）循环节：再来看被模数是固定的，那么能否找到一个 $10^6$ 量级的循环节，然后存下循环节内所有的数，再利用循环节求解，很不幸，循环节是 $200000016$ ，太大了。
• 3）二次剩余 + 通项公式：另一种做法就是通项公式了。首先求出这个序列的通项公式，通项公式求解方法有很多，类似斐波那契数列的推导方法，这里不再累述，可以参考这篇文章：斐波那契数列通项公式推导
• 我们直接给出这个数列的通项公式如下：$$L(n)=(\frac{1+\sqrt{5}}{2}{)}^n+(\frac{1-\sqrt{5}}{2}{)}^n$$
• 那么这里存在 $\sqrt{5}$，所以要对模数 1000000007 找 5 的二次剩余，结果没有二次剩余，所以这种方法也行不通；
• 4）二次域运算：还是从通项公式出发，令 $p=\frac{1+\sqrt{5}}{2}，q=\frac{1-\sqrt{5}}{2}$，则有 $$L(n)=p^n+q^n$$
• 那么对 $L(n{)}^k$ 进行二项展开得到： $$\begin{array}{rl}L(n{)}^k& =(p^n+q^n{)}^k\\ & =C_k^0{p}^{nk}+C_k^1{p}^{n(k-1)}{q}^n+...+C_k^{k-1}{p}^n{q}^{n(k-1)}+C_k^k{q}^{nk}\\ & =\sum _{j=0}^k{C}_k^j{p}^{n(k-j)}{q}^{nj}\end{array}$$
• 那么 ${\sum }_{i=0}^{n}L(i{)}^k$ 就是将上述二项式再进行一次求和，得到：
  $$\sum _{i=0}^{n}L(i{)}^k=\sum _{i=0}^n\sum _{j=0}^k{C}_k^j{p}^{i(k-j)}{q}^{ij}$$
• 然后将 $i$ 和 $j$ 的位置进行一个交换，得到：
  $$\sum _{i=0}^{n}L(i{)}^k=\sum _{j=0}^k{C}_k^j\sum _{i=0}^n{p}^{i(k-j)}{q}^{ij}$$
• 这里的组合数，可以通过递推在均摊 $O(1)$ 的时间复杂度内求得：$$\begin{array}{rl}{C}_k^j& =\frac{k!}{j!(k-j)!}\\ & =\frac{k!(k-j+1)}{j(j-1)!(k-j+1)!}\\ & =\frac{k-j+1}{j}{C}_k^{j-1}\end{array}$$
• 由于 $k$ 的范围远小于 $n$，我们可以枚举 $j=0...k$，然后再想办法计算内层循环，这样 $j$ 在每次循环中就变成了常量，我们只需要关心变量的部分，即：$$s(i)=\sum _{i=0}^n{p}^{i(k-j)}{q}^{ij}$$
• 令 $a(i)=[p^{k-j}{q}^j{]}^i$，则有：$$\begin{array}{rl}a(i)& =[p^{k-j}{q}^j{]}^i=[(\frac{1+\sqrt{5}}{2}{)}^{k-j}(\frac{1-\sqrt{5}}{2}{)}^j{]}^i\\ & =[2^{-k}(1+\sqrt{5}{)}^{k-j}(1-\sqrt{5}{)}^j{]}^i\\ & =(X+Y\sqrt{5}{)}^i\end{array}$$
• 构造矩阵：$$\left[\begin{array}{c}a(i{)}_X\\ a(i{)}_Y\end{array}\right]={\left[\begin{array}{cc}X& 5Y\\ Y& X\end{array}\right]}^n\left[\begin{array}{c}1\\ 0\end{array}\right]$$
• 由于 $\sqrt{5}$ 那部分系数一定是会被消去的，所以我们可以不管，只需要计算有理数部分的域，将求和部分也代入矩阵中，得到：$$\left[\begin{array}{c}s(i{)}_X\\ a(i{)}_X\\ a(i{)}_Y\end{array}\right]={\left[\begin{array}{ccc}1& X& 5Y\\ 0& X& 5Y\\ 0& Y& X\end{array}\right]}^n\left[\begin{array}{c}1\\ 1\\ 0\end{array}\right]$$
• 利用矩阵二分求解 $s(i{)}_X$ 再乘上 $C_k^j$ 加和就是答案了。