// You either die a hero, or you live long enough to see yourself become the villain.

/**
 *   author:  Blind_is_love
 *   created: 21.12.2019, 17:08
**/

#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>

using namespace __gnu_pbds;
using namespace std;

typedef long long ll;

template <class T> ostream& operator << (ostream &os, const vector<T> &v) { for (T i : v) os << i << ' '; return os; }
template <class T> ostream& operator << (ostream &os, const set<T> &v) { for (T i : v) os << i << ' '; return os; }
template <class T, class S> ostream& operator << (ostream &os, const pair<T, S> &v) { os << v.first << ' ' << v.second; return os; }
template <class T, class S> ostream& operator << (ostream &os, const map<T, S> &v) { for (auto i : v) os << '(' << i.first << "=>" << i.second << ')' << ' '; return os; }
template <class T, class S> ostream& operator << (ostream &os, const unordered_map<T, S> &v) { for (auto i : v) os << '(' << i.first << "=>" << i.second << ')' << ' '; return os; }
 
#ifndef ONLINE_JUDGE
#define trace(...) __f(#__VA_ARGS__, __VA_ARGS__)
    template <class Arg1> void __f(const char* name, Arg1&& arg1) { cerr << name << " : " << arg1 << endl; }
    template <class Arg1, class... Args>
    void __f(const char* names, Arg1&& arg1, Args&&... args) {
        const char* sep = strchr(names + 1, ',');
        cerr.write(names, sep - names) << " : " << arg1 << "  ";
        __f(sep + 1, args...);
    }
#else
#define trace(...) 0
#pragma GCC optimize ("O3")
#pragma GCC optimize ("unroll-loops")
#pragma GCC target("avx2,sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")
#define _CRT_SECURE_NO_WARNINGS
#endif // ifndef ONLINE_JUDGE

#define endl "\n"
#define vll vector<int>
#define max_heap priority_queue<int>
#define min_heap priority_queue<int , vector<int>, greater<int> >
#define pb push_back
#define all(container) container.begin(), container.end()
#define pll pair<ll,ll>
#define pii pair<int,int>
#define printArr(name,from,to) for(int x = from; x < to; x++) {cout << name[x] << " ";} cout << endl;
#define fi first
#define se second
#define mp make_pair
#define prime1 304933
#define prime2 15486277 
#define mod2 179424691
#define int ll

typedef tree<ll, null_type, less<ll>, rb_tree_tag,tree_order_statistics_node_update> ordered_set; 

const ll MOD = 1000000007;
const ll MAX = 51;
const double PI = 3.1415926;
const ll INF = 1e18;
const double EPS = 1e-5;
int matrix[MAX][MAX];
int n, k;
int dp[64][MAX][MAX];


void calc(int s){
	for(int i = 0; i < n; ++i) {
		for(int j = 0; j < n; ++j) {
			for (int k = 0; k < n; ++k) {
				dp[s][i][j] += (dp[s - 1][i][k] * dp[s - 1][k][j]);
				dp[s][i][j] %= MOD;
			}
		}
	}
}

void precomp() {
	memset(dp, 0, sizeof(dp));
	
	for(int i = 0; i < n; ++i) {
		for(int j = 0; j < n; ++j) {
			dp[0][i][j] = matrix[i][j];
		}
	}

	for(int i = 1; i < 64; ++i) {
		calc(i);
	}
}

signed main() {
	ios::sync_with_stdio(false); cin.tie(0); cout.tie(0);
	clock_t t1, t2;
	t1 = clock();
	cin >> n >> k;
	// cout << endl;
	for(int i = 0; i < n; ++i) {
		for(int j = 0; j < n; ++j) {
			cin >> matrix[i][j];
		}
	}
	precomp();
	int ans[n][n];

	for(int i = 0; i < n; ++i) {
		for(int j = 0; j < n; ++j) {
			ans[i][j] = (i == j) ? 1 : 0;
		}
	}

	for(int i = 0; i < 64; ++i) {
		if ((1LL<<i) & k) {
			int temp[n][n];
			memset(temp, 0, sizeof(temp));
			for(int h = 0; h < n; ++h) {
				for(int j = 0; j < n; ++j) {
					for(int k = 0; k < n; ++k) {
						temp[h][j] += (ans[h][k] * dp[i][k][j]);
						temp[h][j] %= MOD;
					}
				}
			}
			for(int j = 0; j < n; ++j){
				for(int k = 0; k < n; ++k) {
					ans[j][k] = temp[j][k];
				}
			}
		}
	}
	int out = 0;
	for(int i = 0; i < n; ++i) {
		for(int j = 0; j < n; ++j) {
			out += (ans[i][j]);
			out %= MOD;
		}
	}
	cout << out << endl;
	t2 = clock(); cerr << "time taken: " << (t2-t1)/(long double)CLOCKS_PER_SEC << endl;
	return 0;
}


Measured memory_bytes: 1664