#include<bits/stdc++.h>
#define ll long long
#define MAX 2005
using namespace std;
const ll mod=1e9+7;
ll qmod(ll x){ if(x>=mod) x-=mod;return x; }
ll fpow(ll a,ll b){
	ll ret=1;
	while(b){
		if(b&1) ret=ret*a%mod;
		a=a*a%mod;b>>=1;
	}return ret;
}ll inv(ll x){ return fpow(x,mod-2); }
ll fac[MAX*100],ivf[MAX*100];
ll C(ll n,ll m){
	if(n<m||n<0||m<0) return 0;
	return fac[n]*ivf[m]%mod*ivf[n-m]%mod;
}
int n,m,B;ll a[MAX*100];
ll f1[MAX][MAX],f2[50][MAX],f[MAX],g[MAX];
int main(){
	scanf("%d%d",&n,&m);B=sqrt(n);fac[0]=1;
	for(int i=1;i<=n*m;i++){ scanf("%d",&a[i]),fac[i]=fac[i-1]*i%mod; }
	ivf[n*m]=inv(fac[n*m]);for(int i=n*m-1;i>=0;i--) ivf[i]=ivf[i+1]*(i+1)%mod;
	f1[0][0]=1;
	for(int i=1;i<=n;i++) for(int j=1;j<=n;j++){
		for(int k=0;k<=min(B,j);k++)
			f1[i][j]=qmod(f1[i][j]+f1[i-1][j-k]*ivf[k*m-1]%mod);
//		cout<<i<<" "<<j<<" "<<f1[i][j]*fac[j*m-i]%mod<<endl;
	}		
	f2[0][0]=1;
	for(int i=1;i<=n/(B+1);i++) for(int j=B+1;j<=n;j++){
		for(int k=B+1;k<=j;k++)
			f2[i][j]=qmod(f2[i][j]+f2[i-1][j-k]*ivf[k*m-1]%mod);
	} 
	for(int i=1;i<=n;i++) for(int j=0;j<=min(i,n/(B+1));j++)
		for(int k=0;k<=n;k++){
			if(k*m-j>=0&&(n-k)*m-(i-j)>=0){
//				cout<<i<<" "<<j<<" "<<k<<" "<<f2[j][k]*fac[k*m-j]%mod<<" "<<f1[i-j][n-k]%mod*fac[(n-k)*m-(i-j)]%mod<<" "<<C(n*m-i,k*m-j)<<endl;			
				f[i]=qmod(f[i]+f2[j][k]*fac[k*m-j]%mod*f1[i-j][n-k]%mod*fac[(n-k)*m-(i-j)]%mod*C(n*m-i,k*m-j)%mod*C(i,j)%mod);
			}
		}
	g[0]=1;
	for(int i=1;i<=n*m;i++)
		for(int j=n;j>=0;j--)
			if(j) g[j]=qmod(g[j]+g[j-1]*a[i]%mod);
//	for(int i=1;i<=n;i++) cout<<g[i]<<" "<<f[i]<<endl;
	ll ans=0;
	for(int i=1;i<=n;i++) ans=qmod(ans+g[i]*f[i]%mod);
	printf("%lld",ans);
}