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
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
#include<bits/stdc++.h>

using namespace std;

long long matrix1[2000][2000];
long long matrix2[2000][2000];

int main()
{
    int n,m,p;
    cin>>n>>m>>p;
    for(int x=0;x<m;x++)
        for(int y=0;y<=x;y++)
            matrix1[x][y]=1;
    bool even = 1;
    for(int i=0;i<n;i++)
    {
        if(even)
        {
            for(int x=0;x<m;x++)
            {
                for(int y=0;y<=x;y++)
                {
                    matrix2[x][y]=0;
                    for(int x2=y;x2<m;x2++)
                    {
                        for(int y2=0;y2<=x;y2++)
                        {
                            matrix2[x][y]=(matrix2[x][y] + matrix1[x2][y2]) % p;
                        }
                    }
                }
            }
            // for(int i=0;i<m;i++)
            // {
            //     for(int j=0;j<m;j++)
            //     {
            //         cerr<<setw(10)<<matrix2[j][i];
            //     }
            //     cerr<<endl;
            // }
        }
        else
        {
            for(int x=0;x<m;x++)
            {
                for(int y=0;y<=x;y++)
                {
                    matrix1[x][y]=0;
                    for(int x2=y;x2<m;x2++)
                    {
                        for(int y2=0;y2<=x;y2++)
                        {
                            matrix1[x][y]=(matrix1[x][y]+matrix2[x2][y2]) % p;
                        }
                    }
                }
            }
            // for(int i=0;i<m;i++)
            // {
            //     for(int j=0;j<m;j++)
            //     {
            //         cerr<<setw(10)<<matrix1[j][i];
            //     }
            //     cerr<<endl;
            // }
        }
        //cerr<<endl;
        even = !even;
    }
    if(n%2)
    {
        cout<<matrix2[m-1][0]<<endl;
        return 0;
    }
    else
    {
        cout<<matrix1[m-1][0]<<endl;
        return 0;
    }
}