#include <bits/stdc++.h> using namespace std; #define FOR(i, a, b) for(int i = (a); i <= (b); ++i) #define REP(i, n) for(int i = 0; i < (n); ++i) #define ff first #define ss second using ll = long long; using pii = pair<int, int>; template<class T> int size(T && a) { return (int)(a.size()); } ostream& operator << (ostream &os, string str) { for(char c : str) os << c; return os; } template <class A, class B> ostream& operator << (ostream &os, const pair<A, B> &p) { return os << '(' << p.ff << "," << p.ss << ')'; } template <class T> auto operator << (ostream &os, T &&x) -> decltype(x.begin(), os) { os << '{'; for(auto it = x.begin(); it != x.end(); ++it) os << *it << (it == prev(x.end()) ? "" : " "); return os << '}'; } template <class T> ostream& operator << (ostream &os, vector<vector<T>> vec) { for(auto x : vec) os << "\n " << x; return os; } void dump() {} template <class T, class... Args> void dump(T &&x, Args... args) { cerr << x << "; "; dump(args...); } #ifdef DEBUG # define debug(x...) cerr << "[" #x "]: ", dump(x), cerr << '\n' #else # define debug(...) 0 #endif //-------------------------------------------------- int MOD; int add(int a, int b) { if(a+b >= MOD) return a+b-MOD; else return a+b; } int main() { ios_base::sync_with_stdio(0); cin.tie(0); int n, m; cin >> n >> m >> MOD; swap(n, m); vector<vector<int>> cnt(n, vector<int>(n)); REP(i, n) FOR(j, i, n-1) cnt[i][j] = 1; debug(cnt); REP(rep, m-1) { vector<vector<int>> new_cnt(n, vector<int>(n)); REP(i, n) FOR(j, i, n-1) { REP(i1, n) FOR(j1, i1, n-1) { if(j < i1 or i > j1) continue; new_cnt[i][j] = add(new_cnt[i][j], cnt[i1][j1]); } } swap(cnt, new_cnt); debug(cnt); } int ans = 0; REP(i, n) FOR(j, i, n-1) ans = add(ans, cnt[i][j]); cout << ans << '\n'; }
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 | #include <bits/stdc++.h> using namespace std; #define FOR(i, a, b) for(int i = (a); i <= (b); ++i) #define REP(i, n) for(int i = 0; i < (n); ++i) #define ff first #define ss second using ll = long long; using pii = pair<int, int>; template<class T> int size(T && a) { return (int)(a.size()); } ostream& operator << (ostream &os, string str) { for(char c : str) os << c; return os; } template <class A, class B> ostream& operator << (ostream &os, const pair<A, B> &p) { return os << '(' << p.ff << "," << p.ss << ')'; } template <class T> auto operator << (ostream &os, T &&x) -> decltype(x.begin(), os) { os << '{'; for(auto it = x.begin(); it != x.end(); ++it) os << *it << (it == prev(x.end()) ? "" : " "); return os << '}'; } template <class T> ostream& operator << (ostream &os, vector<vector<T>> vec) { for(auto x : vec) os << "\n " << x; return os; } void dump() {} template <class T, class... Args> void dump(T &&x, Args... args) { cerr << x << "; "; dump(args...); } #ifdef DEBUG # define debug(x...) cerr << "[" #x "]: ", dump(x), cerr << '\n' #else # define debug(...) 0 #endif //-------------------------------------------------- int MOD; int add(int a, int b) { if(a+b >= MOD) return a+b-MOD; else return a+b; } int main() { ios_base::sync_with_stdio(0); cin.tie(0); int n, m; cin >> n >> m >> MOD; swap(n, m); vector<vector<int>> cnt(n, vector<int>(n)); REP(i, n) FOR(j, i, n-1) cnt[i][j] = 1; debug(cnt); REP(rep, m-1) { vector<vector<int>> new_cnt(n, vector<int>(n)); REP(i, n) FOR(j, i, n-1) { REP(i1, n) FOR(j1, i1, n-1) { if(j < i1 or i > j1) continue; new_cnt[i][j] = add(new_cnt[i][j], cnt[i1][j1]); } } swap(cnt, new_cnt); debug(cnt); } int ans = 0; REP(i, n) FOR(j, i, n-1) ans = add(ans, cnt[i][j]); cout << ans << '\n'; } |