#include<bits/stdc++.h>
using namespace std;
typedef long long int ll;
const int maxn=1e2+5;
const ll mod= 1e9+7;
struct matrix{
    ll m[maxn][maxn];
}ans,res;
ll n,N;

matrix mul(matrix a,matrix b){
    matrix tmp;
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            tmp.m[i][j]=0;
        }
    }
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            for(int k=1;k<=n;k++){
                tmp.m[i][j]=(tmp.m[i][j]+(a.m[i][k]*b.m[k][j])%mod)%mod;
            }
        }
    }
    return tmp;
}
void quickpow(ll N,ll n){
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            if(i==j)ans.m[i][j]=1;
            else ans.m[i][j]=0;
        }
    }
    while(N){
        if(N&1)res=mul(res,ans);
        ans=mul(ans,ans);
        N=N>>1;
    }
}
int main(){
    scanf("%lld%lld",&n,&N);
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            scanf("%lld",&res.m[i][j]);
        }
    }
    quickpow(N,n);
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            cout<<res.m[i][j]<<" ";
        }
        cout<<endl;
    }
    return 0;
}