#include<cstdio>
using namespace std;
const int o=40010;
int n,p,fac[o],inv[o],f[o][2],ans=1;
inline int fix(int x){return x+(x>>31&p);}
int main(){
	scanf("%d%d",&n,&p);inv[1]=1;
	for(int i=2;i<=n*2;++i) inv[i]=p-p/i*1ll*inv[p%i]%p;
	for(int i=fac[0]=inv[0]=1;i<=n*2;++i) fac[i]=fac[i-1]*1ll*i%p,inv[i]=inv[i-1]*1ll*inv[i]%p;
	f[1][1]=f[2][0]=1;
	for(int i=1;i<=n;++i){
		for(int j=0,v=f[i][0],t;j<n-i&&v;++j)
			t=v*(j+1ll)%p,f[i+j+2][0]=(f[i+j+2][0]+t*(j+2ll))%p,f[i+j+1][1]=fix(f[i+j+1][1]+t-p),v=v*(2*i-2ll+j)%p;
		for(int j=0,v=f[i][1];j<n-i&&v;++j)
			f[i+j+2][0]=(f[i+j+2][0]+v*(j+1ll))%p,f[i+j+1][1]=fix(f[i+j+1][1]+v-p),v=v*(2*i-1ll+j)%p;
	}
	for(int i=1;i<=n;++i) ans=ans*(2*i-1ll)%p;
	for(int i=2;i<=n;++i) ans=(ans+(p-f[i][0])*(n-i+1ll)%p*fac[n+i-3]%p*inv[i*2-3])%p;
	for(int i=1;i<=n;++i) ans=(ans+(p-f[i][1])*1ll*fac[n+i-2]%p*inv[i*2-2])%p;
	printf("%d",ans);
	return 0;
}