题目链接
存一下矩阵快速幂的板子。= =
AC代码:
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=105;
const int mo=1e9+7;
ll n,m;
struct ahaha
{
ll a[maxn][maxn]; //一定要用long long存矩阵,否则在过程中会爆掉
ahaha()
{
memset(a,0,sizeof a);
}
inline void build()
{ //建造单位矩阵
for(int i=1;i<=n;++i)a[i][i]=1;
}
}a;
ahaha operator *(const ahaha &x,const ahaha &y)
{ //重载运算符
ahaha z;
for(int k=1;k<=n;++k)
for(int i=1;i<=n;++i)
for(int j=1;j<=n;++j)
z.a[i][j]=(z.a[i][j]+x.a[i][k]*y.a[k][j]%mo)%mo;
return z;
}
ll k;
ahaha qk(ahaha a,ll b)
{
ahaha ans;
ans.build();
while(b)
{
if(b&1) ans=ans*a;
a=a*a;
b=b>>1;
}
return ans;
}
int main()
{
scanf("%lld%lld",&n,&m);
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
scanf("%lld",&a.a[i][j]);
}
}
ahaha fin=qk(a,m);
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
if(j!=n) printf("%lld ",fin.a[i][j]);
else printf("%lld\n",fin.a[i][j]);
}
}
return 0;
}