結果
問題 | No.1667 Forest |
ユーザー | Kude |
提出日時 | 2021-09-03 22:28:44 |
言語 | C++17 (gcc 13.2.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 69 ms / 3,000 ms |
コード長 | 1,901 bytes |
コンパイル時間 | 4,222 ms |
コンパイル使用メモリ | 261,376 KB |
実行使用メモリ | 11,924 KB |
最終ジャッジ日時 | 2023-08-21 23:10:38 |
合計ジャッジ時間 | 6,009 ms |
ジャッジサーバーID (参考情報) |
judge11 / judge13 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 69 ms
11,704 KB |
testcase_01 | AC | 68 ms
11,656 KB |
testcase_02 | AC | 67 ms
11,660 KB |
testcase_03 | AC | 28 ms
11,600 KB |
testcase_04 | AC | 67 ms
11,588 KB |
testcase_05 | AC | 51 ms
11,872 KB |
testcase_06 | AC | 41 ms
11,924 KB |
testcase_07 | AC | 36 ms
11,688 KB |
testcase_08 | AC | 31 ms
11,584 KB |
testcase_09 | AC | 29 ms
11,620 KB |
testcase_10 | AC | 28 ms
11,600 KB |
testcase_11 | AC | 27 ms
11,676 KB |
testcase_12 | AC | 27 ms
11,656 KB |
testcase_13 | AC | 28 ms
11,644 KB |
testcase_14 | AC | 27 ms
11,620 KB |
testcase_15 | AC | 27 ms
11,596 KB |
testcase_16 | AC | 27 ms
11,608 KB |
testcase_17 | AC | 27 ms
11,872 KB |
ソースコード
#include<bits/stdc++.h> #include<atcoder/all> using namespace std; using namespace atcoder; #define rep(i,n)for (int i = 0; i < int(n); ++i) #define rrep(i,n)for (int i = int(n)-1; i >= 0; --i) #define all(x) (x).begin(), (x).end() #define rall(x) (x).rbegin(), (x).rend() template<class T> void chmax(T& a, const T& b) {a = max(a, b);} template<class T> void chmin(T& a, const T& b) {a = min(a, b);} using ll = long long; using P = pair<int,int>; using VI = vector<int>; using VVI = vector<VI>; using VL = vector<ll>; using VVL = vector<VL>; using mint = modint; constexpr int FACT_SIZE = 1000000; mint Fact[FACT_SIZE + 1]; mint iFact[FACT_SIZE + 1]; const auto fact_init = [] { Fact[0] = mint::raw(1); for(int i = 1; i <= FACT_SIZE; ++i) { Fact[i] = Fact[i-1] * i; } iFact[FACT_SIZE] = Fact[FACT_SIZE].inv(); for(int i = FACT_SIZE; i; --i) { iFact[i-1] = iFact[i] * i; } return false; }; mint comb(int n, int k) { if (k == 0) return mint::raw(1); assert(n >= 0 && k >= 0); if (k > n) return mint::raw(0); return Fact[n] * iFact[n - k] * iFact[k]; } mint icomb(int n, int k) { return iFact[n] * Fact[n - k] * Fact[k]; } mint fact(int n) {return Fact[n];} mint perm(int n, int k) { assert(0 <= n); return Fact[n] * iFact[n - k]; } mint dp[310][310]; int main() { ios::sync_with_stdio(false); cin.tie(0); int n, mod; cin >> n >> mod; mint::set_mod(mod); fact_init(); dp[0][0] = 1; for(int sz = 1; sz <= n; sz++) for(int c = 1; c <= sz; c++) { mint cnt; for(int i = 1; i < sz && c - 1 <= sz - i; i++) { cnt += dp[i][1] * dp[sz - i][c - 1] * comb(sz - 1, i - 1); } if (c == 1) { if (sz == 1) cnt++; else cnt += mint(sz).pow(sz - 2); } dp[sz][c] = cnt; } rep(i, n) cout << dp[n][n - i].val() << '\n'; }