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
56
57
58
59
60
61
// Jakub Rożek
// Malowanie płotu [B]
// Potyczki 2020
// O( n*m^3)

#include <bits/stdc++.h>
using namespace std;

const int M=5003;

int n,m,kk,kkk;
long long mod,odp;
long long t[2][M][M];

int main()
{
	ios_base::sync_with_stdio(0);
	cin.tie(0);
	
	cin>>n>>m>>mod;
	if(m>=M) return 0;
	
	for(int i=1; i<=m; ++i)
	{
		t[0][i][i-1]=0;
		for(int j=i; j<=m; ++j)
		{
			t[0][i][j]=t[0][i][j-1]+1;
		}
	}
	
	for(int k=1; k<n; ++k)
	{
		kk=k%2;
		kkk=(kk+1)%2;
		
		for(int i=1; i<=m; ++i)
		{
			t[kk][i][i-1]=0;
			for(int j=i; j<=m; ++j)
			{
				t[kk][i][j]=t[kk][i][j-1];
				
				for(int ii=1; ii<=j; ++ii)
				{
					t[kk][i][j]+=t[kkk][ii][m]-t[kkk][ii][max(i,ii)-1]+mod;
					t[kk][i][j]%=mod;
				}
			}
		}
	}
	
	for(int i=1; i<=m; ++i)
	{
		odp+=t[kk][i][m]-t[kk][i][i-1]+mod;
		odp%=mod;
	}
	
	cout<<odp;
	return(0);
}