#include<cstdio>

#include<cstring>

#include<algorithm>

#define ll long long

using namespace std;

const int N=,M=N*N,mod=1e9+;

int upt(int x)

{while(x>=mod)x-=mod;while(x<)x+=mod;return x;}

int pw(int x,int k)

{int ret=;while(k){if(k&)ret=(ll)ret*x%mod;x=(ll)x*x%mod;k>>=;}return ret;}

int n,k,jc[M],jcn[M],dp[N][N];

void init()

{

  int lm=n*k;

  jc[]=;for(int i=;i<=lm;i++)jc[i]=(ll)jc[i-]*i%mod;

  jcn[lm]=pw(jc[lm],mod-);

  for(int i=lm-;i>=;i--)jcn[i]=(ll)jcn[i+]*(i+)%mod;

}

int C(int n,int m)

{

  if(n<||m<||n<m)return ;

  return (ll)jc[n]*jcn[m]%mod*jcn[n-m]%mod;

}

int main()

{

  scanf("%d%d",&n,&k); if(k==){puts("");return ;}

  init();

  dp[][]=;

  for(int i=;i<=n;i++)

    for(int j=;j<=i;j++)

      {

    if(!dp[i][j])continue;int tp=dp[i][j];

    if(i<n)dp[i+][j]=upt(dp[i+][j]+tp);

    if(j<i)

      {

        int ml=C((n-j)*(k-)+n-i-,k-);

        dp[i][j+]=(dp[i][j+]+(ll)ml*tp)%mod;

      }

      }

  printf("%lld\n",(ll)dp[n][n]*jc[n]%mod);

  return ;

}