在讲矩阵快速幂之前，要说清楚矩阵乘法是什么。

矩阵乘法，用来干什么，我这个小蒟蒻并不知道。

但是，有大佬对于他的作用可能有着自己的理解

【看大佬blog】

我这里找的题目是洛谷3390 矩阵快速幂的模板

方法和整数的快速幂是一样的

自己重载一下乘号即可

顺带推荐一个大佬对于矩阵快速幂的理解（%%%songyuchen大佬）
%%%大佬

已经有大佬对于矩阵乘法的方法的讲解了，我就不再多说
直接贴一下矩阵快速幂的代码

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
#define MAX 101
#define MOD (1000000007)
long long n,k;
struct yl//矩阵 
{
	   int n;//大小
	   long long g[MAX][MAX];  
};
yl operator*(yl a,yl b)//定义乘法 
{
	   int n=a.n;
	   yl cool;
	   memset(cool.g,0,sizeof(cool.g));
	   for(int i=1;i<=n;++i)
	   	    for(int j=1;j<=n;++j)
				 for(int k=1;k<=n;++k)
				        cool.g[i][j]=(cool.g[i][j]+1ll*a.g[i][k]*b.g[k][j]%MOD)%MOD;
	   cool.n=n;
	   return cool;
}
void write(yl a)
{
	   int n=a.n;
	   for(int i=1;i<=n;++i)
	   {
	       for(int j=1;j<=n;++j)
		      cout<<a.g[i][j]<<' ';
		   cout<<endl;
	   }
}
yl Pow(yl a,long long b)//a的b次方
{
	   if(b==1)return a;
	   yl s=Pow(a,b/2);
	   s.n=a.n;
	   s=s*s;
	   if(b&1)s=s*a;
	   return s;
}

int main()
{   
       cin>>n>>k;
       yl Six;
       Six.n=n;
       for(int i=1;i<=n;++i)
          for(int j=1;j<=n;++j)
             cin>>Six.g[i][j];
       write(Pow(Six,k));
	   return 0;
}