首先是矩阵乘法
在矩阵乘法中满足以下运算律:
1.(AB)C=a(BC) 2.(A+B)C=AC+BC 3.C(A+B)=CA+CB
在普通的乘法中,一个数乘1还是等于它本身,在矩阵乘法中也有这么一个“1”,它就是单位矩阵
不同于普通乘法中的单位1,对于不同矩阵他们的单位矩阵大小是不同的
对于n∗mn*mn∗m的矩阵,它的单位矩阵大小为m∗mm*mm∗m,对于m∗nm*nm∗n的矩阵,它的单位矩阵大小为n∗nn*nn∗n
也就是说单位矩阵都是正方形的,这是因为只有正方形的矩阵能保证结果和前一个矩阵形状相同
单位矩阵的元素非0即1,从左上角到右下角的对角线上元素皆为1,其他皆为0
矩阵快速幂就是把快速幂和矩阵乘法结合在一起
#include<bits/stdc++.h>
using namespace std;
inline int read(){
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch == '-') f=-1 ; ch=getchar();}
while(ch>='0'&&ch&l