//Author: Kevin
#include<bits/stdc++.h>
//#pragma GCC optimize("O2")
using namespace std;
#define ll long long
#define ull ll
#define pb emplace_back
#define mp make_pair
#define ALL(x) (x).begin(),(x).end()
#define rALL(x) (x).rbegin(),(x).rend()
#define srt(x) sort(ALL(x))
#define rev(x) reverse(ALL(x))
#define rsrt(x) sort(rALL(x))
#define sz(x) (int)(x.size())
#define inf 0x3f3f3f3f
#define pii pair<int,int>
#define lb(v,x) (int)(lower_bound(ALL(v),x)-v.begin())
#define ub(v,x) (int)(upper_bound(ALL(v),x)-v.begin())
#define uni(v) v.resize(unique(ALL(v))-v.begin())
#define longer __int128_t
void die(string S){puts(S.c_str());exit(0);}
const ull mod=1e9+7;
int n,m;
ull A[6][6],B[6][6];
int Lim[500500];
ull Val[6][6][1001000],nVal[6][6][1001000];
ull pw2[2001000];
void solve(int m,int Lim)
{
	if(!m)
	{
		for(int i=0;i<n;i++)
			for(int j=0;j<n;j++)
				for(int x=0;x<Lim;x++)
					Val[i][j][x]=(i==j);
		return ;
	}
	solve(m/2,Lim+m/2);
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
			memset(nVal[i][j],0,sizeof(ull)*Lim);
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
			for(int k=0;k<n;k++)
				for(int x=0;x<Lim;x++)
					nVal[i][k][x]=(nVal[i][k][x]+Val[i][j][x]*Val[j][k][x+m/2])%mod;
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
			memcpy(Val[i][j],nVal[i][j],sizeof(ull)*Lim);
	if(m&1)
	{
		for(int i=0;i<n;i++)
			for(int j=0;j<n;j++)
				memset(nVal[i][j],0,sizeof(ull)*Lim);
		for(int i=0;i<n;i++)
			for(int j=0;j<n;j++)
				for(int k=0;k<n;k++)
					for(int x=0;x<Lim;x++)
						nVal[i][k][x]=(nVal[i][k][x]+(A[j][k]*pw2[x+m-1]+B[j][k])%mod*Val[i][j][x])%mod;
		for(int i=0;i<n;i++)
			for(int j=0;j<n;j++)
				memcpy(Val[i][j],nVal[i][j],sizeof(ull)*Lim);
	}
}
ll tot[500500],tmp[500500];
ll ksm(ll a,ll b)
{
	ll ans=1;
	while(b)
	{
		if(b&1) ans=ans*a%mod;
		b>>=1;
		a=a*a%mod;
	}
	return ans;
}
int main()
{
	ios_base::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	cin>>n>>m;
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
			cin>>A[i][j];
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
			cin>>B[i][j];
	pw2[0]=1;
	for(int i=1;i<2001000;i++)
		pw2[i]=pw2[i-1]*2%mod;
	solve(m,m+1);
	for(int i=0;i<=m;i++)
	{
		memset(tmp,0,sizeof(ll)*(m+1));
		tmp[0]=1;
		for(int j=0;j<=m;j++)
			if(i!=j)
			{
				for(int k=m;k;k--)
					tmp[k]=tmp[k-1];
				tmp[0]=0;
				ll val=(mod-pw2[j])%mod;
				for(int k=1;k<=m;k++)
					tmp[k-1]=(tmp[k-1]+tmp[k]*val)%mod;
			}
		ll val2=1;
		for(int j=0;j<=m;j++)
			if(i!=j)
				val2=val2*(pw2[i]-pw2[j]+mod)%mod;
		val2=ksm(val2,mod-2);
		for(int j=0;j<=m;j++)
			tot[j]=(tot[j]+tmp[j]*val2%mod*Val[0][0][i])%mod;
	}
	for(int j=0;j<=m;j++)
		cout<<tot[j]<<'\n';
	return 0;
}