#include <bits/stdc++.h>
using namespace std;

#define int long long

#define MOD 1000000007

int inv(int x){
	int p=MOD-2;
	int r=1;
	for(int p=MOD-2;p;p>>=1){
		if(p&1){
			r*=x;
			r%=MOD;
		}
		x*=x;
		x%=MOD;
	}
	return r;
}

signed main(){
	ios::sync_with_stdio(false);
	cin.tie(nullptr);
	int N,K,L;
	cin>>N>>K>>L;
	int X[K+1][L+1];
	for(int i=1;i<=K;i++){
		for(int j=1;j<=L;j++){
			cin>>X[i][j];
		}
	}
	int F[N+1];
	F[0]=1;
	for(int i=1;i<=N;i++){
		F[i]=F[i-1]*i%MOD;
	}
	int FI[N+1];
	for(int i=0;i<=N;i++){
		FI[i]=inv(F[i]);
	}
	int dp[N+1][N+1];
	memset(dp,0,sizeof(dp));
	for(int i=0;i<=N;i++){
		dp[0][i]=0;
		dp[i][0]=F[i]*F[i]%MOD;
	}
	for(int i=1;i<=N;i++){
		for(int j=1;j<=i;j++){
			dp[i][j]+=dp[i-1][j]*(i-j)*(i-j);
			dp[i][j]+=dp[i-1][j-1]*j*(i-j)*2;
			if(j>=2){
				dp[i][j]+=j*(j-1)*dp[i-1][j-2];
			}
			dp[i][j]%=MOD;
		}
	}
	for(int i=0;i<=N;i++){
		for(int j=0;j<=i;j++){
			dp[i][j]*=FI[i];
			dp[i][j]%=MOD;
		}
	}
	// for(int i=1;i<=N;i++){
		// for(int j=0;j<=i;j++){
			// vector<int>v(i),w(i);
			// for(int k=0;k<j;k++){
				// v[k]=w[k]=k;
			// }
			// for(int k=j;k<i;k++){
				// v[k]=i+k;
				// w[k]=i*2+k;
			// }
			// int ans=0;
			// do{
				// int ok=1;
				// for(int k=0;k<i;k++){
					// ok&=v[k]!=w[k];
				// }
				// ans+=ok;
			// }
			// while(next_permutation(w.begin(),w.end()));
			// if(ans!=dp[i][j]){
				// cout<<i<<' '<<j<<endl;
				// cout<<"got "<<dp[i][j]<<'\n';
				// cout<<"exp "<<ans<<'\n';
				// exit(1);
			// }
		// }
	// }
	int ans=F[N-L];
	int v[N+1];
	for(int i=2;i<=K;i++){
		memset(v,0,sizeof(v));
		for(int j=1;j<=L;j++){
			v[X[i][j]]++;
			v[X[i-1][j]]++;
		}
		int y=0;
		for(int j=1;j<=N;j++){
			y+=v[j]==0;
		}
		ans*=dp[N-L][y];
		ans%=MOD;
	}
	for(int i=K+1;i<=N;i++){
		ans*=dp[N][N];
		ans%=MOD;
	}
	cout<<ans<<'\n';
	return 0;
}