矩阵快速幂题目小结

1、HDU - 4686 (https://vjudge.net/problem/HDU-4686)

题意：

已知 $a_0=A0,a_i=a_{i-1}\ast AX+AY,b_0=B0,b_i=b_{i-1}\ast BX+BY$，给定 $n,A0,AX,AY,B0,BX,BY$，求 $\sum _{i=0}^{n-1}{a}_i{b}_i$。$(n\le 10^{18})$

解题思路：

令 $c_i=a_i{b}_i,S_{n-1}=\sum _{i=0}^{n-1}{c}_i$，易得递推关系
$$\begin{gathered} c_i=AX\ast BX\ast c_{i-1}+AX\ast BY\ast a_{i-1}+AY\ast BX\ast b_{i-1}+AY\ast BY \\ {S}_{i-1}=S_{i-2}+c_{i-1} \\ 设M_{i-1}=A{M}_{i-2},M_{-1}=(0,A0\ast B0,A0,B0,1{)}^T \\ {M}_{i-1}=(S_{i-1},c_i,a_i,b_i,1{)}^T \\ {M}_{i-2}=(S_{i-2},c_{i-1},a_{i-1},b_{i-1},1{)}^T \end{gathered}$$
不难写出转移矩阵 $A$，则 $M_{n-1}=A^n{M}_{-1}$矩阵快速幂即可。
 
2、HDU - 5667 (https://vjudge.net/problem/HDU-5667)

题意：

$$f_n=\{\begin{array}{ll}1& n=1\\ a^b& n=2\\ a^b{f}_{n-1}^c{f}_{n-2}& otherwise\end{array}$$
给定 $n,a,b,c,p$，求 $f_{n}modp$。$(n\le 10^{18})$

解题思路：

取对数，
$$\begin{gathered} g_n=lo{g}_a{f}_n=\{\begin{array}{ll}0& n=1\\ b& n=2\\ b+c{g}_{n-1}+g_{n-2}& otherwise\end{array} \\ {f}_n=a^{g_n} \end{gathered}$$
矩阵快速幂求出 $g_n$，再快速幂求出 $f_n$。
 
3、HDU - 5451 (https://vjudge.net/problem/HDU-5451)

题意：

给定 $x,M$，求 $\lfloor (5+2\sqrt{6}{)}^{1+2^x}\rfloor modM$。$(x<2^{32},M\le 46337且是素数)$

解题思路：

先考虑 $f(n)=\lfloor (5+2\sqrt{6}{)}^n\rfloor$，这像是二阶线性递推式的通项的一部分，考虑下式
$$\begin{gathered} g(n)=ag(n-1)+bg(n-2) \\ 特征方程r^2-ar-b=0\Delta =\sqrt{a^2+4b} \\ 特征根r=\frac{a\pm \sqrt{a^2+4b}}{2} \\ 若r_1\ne r_2，g(n)=c_1{r}_1^n+c_2{r}_2^n \\ 构造g(n)=(5+2\sqrt{6}{)}^n+(5-2\sqrt{6}{)}^n \\ 得a=10,b=-1,即g(n)=10g(n-1)-g(n-2) \\ 而0<(5-2\sqrt{6}{)}^n\le 1,故f(n)=g(n)-1 \end{gathered}$$
考虑矩阵快速幂求解 $g(1+2^x)$，由于 $1+2^x$ 太大，需要进一步寻找循环节优化，由于 $M$ 不大，暴力求循环节可行，通用方法见下题。
 
4、HDU - 5674 (https://vjudge.net/problem/HDU-5674)

题意：

给定 $a,b,x,y$，求 $f(x^y)$，其中 $f(n)=(a+\sqrt{b}{)}^n+(a-\sqrt{b}{)}^{n}mod10^9+7$。$(x\le 50,y\le 10^{18})$

解题思路：

由上一题特征根法得到通项
$$g(n)=cg(n-1)+dg(n-2),c=2a,d=b-a^2$$
依旧需要求循环节来优化，记 $per$ 为一个周期，对于上述二阶常系数线性递推式$(cmodp\ne 0,dmodp\ne 0)$，这里不加证明给出：
$$\begin{gathered} per=\{\begin{array}{llc}p-1& (\frac{\Delta }{p})=1& \\ (p+1)(p-1)& (\frac{\Delta }{p})=-1& (p为奇素数)\\ p(p-1)& (\frac{\Delta }{p})=0& \end{array} \\ 其中,(\frac{\Delta }{p})为勒让德符号,若p=2,有per=3 \\ 统一化有per=p^2-1(p为素数) \\ 若cmodp=0或dmodp=0,数列很容易得到，不再讨论 \end{gathered}$$
 
5、CodeForces - 719E (https://vjudge.net/problem/CodeForces-719E)

题意：

给定一个长度为 $n$ 的数组 $a$，再有 $m$ 次操作，每次操作为：① $1,l,r,x$，将 $a_l,a_{l+1},...,a_r$ 加上 $x$，② $2,l,r$，询问 $\sum _{i=l}^{r}f(a_i)$，其中 $f(n)$ 为斐波那契数列。$(n,m\le 10^5)$

解题思路：

考虑用矩阵来得到斐波那契数列的值，令 $M_i=A{M}_{i-1}$，

其中 $M_i=\left[\begin{array}{c}f(i+1)\\ f(i)\end{array}\right],A=\left[\begin{array}{cc}1& 1\\ 1& 0\end{array}\right],M_0=\left[\begin{array}{c}1\\ 0\end{array}\right],M_n=A^n{M}_0$

先看求和， $\sum _{i=l}^r{M}_{a_i}=\sum _{i=l}^r{A}^{a_i}{M}_0=(\sum _{i=l}^r{A}^{a_i})M_0$

再看更新，$\sum _{i=l}^r{M}_{a_i+x}=(\sum _{i=l}^r{A}^{a_i+x})M_0=(\sum _{i=l}^r{A}^{a_i})A^x{M}_0$

线段树维护 $\sum _{i=l}^r{A}^{a_i}$，更新操作为区间乘矩阵 $A^x$。