#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];
	for(int i=2;i<=K;i++){
		set<int>s;
		for(int j=1;j<=L;j++){
			s.insert(X[i][j]);
			s.insert(X[i-1][j]);
		}
		ans*=dp[N-L][L*2-s.size()];
		ans%=MOD;
	}
	for(int i=K+1;i<=N;i++){
		ans*=dp[N][N];
		ans%=MOD;
	}
	cout<<ans<<'\n';
	return 0;
}