矩阵快速幂

【模板】矩阵快速幂

去年这个时候也是在学矩阵快速幂，但去年抄的$AC$代码自己都看不懂了。。。
神仙题解
主要是用来优化递推的，对决策递推式构建一个矩阵然后快速幂就好了
不知道为什么写的函数乘就错了，写重载运算符就对了，初始$build$函数最好写在结构体里，或者记得加上$\&$
记得$longlong$

#include <bits/stdc++.h>
#define ll long long
using namespace std;
const int N=150;
const int Mod=1e9+7;
ll n;
ll k;
inline ll read(){
	ll cnt=0,f=1;char c=getchar();
	while(!isdigit(c)){if(c=='-')f=-f;c=getchar();}
	while(isdigit(c)){cnt=(cnt<<1)+(cnt<<3)+(c^48);c=getchar();}
	return cnt*f;
}
struct matri{
	ll d[N][N];
	matri() {memset(d,0,sizeof(d));}
	inline void build(){for(int i=1;i<=n;++i) d[i][i]=1;}
}a;
inline void init(){
	n=read();k=read();
	for(int i=1;i<=n;++i)
		for(register int j=1;j<=n;++j)
			a.d[i][j]=read();
}
matri operator *(const matri &a,const matri &b){
	matri z;
	for(int k=1;k<=n;++k)
		for(int i=1;i<=n;++i)
			for(register int j=1;j<=n;++j)
				z.d[i][j]=(z.d[i][j]+a.d[i][k]*b.d[k][j]%Mod)%Mod;
	return z;
}
signed main(){
	init();matri ans;ans.build();
	for(;k;k>>=1){
		if(k&1) ans=ans*a;
		a=a*a;
	}
	for(int i=1;i<=n;++i){
		for(register int j=1;j<=n;++j) cout<<ans.d[i][j]<<' ';
		cout<<endl;
	}
	return 0;
}