1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
#include <bits/stdc++.h>
using namespace std;
void imax(int &a, int b){
	a=max(a, b);
}
void imin(int &a, int b){
	a=min(a, b);
}
void lmax(long long &a, long long b){
	a=max(a, b);
}
void lmin(long long &a, long long b){
	a=min(a, b);
}
/*
	WARNING: I'm using strange bracket style!
*/
vector <vector <long long> > dp;
vector <long long> x;
long long MOD;
int main()
	{
	ios::sync_with_stdio(0);
	cin.tie(0);
	cout.tie(0);
	int n, m;
	cin>>n>>m>>MOD;
	for (int i=0; i<=m; i++)
		x.push_back(0);
	for (int i=0; i<=n; i++)
		dp.push_back(x);
	for (int i=1; i<=m; i++)
		dp[1][i]=i;
	for (int i=2; i<=n; i++)
		{
		long long S=0;
		for (int j=1; j<=m; j++)
			(S+=dp[i-1][j])%=MOD;
		for (int j=1; j<=m; j++)
			dp[i][j]=(j*S)%MOD;
		for (int j=1; j<=m; j++)
			x[j]=(dp[i-1][j]+x[j-1])%MOD;
		for (int j=1; j<=m; j++)
			dp[i][j]-=(x[m-j]*j)%MOD;
		for (int j=1; j<=m; j++)
			(x[j]+=x[j-1])%=MOD;
		for (int j=1; j<=m; j++)
			dp[i][j]=((dp[i][j]-x[j-1])%MOD+MOD)%MOD;
		}
	long long RES=0;
	for (int j=1; j<=m; j++)
		RES+=dp[n][j];
	cout<<RES%MOD<<"\n";
	return 0;
	}