#include <bits/stdc++.h>
using namespace std;
#define int long long
const int N=1505,M=105,T=40,G=1.5e5+5,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][N],g[N][N],rs[N][N],dp[2][G],n,m,jc[G],ny[G],p[G],res;
int C(int a,int b){
	return jc[a]*ny[b]%mods*ny[a-b]%mods;
}
signed main(){
	cin>>n>>m;
	n*=m;
	jc[0]=ny[0]=1;
	for(int i=1;i<=n;i++)jc[i]=jc[i-1]*i%mods,ny[i]=pows(jc[i],mods-2);
	for(int i=1;i<=n;i++)cin>>p[i];
	f[0][0]=g[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-(j-1)*(m-1)&&k<m-1+T*m;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;
			}
		}
	}
	for(int j=1;j<=n/m&&j<=T;j++){
		for(int i=j*(m-1);i<=n;i+=m){
			for(int k=m-1+T*m;k<=i-(j-1)*(m-1);k+=m){
				g[(i-j*(m-1))/m][j]+=g[(i-k-(j-1)*(m-1))/m][j-1]*C(i,k);
				g[(i-j*(m-1))/m][j]%=mods;
			}
		}
	}
	for(int j=1;j<=n/m;j++){
		for(int a=0;a<=j&&a<=T;a++){
			for(int b=a*(m-1);b<=n-(j-a)*(m-1);b+=m){
				rs[(n-j-j*(m-1))/m][j]+=g[(b-a*(m-1))/m][a]*f[(n-j-b-(j-a)*(m-1))/m][j-a]%mods*C(j,a)%mods*C(n-j,b);
				rs[(n-j-j*(m-1))/m][j]%=mods;
			}
		}
	}
	for(int i=0;i<=n;i++)dp[0][i]=1;
	return 0;
	for(int i=1;i<=n/m;i++){
		dp[1][0]=0;
		for(int j=1;j<=n;j++){
			dp[1][j]=dp[1][j-1]+dp[0][j-1]*p[j];
			dp[1][j]%=mods;
		}
		res+=dp[1][n]*rs[n/m-i][i];
		res%=mods;
		swap(dp[1],dp[0]);
	}
	cout<<res;
}