#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using ull = unsigned long long;
using ld = long double;
using i128 = __int128;
using u128 = unsigned __int128;
template <typename T>
void chkmax(T &x, const T &y) {
  if (x < y) x = y;
}
template <typename T>
void chkmin(T &x, const T &y) {
  if (y < x) x = y;
}
namespace barrett {
int p;
i128 w;
void init(int _p) {
  p = _p;
  w = -1ull / p;
}
int reduce(ll x) {
  ll t = x - (w * x >> 64) * p;
  return t >= p ? t - p : t;
}
}  // namespace barrett
int mod;
void inc(int &x, int y) { x += y, x >= mod && (x -= mod); }
int mul(int x, int y) { return barrett::reduce((ll)x * y); }
int qpow(int x, int y) {
  int ret = 1;
  for (; y; y >>= 1, x = (ll)x * x % mod)
    if (y & 1) ret = (ll)ret * x % mod;
  return ret;
}
constexpr int MAXN = 2e4 + 10;
int n, f[2][MAXN][3], ext[MAXN * 2][2], fac[MAXN * 2], ifac[MAXN * 2];
int A(int x, int y) { return mul(fac[x], ifac[x - y]); }
int main() {
  ios::sync_with_stdio(false);
  cin.tie(nullptr);
  cin >> n >> mod;
  barrett::init(mod);
  fac[0] = 1;
  for (int i = 1; i <= n * 2; ++i) fac[i] = (ll)fac[i - 1] * i % mod;
  ifac[n * 2] = qpow(fac[n * 2], mod - 2);
  for (int i = n * 2 - 1; i >= 0; --i) ifac[i] = (ll)ifac[i + 1] * (i + 1) % mod;
  int p = 0;
  f[0][0][2] = 1;
  for (int i = 2; i <= n * 2; ++i) {
    p ^= 1;
    memset(f[p], 0, sizeof(int) * (min(i, n * 2 - i) + 1) * 3);
    int up = min(i - 1, n * 2 - i + 1);
    for (int j = 0; j <= up; ++j) {
      if (f[!p][j][0]) {
        inc(f[p][j + 1][0], f[!p][j][0]);
        inc(f[p][j][1], f[!p][j][0]);
      }
      if (f[!p][j][1]) inc(f[p][j][2], f[!p][j][1]);
      if (f[!p][j][2]) {
        inc(f[p][j + 1][2], f[!p][j][2]);
        if (i + j <= n * 2)
          inc(ext[i + j][1], mul(f[!p][j][2], A(n * 2 - i, j)));
        if (j && i + j - 1 <= n * 2) {
          inc(ext[i + j - 1][0], mul(mul(f[!p][j][2], j), A(n * 2 - i, j - 1)));
        }
      }
    }
    inc(f[p][0][0], ext[i][0]);
    inc(f[p][0][1], ext[i][1]);
  }
  int all = 1;
  for (int i = 1; i <= n * 2; i += 2) all = (ll)all * i % mod;
  cout << (all - f[p][0][1] + mod) % mod << "\n";
  return 0;
}
// g++ festival.cpp -o festival -Wall -Wextra -Wshadow -O2 -std=c++14 -fsanitize=address,undefined