#include<bits/stdc++.h>
using namespace std;
int mod, n;
typedef long long ll;
int add(int a,int b){a+=b;return a>=mod?a-mod:a;}
int sub(int a,int b){a-=b;return a<0?a+mod:a;}
int mul(int a,int b){return (ll)a*b%mod;}
int qpow(int a,int b){int ret=1;for(;b;b>>=1,a=mul(a,a))if(b&1)ret=mul(ret,a);return ret;}
/* math */
const int N = 5010;
int f[2][N];
int g[2][N];
int h1[2][N];
int h2[2][N];

inline void Do0(int p,int i){
	for(int j=0;j<=n;j++)f[i][j] = g[i][j] = h1[i][j] = h2[i][j] = 0;
	for(int j=0;j<=n;j++){
		f[i][j] = add(f[i][j], f[i^1][j]);
		// g[i][j] = add(g[i][j], g[i^1][j]);
		if(j>=p){
			f[i][j] = add(f[i][j], f[i][j-p]);
			// g[i][j] = add(g[i][j], g[i][j-p]);// + 1
			// g[i][j] = add(g[i][j], f[i][j-p]);
		}
	}
}

inline void Do1(int p,int i){
	for(int j=0;j<=n;j++)f[i][j] = g[i][j] = h1[i][j] = h2[i][j] = 0;
	for(int j=0;j<=n;j++){
		f[i][j] = add(f[i][j], f[i^1][j]);
		// g[i][j] = add(g[i][j], g[i^1][j]);
		if(j>=p){
			f[i][j] = add(f[i][j], f[i][j-p]);
			// g[i][j] = add(g[i][j], g[i][j-p]);// + 1
			// g[i][j] = add(g[i][j], f[i][j-p]);
		}
	}
}
inline void Do2(int p,int i){
	for(int j=0;j<=n;j++)f[i][j] = g[i][j] = h1[i][j] = h2[i][j] = 0;
	for(int j=0;j<=n;j++){
		f[i][j] = add(f[i][j], f[i^1][j]);
		// g[i][j] = add(g[i][j], g[i+1][j]);
		if(j>=p){
			f[i][j] = add(f[i][j], f[i][j-p]);
			// g[i][j] = add(g[i][j], g[i][j-p]);// + 1
			// g[i][j] = add(g[i][j], f[i][j-p]);
		}
	}
	for(int j=0;j<=n;j++){
		h1[i][j] = add(h1[i][j], 1ll*f[i^1][j]*(j/p-1)%mod);
		// h2[i][j] = add(h2[i][j], 1ll*g[i+1][j]*(j/p)%mod);
		if(j>=p){
			h1[i][j] = add(h1[i][j], h1[i][j-p]);
			// h2[i][j] = add(h2[i][j], h2[i][j-p]);// + 1
			// h2[i][j] = add(h2[i][j], h1[i][j-p]);
		}
	}
	for(int j=0;j<=n;j++){
		f[i][j] = sub(mul(f[i][j], j/p),h1[i][j]);
		// g[i][j] = sub(mul(g[i][j], j/p),h2[i][j]);
	}
}

int main()
{
	cin >> n >> mod;
	int _d = (n-1)/2;
	bool flg=0;
	if(n%2==0){ flg=1; }
	// cout << _d << endl;
	int cur = 0;
	f[cur][0] = 1;
	// for(int st=1;st<=n;st++){
	// 	f[_d+1][n-st+1]=1;//max=st-1
	// 	g[_d+1][n-st+1]=st-1;
	// 	// cout << st << ":" << n-st+1 << endl;
	// }
	for(int i=_d;i;i--){
		int p = i;
		Do0(p,cur^=1);
		if(i==_d && flg)Do2(p+1,cur^=1);
		else Do1(p+1,cur^=1);
		// if(i==_d && flg==0){
		// 	for(int j=0;j<=n;j++){
		// 		f[i][j] = add(f[i][j], f[i+1][j]);
		// 		g[i][j] = add(g[i][j], g[i+1][j]);
		// 		if(j>=p){
		// 			f[i][j] = add(f[i][j], f[i][j-p]);
		// 			g[i][j] = add(g[i][j], g[i][j-p]);// + 1
		// 			g[i][j] = add(g[i][j], f[i][j-p]);
		// 		}
		// 	}
		// }else{
		// 	for(int j=0;j<=n;j++){
		// 		f[i][j] = add(f[i][j], f[i+1][j]);
		// 		g[i][j] = add(g[i][j], g[i+1][j]);
		// 		if(j>=p){
		// 			f[i][j] = add(f[i][j], f[i][j-p]);
		// 			g[i][j] = add(g[i][j], g[i][j-p]);// + 1
		// 			g[i][j] = add(g[i][j], f[i][j-p]);
		// 		}
		// 	}
		// 	for(int j=0;j<=n;j++){
		// 		h1[i][j] = add(h1[i][j], 1ll*f[i+1][j]*(j/p)%mod);
		// 		h2[i][j] = add(h2[i][j], 1ll*g[i+1][j]*(j/p)%mod);
		// 		if(j>=p){
		// 			h1[i][j] = add(h1[i][j], h1[i][j-p]);
		// 			h2[i][j] = add(h2[i][j], h2[i][j-p]);// + 1
		// 			h2[i][j] = add(h2[i][j], h1[i][j-p]);
		// 		}
		// 	}
		// 	for(int j=0;j<=n;j++){
		// 		f[i][j] = sub(mul(f[i][j], j/p),h1[i][j]);
		// 		g[i][j] = sub(mul(g[i][j], j/p),h2[i][j]);
		// 	}
		// }
	}
	int ans = 0;
	if(n==2){cout << 3 << endl;return 0;}
	for(int i=0;i<=n;i++){
		// cout << f[cur][i] << " " ;
		// cout << n-i << endl;
		// cout << f[1][i] << " " << g[1][i] << endl;
		//p >= i;
		ans = add(ans, mul(f[cur][i], n-i));
		// ans = sub(ans, g[cur][i]);
	}
	cout << ans << endl;
}

Measured memory_bytes: 384