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

using ll=long long;
using pii=pair<int,int>;
#define pb push_back
#define all(a) a.begin(),a.end()
#define sz(a) ((int)a.size())

const int N=2005,M=1000005,mod=1000000007;
int add(int x, int y){x+=y; if(x>=mod) x-=mod; return x;}
int sub(int x, int y){x-=y; if(x<0) x+=mod; return x;}
int mul(int x, int y){return 1ll*x*y%mod;}
int Pow(int x, ll y=mod-2){int res=1; for(; y; y>>=1,x=mul(x,x)) if(y&1) res=mul(res,x); return res;}

int n,m,k,pw[M],dp[N][N],fac[N],inv[N],ifac[N],kval[N],inv2[N];

signed main(){
    ios_base::sync_with_stdio(0),cin.tie(0);
    fac[0]=inv[1]=ifac[0]=1;
    for(int i=1; i<N; ++i) fac[i]=mul(fac[i-1],i);
    for(int i=2; i<N; ++i) inv[i]=mul(inv[mod%i],mod-mod/i);
    for(int i=1; i<N; ++i) ifac[i]=mul(ifac[i-1],inv[i]);
    inv2[0]=1,inv2[1]=Pow(2);
    for(int i=2; i<N; ++i) inv2[i]=mul(inv2[i-1],inv2[1]);
    cin >> n >> m >> k;
    pw[0]=1;
    for(int i=1; i<M; ++i) pw[i]=mul(pw[i-1],k);
    kval[0]=1;
    for(int i=1; i<N; ++i) kval[i]=mul(kval[i-1],k+1-i);
    int tot=0;
    for(int i=0; i*2<=m; ++i){
        int tmp=mul(mul(kval[m-i],fac[m]),mul(ifac[i],ifac[m-i*2]));
        tmp=mul(tmp,inv2[i]);
        tot=add(tot,tmp);
        //cout << tot << "\n";
    }

    int res=0;
    for(int x=1; x*2<=m; ++x){
        int de=0;
        for(int i=x; i*2<=m; ++i){
            int tmp=mul(mul(kval[m-i],fac[m-x*2]),mul(ifac[i-x],ifac[m-i*2]));
            //tmp=mul(mul(tmp,fac[i]),mul(ifac[i-x],ifac[x]));
            tmp=mul(tmp,inv2[i-x]);
            de=add(de,tmp);
        }
        //cout << x << ' ' << de << "\n";
        tot=sub(tot,de);
        for(int i=1; (i-1)*(m-x)+m<=n; ++i){
            if(i&1) res=add(res,mul(de,mul(pw[n-((i-1)*(m-x)+m)],n-((i-1)*(m-x)+m)+1)));
            else res=sub(res,mul(de,mul(pw[n-((i-1)*(m-x)+m)],n-((i-1)*(m-x)+m)+1)));
        }
    }
    res=add(res,mul(tot,mul(pw[n-m],n-m+1)));
    cout << res << "\n";
}