This submission is migrated from previous version of oj.uz, which used different machine for grading. This submission may have different result if resubmitted.
#include <iostream>
#define MOD 1000000007
#define MAXN 10000
using namespace std;
long long dp[MAXN + 1][2];
int val[MAXN];
long long fact[MAXN + 1];
long long invFact[MAXN + 1];
int m, n;
long long makeInv(long long x) {
  int pow;
  long long ans;
  ans = 1;
  pow = MOD - 2;
  while ( pow ) {
    if ( pow % 2 ) {
      ans = (ans * x) % MOD;
    }
    pow /= 2;
    x = (x * x) % MOD;
  }
  return ans;
}
void makeFact() {
  long long ans;
  int i;
  ans = 1;
  fact[0] = invFact[0] = 1;
  for ( i = 1; i <= MAXN; i++ ) {
    ans = (ans * i) % MOD;
    fact[i] = ans;
    invFact[i] = makeInv(ans);
  }
}
long long makeCombi(int n, int choose) {
  if ( n >= choose && choose >= 0 ) {
    return (((fact[n] * invFact[choose]) % MOD) * invFact[n - choose]) % MOD;
  }
  return 0;
}
long long makeChoose(int n, int choose) {
  if ( n >= choose && choose >= 0 ) {
    return (fact[n] * invFact[n - choose]) % MOD;
  }
  return 0;
}
long long calcAns(int a, int b, int x) { ///from a to b(x from b have pair); we choose from a increasing, other does not matter
  if ( a >= x && x >= 0 && b >= x ) {
    return (makeCombi(a, x) * makeChoose(b, x)) % MOD;
  }
  return 0;
}
void smh() {
  int currNr, i, pos, i2, i3;
  //makeFact();
  //scanf("%d%d", &n, &m);
  n *= 2;
  for ( i = 0; i < n; i++ ) {
    val[i % m]++;
  }
  pos = 0;
  dp[0][0] = 1;
  currNr = 0;
  for ( i = 0; i < m; i++ ) {
    pos ^= 1;
    for ( i2 = 0; i2 <= MAXN; i2++ ) {
      dp[i2][pos] = 0;
    }
    for ( i2 = 0; i2 <= val[i]; i2++ ) { //how many we pair
      for ( i3 = 0; i3 <= currNr; i3++ ) { //how many are used
        dp[i3 + 2 * i2][pos] = (dp[i3 + 2 * i2][pos] + (dp[i3][pos ^ 1] * calcAns(currNr - i3, val[i], i2)) % MOD) % MOD;
      }
    }
    currNr += val[i];
  }
  long long pow;
  pow = 1;
  for ( i = 1; i <= n / 2; i++ ) {
    pow = ((long long) pow * 2) % MOD;
  }
  printf("%lld\n", (((dp[n][pos] * fact[n / 2]) % MOD) * pow) % MOD);
}
int main() {
  makeFact();
  scanf("%d%d", &n, &m);
  smh();
}
Compilation message (stderr)
Main.cpp: In function 'int main()':
Main.cpp:109:8: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
  109 |   scanf("%d%d", &n, &m);
      |   ~~~~~^~~~~~~~~~~~~~~~| # | Verdict  | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict  | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict  | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict  | Execution time | Memory | Grader output | 
|---|
| Fetching results... | 
| # | Verdict  | Execution time | Memory | Grader output | 
|---|
| Fetching results... |