題目來源：https://www.51nod.com/onlineJudge/questionCode.html#!problemId=1137
【題意】
中文題意不在敘述，只是讓求一個矩陣的乘法而已。
【思路】
終於可以做矩陣快速冪的題了，這個專題一直被我拖到現在，作為一個弱弱內心無比難受，終於可以把它學了，十分開心。
此題非常裸，直接過載一個乘號運算子就可以了，只不過矩陣乘法要知道怎麼回事。不知道的可以看下這位大佬的部落格，知道的也就可以看下，加深理解：矩陣乘法
並且初始矩陣這個概念也要知道，這些可以看下大佬的矩陣快速冪視訊
【程式碼】

#include<set>
#include<map>
 

#include<stack>
#include<cmath>
#include<queue>
#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<limits.h>
#include<algorithm>
#define mem(a,b) memset(a,b,sizeof(a))
using namespace std;
const int mod=1e9+7;
typedef
 
 unsigned long long ll;
typedef long long LL;
const int N=100+10;
LL n,m;
struct mat
{
    LL a[N][N];
};
mat operator*(mat s,mat t)
{
    mat r;
    LL x;
    mem(r.a,0);
    for(int i=1; i<=n; i++)
        for(int j=1; j<=n; j++)
        {
            x=0;
            for(int k=1; k<=n; k++)
                x+=(s.a[i][k]*t.a[k][j])%mod;
            r.a[i][j]=x%mod;
        }
    return
 
 r;
}
void init(mat& s)
{
    for(int i=1; i<=n; i++)
        s.a[i][i]=1;
}
mat pow_quick(mat r)
{
    mat ans;
    init(ans);
    while(m)
    {
        if(m&1)
            ans=ans*r;
        r=r*r;
        m>>=1;
    }
    return ans;
}
int main()
{

    mat t;
    scanf("%lld%lld",&n,&m);
    for(int i=1; i<=n; i++)
        for(int j=1; j<=n; j++)
            scanf("%lld",&t.a[i][j]);
    t=pow_quick(t);
    for(int i=1; i<=n; i++)
    {
        printf("%lld",t.a[i][1]);
        for(int j=2; j<=n; j++)
            printf(" %lld",t.a[i][j]);
        printf("\n");
    }
}