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
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
// mal-malowanie-plotu.cpp : Ten plik zawiera funkcję „main”. W nim rozpoczyna się i kończy wykonywanie programu.
//

#include <iostream>

constexpr int MAXM = 10000007;

long long beg[2][MAXM];
long long end[2][MAXM];

int main()
{
	int n, q = 0;
	long long p, md, m, i, lsm, j, k, tmp, p2;
	std::cin >> n >> m >> md;
	n--;
	for (i = 1; i <= m; i++)
	{
		end[1][i] = i;
		beg[1][i] = m - i + 1;
	}
	lsm = (m * (m + 1) / 2)%md;
	while (n--)
	{
		p = 0;
		p2 = 0;
		/*
		if (q)
		{
			for (i = 1; i <= m; i++)
			{
				p += dp[!q][i];
				p %= md;
				dp[q][i] = (dp[q][i - 1] + p) % md;
			}
		}
		else
		{
			for (i = m; i > 0; i--)
			{
				p += dp[!q][i];
				p %= md;
				dp[q][i] = (dp[q][i + 1] + p) % md;
			}
		}*/
		for (i = 1; i <= m; i++)
		{
			end[q][i] = (i * lsm)%md;
			end[q][i] = ((end[q][i]-p2)%md+md)%md;
			p += end[!q][i];
			p %= md;
			p2 += p;
			p2 %= md;
		}
		p = p2 = 0;
		for (i = m; i > 0; i--)
		{
			end[q][i] = ((end[q][i] - p * i) % md + md) % md;
			end[q][i] %= md;
			p += beg[!q][i];
			p %= md;
		}
		/*
		for (i = 1; i <= m; i++)
		{
			beg[q][i] = end[q][i] = 0;
		}
		
		for (i = 1; i <= m; i++)
		{
			for (j = 1; j <= i; j++)
			{
				// (j, i)
				tmp = lsm;
				for (k = 1; k < j; k++)
				{
					tmp -= end[!q][k];
				}
				for (k = i+1; k <= m; k++)
				{
					tmp -= beg[!q][k];
				}
				end[q][i] += tmp;
				beg[q][j] += tmp;
			}
		}*/
		lsm = 0;
		for (size_t i = 1; i <= m; i++)
		{
			beg[q][i] = end[q][m - i + 1];
			lsm += beg[q][i];
			//lsm %= md;
		//	std::cout << beg[q][i] << ' ';
		}
	/*	std::cout << '\n';
		for (size_t i = 1; i <= m; i++)
		{
			std::cout << end[q][i] << ' ';
		}
		std::cout << '\n';*/
		q = !q;
	}
	p = 0;
	q = !q;/*
	for (i = 1; i <= m; i++)
	{
		p += dp[q][i];
		p %= md;
	}
	std::cout << p << '\n';*/
	std::cout << lsm%md << '\n';
}

// Uruchomienie programu: Ctrl + F5 lub menu Debugowanie > Uruchom bez debugowania
// Debugowanie programu: F5 lub menu Debugowanie > Rozpocznij debugowanie

// Porady dotyczące rozpoczynania pracy:
//   1. Użyj okna Eksploratora rozwiązań, aby dodać pliki i zarządzać nimi
//   2. Użyj okna programu Team Explorer, aby nawiązać połączenie z kontrolą źródła
//   3. Użyj okna Dane wyjściowe, aby sprawdzić dane wyjściowe kompilacji i inne komunikaty
//   4. Użyj okna Lista błędów, aby zobaczyć błędy
//   5. Wybierz pozycję Projekt > Dodaj nowy element, aby utworzyć nowe pliki kodu, lub wybierz pozycję Projekt > Dodaj istniejący element, aby dodać istniejące pliku kodu do projektu
//   6. Aby w przyszłości ponownie otworzyć ten projekt, przejdź do pozycji Plik > Otwórz > Projekt i wybierz plik sln