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
#include<bits/stdc++.h>
using namespace std;
long long sumypre[2][10000010];
long long sumysuf[2][10000010];
int main()
{
	int n, m, i, j, p, szt;
	scanf("%d%d%d", &n, &m, &p);
	for(i=1;i<=m;i++)
	{
		for(j=i;j<=m;j++)
		{
			sumypre[1][j]+=1;
			sumysuf[1][i]+=1;
		}
	}
	
	for(i=1;i<=m;i++)
		sumypre[1][i] = (sumypre[1][i-1] + sumypre[1][i])%p;
	for(i=m;i>=1;i--)
		sumysuf[1][i] = (sumysuf[1][i+1] + sumysuf[1][i])%p;

	
	for(szt=2;szt<=n;szt++)
	{
		for(i=1;i<=m;i++)
		{
			sumypre[szt%2][i] = 0;
			sumysuf[szt%2][i] = 0;
		}
		for(i=1;i<=m;i++)
		{
			for(j=i;j<=m;j++)
			{
				//ile jest mozliwosci dla takiego
				long long w=0;
				w = sumypre[(szt-1)%2][m];
				w+=2*p;
				w-=sumypre[(szt-1)%2][i-1];
				w-=sumysuf[(szt-1)%2][j+1];
				w%=p;
				sumypre[szt%2][j]+=w;
				sumypre[szt%2][j]%=p;
				
				sumysuf[szt%2][i]+=w;
				sumysuf[szt%2][i]%=p;
				
			}
		}
		
		
		for(i=1;i<=m;i++)
			sumypre[szt%2][i] = (sumypre[szt%2][i-1] + sumypre[szt%2][i])%p;
		for(i=m;i>=1;i--)
			sumysuf[szt%2][i] = (sumysuf[szt%2][i+1] + sumysuf[szt%2][i])%p;
			
	}
	printf("%lld", sumypre[n%2][m]%p);
}