#include <cstdio>
#include <iostream>
#include <algorithm>
using namespace std;
const int M = 105;
const int MOD = 1e9+7;
#define int long long
int read()
{
	int x=0,f=1;char c;
	while((c=getchar())<'0' || c>'9') {if(c=='-') f=-1;}
	while(c>='0' && c<='9') {x=(x<<3)+(x<<1)+(c^48);c=getchar();}
	return x*f;
}
int n,c,k,ans,sub,h[M],f[M][M],g[M][M];
pair<int,int> b[M*M];
void add(int &x,int y) {x=(x+y)%MOD;}
signed main()
{
	n=read();c=read();k=read();
	for(int i=1;i<=n;i++) for(int j=1;j<=k;j++)
	{
		int x=read();add(ans,x);
		b[(i-1)*k+j]=make_pair(x,i);
	}
	sort(b+1,b+1+n*k);
	for(int w=1;w<=n*k;w++)
	{
		for(int i=0;i<=n;i++) for(int j=0;j<=n;j++)
			f[i][j]=g[i][j]=0;
		f[0][c]=1;
		for(int i=1;i<=n;i++) for(int j=0;j<=c;j++)
		for(int t=0;t<2;t++)
		{
			int vl=t?h[i]:k-h[i],tj=j+2*t;
			int f1=f[i-1][j]*vl%MOD,f2=g[i-1][j]*vl%MOD;
			if(tj) tj--;
			else f2=(f2+f1)%MOD;
			tj=min(tj,c);
			add(f[i][tj],f1);add(g[i][tj],f2);
		}
		for(int i=0;i<=n;i++)
			add(sub,g[n][i]*(b[w].first-b[w-1].first));
		h[b[w].second]++;
	}
	for(int i=1;i<n;i++) ans=ans*k%MOD;
	printf("%lld\n",(ans-sub+MOD)%MOD);
}