#include<bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
int n,m,k,f[2005],g[2005],ans,pw[2005],ipw[2005],pwk[1000005];
int fac[2005],inv[2005],dn[2005][2005];
int qpow(int a,int b)
{
	if(b==0) return 1;
	int g=qpow(a,b/2);
	g=1ll*g*g%mod;
	if(b&1) g=1ll*g*a%mod;
	return g;
}
int calc(int x)
{
	if(x>n) return 0;
	return 1ll*(n-x+1)*pwk[n-x]%mod;
}
int C(int n,int m)
{
	if(n<m) return 0;
	return 1ll*dn[k-n+1][m]*inv[m]%mod;
}
int main()
{
	cin>>n>>m>>k;
	fac[0]=inv[0]=pw[0]=ipw[0]=pwk[0]=1;
	for(int i=1;i<=m;i++)
	{
		fac[i]=1ll*fac[i-1]*i%mod;
		inv[i]=qpow(fac[i],mod-2);
		pw[i]=1ll*pw[i-1]*2%mod;
		ipw[i]=qpow(pw[i],mod-2);
	}
	for(int i=1;i<=n;i++)
	{
		pwk[i]=1ll*pwk[i-1]*k%mod;
	}
	for(int i=1;i<=m;i++)
	{
		dn[i][0]=1;
		dn[i][1]=(k-i+1);
		for(int j=1;j<=m;j++)
		{
			dn[i][j+1]=1ll*dn[i][j]*(k-i+1-j)%mod;
		}
	}
	for(int i=0;i*2<=m;i++)
	{
		if(i==0) g[i]=calc(m);
		else
		{
			int c=1;
			for(int j=m;j<=n;j+=m-i)
			{
				g[i]+=1ll*c*calc(j)%mod;
				if(g[i]>=mod) g[i]-=mod;
				c=mod-c;
			}
		}
		int w=1ll*C(k,i)*fac[i]%mod,sum=0;
		for(int j=0;i*2+j*2<=m;j++)
		{
			sum+=1ll*C(k-i,j)*C(k-i-j,m-i*2-j*2)%mod*fac[m-i*2]%mod*ipw[j]%mod;
			if(sum>=mod) sum-=mod;
		}
		f[i]=1ll*w*sum%mod;
	}
	for(int i=1;i*2<=m;i++)
	{
		f[0]-=f[i];
		if(f[0]<0) f[0]+=mod;
	}
	for(int i=0;i*2<=m;i++)
	{
		ans+=1ll*f[i]*g[i]%mod;
		if(ans>=mod) ans-=mod;
	}
	cout<<ans;
}