// Slower C++ solution:
#include<bits/stdc++.h>

using namespace std;

long long n,k,mod=1e9+7;

struct matrix{

	long long x[100][100];	

}mat;

matrix matrixidentity(long long sz)

{

	matrix res;

	for(int i=0;i<=sz;i++)

	{

		for(int j=0;j<=sz;j++)

		{

			res.x[i][j] = 0;

			if(i==j) res.x[i][j] = 1;

		}

	}

	return res;

}

matrix multiplication(matrix a,matrix b,int sz,long long md)

{

	matrix res;

	for(int i=0;i<sz;i++)

	{

		for(int j=0;j<sz;j++)

		{

			res.x[i][j] = 0;

			for(int k=0;k<sz;k++)

			{

				res.x[i][j]+=(a.x[i][k]%md*b.x[k][j]%md)%md;

				res.x[i][j]%=md;

			}

		}

	}

	return res;

}

matrix matpow(matrix a,long long p,long long md)

{

	matrix res;

	res = matrixidentity(n);

	while(p)

	{

		if(p%2) res = multiplication(res,a,n,md);

		

		a = multiplication(a,a,n,md);

		p/=2;

	}

	

	return res;

}

int main()

{

	scanf("%lld %lld",&n,&k);

	for(int i=0;i<n;i++)

	{

		for(int j=0;j<n;j++)

		{

			scanf("%lld",&mat.x[i][j]);

		}

	}

	mat = matpow(mat,k,mod);

	//mat = multiplication(mat,mat,n,mod);

	long long ans = 0;

//	printf("\n");

	for(int i=0;i<n;i++)

	{

		for(int j=0;j<n;j++)

		{

			ans+=mat.x[i][j]%mod;

			ans%=mod;

			//printf("%lld ",mat.x[i][j]);

		}

		//printf("\n");

	}

	printf("%lld\n",ans);

	return 0;

}

// Faster C++ solution for the same problem (measured speedup: 1.4976965307283907x):
#include <iostream>

#include <vector>



using namespace std;



	

struct Matrix{

  static const int maxn = 105;

  long long int A[maxn][maxn];

  int n;

 

  Matrix(int n_): n(n_) {

    for(int i = 0; i < n; ++i)

      for(int j = 0; j < n; ++j)

        A[i][j] = 0;

  }

 

  Matrix operator * (const Matrix& o) const {

    Matrix ret(n);

    for(int i = 0; i < n; ++i)

      for(int j = 0; j < n; ++j)

        for(int k = 0; k < n; ++k){

          long long int c = A[i][k] * o.A[k][j] % 1000000007;

          ret.A[i][j] = (ret.A[i][j] + c) % 1000000007;

        }

    return ret;

  }

 

  static Matrix identity(int n){

    Matrix ret(n);

    for(int i = 0; i < n; ++i)

      ret.A[i][i] = 1;

    return ret;

  }

 

  Matrix operator ^ (long long k) const {

    Matrix ret = identity(n);

    Matrix a = *this;

    while(k){

      if(k & 1)ret = ret * a;

      a = a * a;

      k /= 2;

    }

    return ret;

  }

 

};



int main(void) {

    long long int N, K;

    cin >> N >> K;

    Matrix matrix(N);

    for(int i=0; i<N; i++) {

        for(int j=0; j<N; j++) {

            cin >> matrix.A[i][j];

        }

    }

    // for(int k=0; k<K-1; k++) {

    //     Matrix matrix_tmp(N, vector<int>(N));

    //     for(int i=0; i<N; i++) {

    //         for(int j=0; j<N; j++) {

    //             for(int l=0; l<N; l++) {

    //                 matrix_tmp[i][j] += matrix[i][l] *matrix[l][j];

    //             }

    //         }

    //     }

    //     matrix = matrix_tmp;

    // }

    // for(int i=0; i<N; i++) {

    //     for(int j=0; j<N; j++) {

    //         cout << matrix[i][j] << " ";

    //     }

    //     cout << endl;

    // }

    matrix = matrix^(K);

    int ans = 0;

    for(int i=0; i<N; i++) {

        for(int j=0; j<N; j++) {

            ans = (ans + matrix.A[i][j]) % 1000000007;

        }

    }

    cout << ans << endl;

    return 0;

}