#include<bits/stdc++.h>
using namespace std;
const int N=4e5+5;

int n,P,fc[N],ifc[N];
int o,w,f[2][N][3],ad[N][3],ans,all;

//calc
inline void inc(int &x,int y) { x+=y-P; x+=(x>>31)&P; }
inline void dec(int &x,int y) { x-=y; x+=(x>>31)&P; }
inline int pls(int x,int y) { return x+=y-P, x+=(x>>31)&P, x; }
inline int power(int a,int b) { int res=1; for(;b;b>>=1,a=1ll*a*a%P) if(b&1) res=1ll*res*a%P; return res; }

int main() {

	cin>>n>>P;
	
	fc[0]=ifc[0]=all=1; 
	for(int i=1;i<N;i++) fc[i]=1ll*i*fc[i-1]%P;
	ifc[N-1]=power(fc[N-1],P-2);
	for(int i=N-2;i>=1;i--) ifc[i]=1ll*ifc[i+1]*(i+1)%P;
	
	int *F, m=2*n; o=0; w=1; f[o][0][0]=1; 
	for(int i=0;i<=m;i++,swap(o,w)) {
		for(int j=0;j<3;j++) inc(f[o][0][j],ad[i][j]);
		if(i==m) break;
		for(int k=0;k<=min(n,i);k++) {
			F = f[o][k];
			/*
			cout<<endl;
			cout<<i<<" "<<k<<endl;
			cout<<0<<" "<<F[0]<<endl;
			cout<<1<<" "<<F[1]<<endl;
			cout<<2<<" "<<F[2]<<endl;
			*/
			inc(f[w][k+1][1],pls(F[0],F[1]));
			inc(f[w][k+1][2],F[2]);
			inc(f[w][k][0],F[2]);
			if(k && i+k<=m) inc(ad[i+k][0],1ll*F[1]*fc[m-i-1]%P*ifc[m-i-k]%P);
			if(k>1 && i+k<=m) inc(ad[i+k-1][2],1ll*F[1]*(k-1)%P*fc[m-i-1]%P*ifc[m-i-k+1]%P);
			F[0]=F[1]=F[2]=0;
			//f[o][k][0]=f[o][k][1]=f[o][k][2]=0;
			//if(i==3) cout<<endl<<k<<" "<<f[w][0][0]<<endl<<endl;
		}
	}
	
	for(int i=1;i<=m;i+=2) all=1ll*all*i%P;
	dec(ans=all,f[o][0][0]);
	cout<<ans<<endl;
	
	return 0;
}

/*
2 1000000007
3 100000007
4 100000007
15 999999937 
*/