#include <bits/stdc++.h>
using namespace std;
#define int long long
const int N=1505,M=105,mods=1e9+7;
int pows(int a,int b){
	if(b==0)return 1;
	int res=pows(a,b>>1);
	res=res*res%mods;
	if(b&1)res=res*a%mods;
	return res;
}
int f[N][M],n,m,C[N][N],cnt[N],p[N],res;
signed main(){
	cin>>n>>m;
	n*=m;
	for(int i=1;i<=n;i++)cin>>p[i];
	for(int i=0;i<=n;i++)C[i][0]=1;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			C[i][j]=(C[i-1][j-1]+C[i-1][j])%mods;
		}
	}
	f[0][0]=1;
	for(int j=1;j<=n/m;j++){
		for(int i=j*(m-1);i<=n;i+=m){
			for(int k=m-1;k<=i;k+=m){
				f[(i-j*(m-1))/m][j]+=f[(i-k-(j-1)*(m-1))/m][j-1]*C[i][k];
				f[(i-j*(m-1))/m][j]%=mods;
			}
		//	cout<<i<<" "<<j<<" "<<f[i][j]<<endl;
		}
	}
	for(int i=1;i<(1<<n);i++){
		cnt[i]=cnt[i&-i^i]+1;
		int sm=1;
		for(int j=1;j<=n;j++){
			if(i>>j-1&1^1)continue;
			sm=sm*p[j]%mods;
		}
		if((n-cnt[i]-cnt[i]*(m-1))%m==0)res+=sm*f[(n-cnt[i]-cnt[i]*(m-1))/m][cnt[i]];
		res%=mods;
	//	cout<<i<<" "<<sm<<" "<<cnt[i]<<" "<<res<<endl;
	}
	cout<<res;
}