本次学习的内容是 矩阵快速幂
分为矩阵乘法和快速幂两个部分;
介绍一下单位矩阵
左上至右下对角线为1,其余部分为0;
矩阵乘法基本步骤
题目大意
n*n矩阵,求其k次方结果;
#include<iostream>
#include<cstdio>
using namespace std;
struct matrix{
long long m[105][105];
}matri;
const int q=1e9+7;
long long n,k;
matrix ksm(matrix a,matrix b)//矩阵乘法
{
matrix c;
for(int i=1;i<=n;i++) //分别行,列推移
for(int j=1;j<=n;j++)
{
c.m[i][j]=0;
for(int p=1;p<=n;p++)
c.m[i][j]=(c.m[i][j]+a.m[i][p]*b.m[p][j])%q;
}
return c;
}
matrix ready(matrix a,long long k)//快速幂处理
{
k--;
matrix ko=a,vic=a;
for(;k;k>>=1,vic=ksm(vic,vic))
if(k&1)
ko=ksm(ko,vic);
return ko;
}
int main()
{
matrix A,C;
scanf("%lld%lld",&n,&k);
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
scanf("%lld",&A.m[i][j]);
C=ready(A,k);
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
printf("%lld ",C.m[i][j]);
printf("\n");
}
return 0;
}