#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[50005][1505],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;
	}
	g[0][0]=1;
	for(int i=1;i<=n*m;i++)
		for(int j=1;j<=n;j++)
			for(int k=m-1;k<=i;k+=m)g[i][j]=(g[i][j]+1ll*g[i-k][j-1]*C(i,k))%mod;
	for(int i=1;i<=n;i++)ans=(ans+1ll*g[n*m-i][i]*f[i]%mod)%mod;
	cout<<ans;
	return 0;
}