#include<bits/stdc++.h>
#define Gc() getchar()
#define Me(x,y) memset(x,y,sizeof(x))
#define Mc(x,y) memcpy(x,y,sizeof(x))
#define d(x,y) ((m)*(x-1)+(y))
#define R(n) (rnd()%(n)+1)
#define Pc(x) putchar(x)
#define LB lower_bound
#define UB upper_bound
#define fi first
#define se second
#define eb emplace_back
#define all(x) x.begin(),x.end()
using namespace std;using ll=long long;using db=double;using lb=long db;using ui=unsigned;using ull=unsigned long long;using pii=pair<int,int>;
const int N=5e5+5,M=(1<<16)+5,K=1000+5,mod=1e9+7,Mod=mod-1;const db eps=1e-9;const int INF=1e9+7;mt19937 rnd(263082);
#define Tp template<typename T>
#define Ts template<typename T,typename... Ar>
namespace Debug{
	Tp void _debug(char* f,T t){cerr<<f<<'='<<t<<endl;}
	Ts void _debug(char* f,T x,Ar... y){while(*f!=',') cerr<<*f++;cerr<<'='<<x<<",";_debug(f+1,y...);}
	#ifdef LOCAL
	#define gdb(...) _debug((char*)#__VA_ARGS__,__VA_ARGS__)
	#else 
	#define gdb(...) void()
	#endif
}using namespace Debug;


const int P=1e9+7;
using polyM=vector<ll>;
namespace PolyM{
	const int N=(1<<21)+5;
	using cp=complex<db>;
	int k,tr[N];cp pw[N];
	const db pi=acos(-1);const int Base=(1<<15);
	ll mpow(ll x,int y=P-2){ll ans=1;while(y) y&1&&(ans=ans*x%P),y>>=1,x=x*x%P;return ans;}
	void init(int n){
		for(k=2;k<=n;k<<=1);
		for(int i=0;i<k;i++) tr[i]=(tr[i>>1]>>1)|(i&1?k/2:0);
		for(int i=k/2;i<k;i++) pw[i]=cp(cos(2*pi/k*(i-k/2)),sin(2*pi/k*(i-k/2)));
		for(int i=k/2-1;i;i--) pw[i]=pw[i<<1];
	}
	void fft(cp *A,int k,int flag){
		for(int i=0;i<k;i++) if(i<tr[i]) swap(A[i],A[tr[i]]);
		for(int i=2;i<=k;i<<=1){
			cp *e=pw+i/2;
			for(int j=0;j<k;j+=i){
				for(int h=j;h<j+i/2;h++){
					cp x=A[h+i/2]*e[h-j];
					A[h+i/2]=A[h]-x;A[h]+=x;
				}
			}
		}
		if(flag) return;
		reverse(A+1,A+k);
		for(int i=0;i<k;i++) A[i]/=k;
	}
	ll round(db x){return ll(x+0.5)%P;}
	namespace Public{
		polyM operator *(const polyM &A,const polyM &B){
			static cp a[N],b[N],c[N],d[N];
			init(A.size()+B.size());
			for(int i=0;i<A.size();i++) a[i]=cp(A[i]/Base,A[i]%Base);fill(a+A.size(),a+k,0);
			for(int i=0;i<B.size();i++) b[i]=cp(B[i]/Base,B[i]%Base);fill(b+B.size(),b+k,0);
			fft(a,k,1);fft(b,k,1);
			for(int i=0;i<k;i++){
				cp x=conj(a[(k-i)%k]);
				c[i]=(x+a[i])*cp(0.5,0)*b[i];
				d[i]=(x-a[i])*cp(0,0.5)*b[i];
			}
			fft(c,k,0);fft(d,k,0);
			polyM C(A.size()+B.size()-1);
			for(int i=0;i<C.size();i++) C[i]=((round(c[i].real())*Base+round(c[i].imag()+d[i].real()))*Base+round(d[i].imag()))%P;
			return C;
		}
		polyM& operator *=(polyM &A,const polyM &B){
			return A=A*B;
		}
		polyM operator *(const polyM &A,const int &B){
			polyM C=A;for(auto &i:C) (i*=B)%=P;
			return C;
		}
		polyM& operator *=(polyM &A,const int &B){
			return A=A*B;
		}
		polyM operator +(const polyM &A,const polyM &B){
			polyM C(max(A.size(),B.size()));
			for(int i=0;i<A.size();i++) C[i]=A[i];
			for(int i=0;i<B.size();i++) (C[i]+=B[i])%=P;
			return C;
		}
		polyM& operator +=(polyM &A,const polyM &B){
			return A=A+B;
		}
		polyM operator -(const polyM &A,const polyM &B){
			polyM C(max(A.size(),B.size()));
			for(int i=0;i<A.size();i++) C[i]=A[i];
			for(int i=0;i<B.size();i++) (C[i]+=P-B[i])%=P;
			return C;
		}
		polyM& operator -=(polyM &A,const polyM &B){
			return A=A-B;
		}
		polyM operator %(const polyM &A,const int &B){
			polyM C=A;
			if(C.size()>B) C.resize(B);
			return C;
		}
		polyM& operator %=(polyM &A,const int B){
			return A=A%B;
		}
		polyM inve(polyM A,int n=-1){
			if(n==-1) n=A.size();
			if(n==1) return polyM({(int)mpow(A[0])});
			polyM A0=inve(A%(n+1>>1),n+1>>1);
			return A0*(polyM({2})-A*A0%n)%n;
		}
	}
}using namespace PolyM::Public;

namespace Pow{
	ll p1[N],p2[N];
	const int B=(1<<15);
	void init(int x){
		for(int i=p1[0]=1;i<=B;i++) p1[i]=p1[i-1]*x%mod;
		for(int i=p2[0]=1;i<=B;i++) p2[i]=p2[i-1]*p1[B]%mod;
	}
	ll qry(ll x){
		x=(x%Mod+Mod)%Mod;
		return p1[x&(1<<15)-1]*p2[x>>15]%mod;
	}
}



int n,m;
ll mpow(ll x,int y=mod-2){ll ans=1;while(y) y&1&&(ans=ans*x%mod),y>>=1,x=x*x%mod;return ans;}
using mat=array<array<ll,6>,6>;
mat operator *(const mat &A,const mat &B){
	mat C;for(int i=0;i<n;i++) C[i].fill(0);
	for(int h=0;h<n;h++) for(int i=0;i<n;i++) for(int j=0;j<n;j++) (C[i][j]+=A[i][h]*B[h][j])%=mod;
	return C;
}
mat C[M*2];
mat A,B;
polyM f;
void transdot(){
	for(int i=0;i<2*m;i++) for(int x=0;x<n;x++) for(int y=0;y<n;y++) C[i][x][y]=(A[x][y]*Pow::qry(i)+B[x][y])%mod;
	for(int i=m-1;i;i--) C[i-1]=C[i-1]*C[i];
	mat w;for(int i=0;i<n;i++) w[i].fill(0),w[i][i]=1;
	f.resize(m+1);
	for(int i=0;i<m;i++){
		f[i]=(C[i]*w)[0][0];
		w=w*C[i+m];
	}
	f[m]=w[0][0];
}
polyM g;
polyM CZT(polyM f){
	polyM g(m+1);
	for(int i=0;i<=m;i++) for(int j=0;j<=m;j++) (g[i]+=f[j]*Pow::qry(-1ll*i*j))%=mod;
	return g;
}
void Solve(){
	scanf("%d%d",&n,&m);
	for(int i=0;i<n;i++) for(int j=0;j<n;j++) scanf("%lld",&A[i][j]);
	for(int i=0;i<n;i++) for(int j=0;j<n;j++) scanf("%lld",&B[i][j]);
	Pow::init(2);
	transdot();
	ll mul=1;for(int i=1;i<=m;i++) (mul*=mod+1-Pow::qry(i))%=mod;
	mul=mpow(mul);
	for(int i=0;i<=m;i++){
		(f[i]*=mul)%=mod;
		(f[i]*=Pow::qry(-1ll*i*m))%=mod;
		(mul*=mpow(mod+1-Pow::qry(-i-1)))%=mod;
		(mul*=mod+1-Pow::qry(m-i))%=mod;
	}
	g=polyM({1});
	for(int i=0;i<=m;i++) g*=polyM({mod-Pow::qry(i),1});
	gdb(g[0]);
	for(int i=0;i<=m;i++) (f[i]*=mod-Pow::qry(-i))%=mod;
	f=CZT(f)*g;
	for(int i=0;i<=m;i++) printf("%lld\n",f[i]);
}
int main(){
	int t=1;
	// scanf("%d",&t);
	while(t--) Solve();
	cerr<<clock()*1.0/CLOCKS_PER_SEC<<'\n';
}