#include <bits/stdc++.h>
using namespace std;
#define For(Ti,Ta,Tb) for(auto Ti=(Ta);Ti<=(Tb);++Ti)
#define Dec(Ti,Ta,Tb) for(auto Ti=(Ta);Ti>=(Tb);--Ti)
#define debug(...) fprintf(stderr,__VA_ARGS__)
#define range(Tx) begin(Tx),end(Tx)
using ll=long long;
using i128=__int128;
const int N=2e4+5;
int n,mod;
ll fac[N*2],inv[N*2],ifac[N*2];
void init_fac(){
	fac[0]=1;
	For(i,1,N*2-1) fac[i]=fac[i-1]*i%mod;
	inv[1]=1;
	For(i,2,N*2-1) inv[i]=(mod-mod/i)*inv[mod%i]%mod;
	ifac[0]=ifac[1]=1;
	For(i,2,N*2-1) ifac[i]=ifac[i-1]*inv[i]%mod;
}
void reduce(ll& x){
	x=(x>=mod?x-mod:x);
}
ll f[2][N][3];
i128 g[N*2][3];
// 0：没有劣势
// 1：没有劣势，有左端点
// 2：有劣势，有左端点
int main(){
	cin.tie(nullptr)->sync_with_stdio(false);
	cin>>n>>mod;
	init_fac();
	auto ff=[&](int n,int m){
		return n<m?0:fac[n]*ifac[n-m];
	};
	f[0][0][0]=1;
	For(i,0,n*2-1){
		memset(f[~i&1],0,sizeof(ll)*min(n,i+1)*3);
		For(k,0,min({n,i,n*2-i})){
			// 左括号
			reduce(f[~i&1][k+1][1]+=f[i&1][k][0]);
			reduce(f[~i&1][k+1][1]+=f[i&1][k][1]);
			reduce(f[~i&1][k+1][2]+=f[i&1][k][2]);
			// 右括号
			reduce(f[~i&1][k][0]+=f[i&1][k][2]);
			if(k){
				if(k>1){
					g[i+k-1][2]+=i128(k-1)*f[i&1][k][1]*ff(n*2-i-1,k-2);
				}
				g[i+k][0]+=i128(f[i&1][k][1])*ff(n*2-i-1,k-1);
			}
		}
		reduce(f[~i&1][0][0]+=g[i+1][0]%mod);
		reduce(f[~i&1][0][2]+=g[i+1][2]%mod);
	}
	ll all=1;
	for(int i=n*2-1;i>=1;i-=2){
		(all*=i)%=mod;
	}
	cout<<(all-f[0][0][0]+mod)%mod<<'\n';
	return 0;
}