#include<bits/stdc++.h>
using namespace std;
using ll=long long;
const int N=1505,M=105,mod=1e9+7;
int n,m,a[N*M];
ll ans,f[N][N],F[N],g[N],fc[N*M],inv[N*M];
ll qp(ll a,int b){ll ans=1;for(;b;b>>=1,a=a*a%mod)if(b&1)ans=ans*a%mod;return ans;}
void init(int n) {
	fc[0]=1;for(int i=1;i<=n;i++)fc[i]=fc[i-1]*i%mod;
	inv[n]=qp(fc[n],mod-2);for(int i=n-1;~i;i--)inv[i]=inv[i+1]*(i+1)%mod;
}
int main() {
	freopen("study.in","r",stdin);
	freopen("study.out","w",stdout);
	ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
	cin>>n>>m,F[0]=f[0][0]=1,init(n*m);
	for(int i=1;i<=n*m;i++)cin>>a[i];
	const int T=8;
	for(int i=1;i<=n;i++)g[i]=inv[i*m-1];
	for(int i=1;i<=n;i++)for(int j=1;j<=n;j++) {
		__int128 F=0;
		for(int k=1;k<=j;k++)F+=f[i-1][j-k]*g[k];
		f[i][j]=F%mod;
	}
	for(int i=1;i<=n*m;i++) {
		const int v=a[i];
		for(int j=n;j>=T;j-=T) {
			F[j]=(F[j-1]*v+F[j])%mod;
			F[j-1]=(F[j-2]*v+F[j-1])%mod;
			F[j-2]=(F[j-3]*v+F[j-2])%mod;
			F[j-3]=(F[j-4]*v+F[j-3])%mod;
			F[j-4]=(F[j-5]*v+F[j-4])%mod;
			F[j-5]=(F[j-6]*v+F[j-5])%mod;
			F[j-6]=(F[j-7]*v+F[j-6])%mod;
			F[j-7]=(F[j-8]*v+F[j-7])%mod;
		}
		for(int j=n%T;j;j--)F[j]=(F[j-1]*v+F[j])%mod;
	}
	for(int i=1;i<=n;i++)ans=(ans+F[i]*f[i][n]%mod*fc[n*m-i])%mod;
	cout<<ans<<'\n';
	return 0;
}