#include<bits/stdc++.h>
using namespace std;

const int N=1505,M=1.5e5+5;
const int mod=1e9+7;

int n,m,ans;
int a[M],f[N],g[N][N],h[N][N],fact[M],infact[M];

int qpow(int a,int n)
{
	int ret=1;
	while(n)
	{
		if(n&1)ret=1ll*ret*a%mod;
		a=1ll*a*a%mod;
		n>>=1;
	}
	return ret;
}

void prepare()
{
	fact[0]=1;
	for(int i=1;i<M;i++)fact[i]=1ll*fact[i-1]*i%mod;
	infact[M-1]=qpow(fact[M-1],mod-2);
	for(int i=M-2;i>=0;i--)infact[i]=1ll*infact[i+1]*(i+1)%mod;
}

int C(int n,int m){return 1ll*fact[n]*infact[m]%mod*infact[n-m]%mod;}

int main()
{
	prepare();
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n*m;i++)scanf("%d",&a[i]);
	f[0]=1;
	for(int i=1;i<=n*m;i++)
		for(int j=min(i,n);j>=1;j--)
			f[j]=(f[j]+1ll*a[i]*f[j-1])%mod;
	if(m==1)
	{
		for(int i=1;i<=n;i++)ans=(ans+1ll*qpow(i,n-i)*f[i])%mod;
		cout<<ans;
		return 0;
	}
	int B=sqrt(n);
	g[0][0]=1;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			for(int k=1;k<=B;k++)
				g[i][j]=(g[i][j]+1ll*g[i-k][j-1]*infact[k*m-1])%mod;
	h[0][0]=1;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n/B;j++)
			for(int k=B+1;k<=n;k++)
				h[i][j]=(h[i][j]+1ll*h[i-k][j-1]*infact[k*m-1])%mod;
	for(int i=1;i<=n;i++)
	{
		int sum=0;
		for(int j=0;j<=n;j++)
			for(int k=0;k<=i&&k<=n/B;k++)
				sum=(sum+1ll*h[j][k]*g[n-j][i-k]%mod*C(i,k))%mod;
		ans=(ans+1ll*sum*fact[n*m-i]%mod*f[i])%mod;
	}
	cout<<ans;
	return 0;
}