基础数论算法（八） 矩阵乘法与线性齐次递推公式的快速求值

矩阵乘法可是zhx钦定说考的可能性挺大的东西！

---

###矩阵乘法
一个n*m的矩阵乘一个m*p的矩阵会得到一个n*p的矩阵。
因此有一道区间dp模板叫矩阵链乘法，不过这个不是重点。
矩阵乘法如果要简单的说，就是行乘每一列放在对应的行上。这么说可能不太清楚，不过我觉得代码比一堆$\sum$清晰的多。

void multi(int m,int n,int p){
    for(int i=0;i<m;++i)
        for(int j=0;j<p;++j)
            for(int k=0;k<n;++k)
                c[i][j]+=a[i][k]*b[k][j];            
}

---

###特征根方程
为什么要提特征根方程？因为觉得很有趣如果解出来是整数说不定可以少些五行矩阵乘法。
对于f(n)=af(n-1)+bf(n-2)这种递推公式，我们可以用特征根方程求出通项。顺便一提，这种方法写大题的时候不能用，老老实实待定系数吧。
首先，把原通项写成f(n)-af(n-1)-bf(n-2)=0的形式。
接下来，解一个方程：$x^2-ax-b=0$，假设解为$x_1,x_2$
那么我们称这两个解是这个数列的特征根。
这里要注意，如果无实数解，那么这个数列一般是具有周期性的。
如果$x_1\ne x_2$，那么$a_n=c_1{x}_1^n+c_2{x}_2^n$，其中$c_1,c_2$带进去两个数就知道了。
例如：f(n)=2f(n-1)+3f(n-2),f(1)=5,f(2)=7
x 2 − 2 x − 3 = 0 ⇒ x 1 = 3 , x 2 = − 1 x^2-2x-3=0 \Rightarrow x_1=3,x_2=-1 x