///////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////
//                       _oo0oo_                         //
//                      o8888888o                        //
//                      88" . "88      ------ hzt1       //
//                      (| -_- |)                        //
//                      0\  =  /0                        //
//                    ___/`---'\___                      //
//                  .' \|     |// '.                     //
//                 / \|||  :  |||// \                    //
//                / _||||| -:- |||||- \                  //
//               |   | \  -  /// |     |                 //
//               | \_|  ''\---/''  |_/ |                 //
//               \  .-\__  '-'  ___/-. /                 //
//             ___'. .'  /--.--\  `. .'___               //
//          ."" '<  `.___\_<|>_/___.' >' "".             //
//         | | :  `- \`.;`\ _ /`;.`/ - ` : | |           //
//         \  \ `_.   \_ __\ /__ _/   .-` /  /           //
//     =====`-.____`.___ \_____/___.-`___.-'=====        //
//                       `=---='                         //
//                                                       //
//                                                       //
//     ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~       //
//                                                       //
//                 God-He Bless All.                     //
//           This Code Will Never Explode.               //
//                                                       //
//                                                       //
///////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////
#include<cstdio>
#include<cstring>
#include<cctype>
#include<algorithm>
#define rep(i,a,b) for(int i=(a),i##_end_=(b);i<=i##_end_;i++)
#define dwn(i,a,b) for(int i=(a),i##_end_=(b);i>=i##_end_;i--)
using namespace std;
const int Size=1<<16;
char buffer[Size],*head,*tail;
inline char Getchar() {
    if(head==tail) {
        int l=fread(buffer,1,Size,stdin);
        tail=(head=buffer)+l;
    }
    if(head==tail) return -1;
    return *head++;
}
inline int read() {
    int x=0,f=1;char c=getchar();
    for(;!isdigit(c);c=getchar()) if(c=='-') f=-1;
    for(;isdigit(c);c=getchar()) x=x*10+c-'0';
    return x*f;
}
typedef long long ll;
const int maxn=2010;
const int maxm=4000010;
const int mod=1000000007;
int fac[maxm],ifac[maxm];
int C(int a,int b) {if(a<b) return 0;return (ll)fac[a]*ifac[b]%mod*ifac[a-b]%mod;}
void init(int n) {
	ifac[0]=fac[0]=ifac[1]=1;
	rep(i,2,n) ifac[i]=(ll)(mod-mod/i)*ifac[mod%i]%mod;
	rep(i,1,n) fac[i]=(ll)fac[i-1]*i%mod,ifac[i]=(ll)ifac[i-1]*ifac[i]%mod;
}
int f[maxn*2][maxn];
int main() {
	int n=read(),k=read();init(n*k);
	if(k==1) {puts("1");return 0;}
	f[2*n+1][0]=1;
	dwn(i,2*n+1,2) rep(j,0,2*n-i+1) {
		int& ans=f[i][j];if(!ans) continue;
		int a=(2*n-i+1+j)/2,b=(2*n-i+1-j)/2;
		if(j) (f[i-1][j-1]+=ans)%=mod;
		if(a<n) (f[i-1][j+1]+=(ll)ans*C((k-1)*(a+1)+b-1,k-2)%mod*(n-a)%mod)%=mod;
	}
	printf("%d\n",f[1][0]);
	return 0;
}
/*
a + b = n-i+1
a - b = j

a = (n-i+1+j)/2
b = (n-i+1-j)/2

(k-1)*a + b + (k-2)
*/


Measured memory_bytes: 62976