#include <bits/stdc++.h> #include <math.h> using namespace std; #define endl "\n" #define mp make_pair #define st first #define nd second #define pii pair<int, int> #define pb push_back #define _upgrade ios_base::sync_with_stdio(0), cout.setf(ios::fixed), cout.precision(10), cin.tie(0), cout.tie(0); #define REP(i, n) for (int i = 0; i < (n); ++i) #define FWD(i, a, b) for (int i = (a); i < (b); ++i) #define rep(i, n) for (int32_t i = 0; i < (n); ++i) #define fwd(i, a, b) for (int i = (a); i < (b); ++i) #define all(c) (c).begin(), (c).end() #define sz(X) (int)((X).size()) #define what(x) cerr << #x << " is " << x << endl; ostream &operator<<(ostream &out, string str) { for (char c : str) out << c; return out; } template <class L, class R> ostream &operator<<(ostream &out, pair<L, R> p) { return out << "(" << p.st << ", " << p.nd << ")"; } template <class T> auto operator<<(ostream &out, T a) -> decltype(a.begin(), out) { out << '{'; for (auto it = a.begin(); it != a.end(); it = next(it)) out << (it != a.begin() ? ", " : "") << *it; return out << '}'; } void dump() { cerr << "\n"; } template <class T, class... Ts> void dump(T a, Ts... x) { cerr << a << ", "; dump(x...); } #define DEBUG true #define debug(...) \ if (DEBUG) \ cerr << "[" #__VA_ARGS__ "]: ", dump(__VA_ARGS__) #define int long long int p; int inf; int add(int a, int b) { return (a + b + inf) % p; } int mult(int a, int b) { return (a * b) % p; } int32_t main() { int n, m; cin >> n >> m >> p; inf = p * p * 3; int res = 1; vector<int> X(m + 1); rep(i, n) { vector<int> cur(m + 1); int ac = 0; for (int j = 1; j <= m; j++) { int l = add(res, -X[m - j]); cur[j] = add(mult(j, l), -ac); cur[j] = add(cur[j], cur[j - 1]); ac = add(ac, X[j]); } X = cur; res = X.back(); // debug(X, res); } cout << res << endl; }
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 | #include <bits/stdc++.h> #include <math.h> using namespace std; #define endl "\n" #define mp make_pair #define st first #define nd second #define pii pair<int, int> #define pb push_back #define _upgrade ios_base::sync_with_stdio(0), cout.setf(ios::fixed), cout.precision(10), cin.tie(0), cout.tie(0); #define REP(i, n) for (int i = 0; i < (n); ++i) #define FWD(i, a, b) for (int i = (a); i < (b); ++i) #define rep(i, n) for (int32_t i = 0; i < (n); ++i) #define fwd(i, a, b) for (int i = (a); i < (b); ++i) #define all(c) (c).begin(), (c).end() #define sz(X) (int)((X).size()) #define what(x) cerr << #x << " is " << x << endl; ostream &operator<<(ostream &out, string str) { for (char c : str) out << c; return out; } template <class L, class R> ostream &operator<<(ostream &out, pair<L, R> p) { return out << "(" << p.st << ", " << p.nd << ")"; } template <class T> auto operator<<(ostream &out, T a) -> decltype(a.begin(), out) { out << '{'; for (auto it = a.begin(); it != a.end(); it = next(it)) out << (it != a.begin() ? ", " : "") << *it; return out << '}'; } void dump() { cerr << "\n"; } template <class T, class... Ts> void dump(T a, Ts... x) { cerr << a << ", "; dump(x...); } #define DEBUG true #define debug(...) \ if (DEBUG) \ cerr << "[" #__VA_ARGS__ "]: ", dump(__VA_ARGS__) #define int long long int p; int inf; int add(int a, int b) { return (a + b + inf) % p; } int mult(int a, int b) { return (a * b) % p; } int32_t main() { int n, m; cin >> n >> m >> p; inf = p * p * 3; int res = 1; vector<int> X(m + 1); rep(i, n) { vector<int> cur(m + 1); int ac = 0; for (int j = 1; j <= m; j++) { int l = add(res, -X[m - j]); cur[j] = add(mult(j, l), -ac); cur[j] = add(cur[j], cur[j - 1]); ac = add(ac, X[j]); } X = cur; res = X.back(); // debug(X, res); } cout << res << endl; } |