結果
問題 | No.1667 Forest |
ユーザー | siro53 |
提出日時 | 2021-09-03 23:39:17 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 744 ms / 3,000 ms |
コード長 | 3,029 bytes |
コンパイル時間 | 2,003 ms |
コンパイル使用メモリ | 209,964 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-05-09 07:49:12 |
合計ジャッジ時間 | 7,261 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 744 ms
5,248 KB |
testcase_01 | AC | 732 ms
5,376 KB |
testcase_02 | AC | 706 ms
5,376 KB |
testcase_03 | AC | 10 ms
5,376 KB |
testcase_04 | AC | 703 ms
5,376 KB |
testcase_05 | AC | 409 ms
5,376 KB |
testcase_06 | AC | 238 ms
5,376 KB |
testcase_07 | AC | 142 ms
5,376 KB |
testcase_08 | AC | 68 ms
5,376 KB |
testcase_09 | AC | 42 ms
5,376 KB |
testcase_10 | AC | 23 ms
5,376 KB |
testcase_11 | AC | 7 ms
5,376 KB |
testcase_12 | AC | 4 ms
5,376 KB |
testcase_13 | AC | 4 ms
5,376 KB |
testcase_14 | AC | 5 ms
5,376 KB |
testcase_15 | AC | 4 ms
5,376 KB |
testcase_16 | AC | 4 ms
5,376 KB |
testcase_17 | AC | 4 ms
5,376 KB |
ソースコード
#line 1 "main.cpp" #pragma region Macros #include <bits/stdc++.h> using namespace std; template <class T> inline bool chmax(T &a, T b) { if(a < b) { a = b; return 1; } return 0; } template <class T> inline bool chmin(T &a, T b) { if(a > b) { a = b; return 1; } return 0; } #ifdef DEBUG template <class T, class U> ostream &operator<<(ostream &os, const pair<T, U> &p) { os << '(' << p.first << ',' << p.second << ')'; return os; } template <class T> ostream &operator<<(ostream &os, const vector<T> &v) { os << '{'; for(int i = 0; i < (int)v.size(); i++) { if(i) { os << ','; } os << v[i]; } os << '}'; return os; } void debugg() { cerr << endl; } template <class T, class... Args> void debugg(const T &x, const Args &... args) { cerr << " " << x; debugg(args...); } #define debug(...) \ cerr << __LINE__ << " [" << #__VA_ARGS__ << "]: ", debugg(__VA_ARGS__) #define dump(x) cerr << __LINE__ << " " << #x << " = " << (x) << endl #else #define debug(...) (void(0)) #define dump(x) (void(0)) #endif struct Setup { Setup() { cin.tie(0); ios::sync_with_stdio(false); cout << fixed << setprecision(15); } } __Setup; using ll = long long; #define ALL(v) (v).begin(), (v).end() #define RALL(v) (v).rbegin(), (v).rend() #define FOR(i, a, b) for(int i = (a); i < int(b); i++) #define REP(i, n) FOR(i, 0, n) const int INF = 1 << 30; const ll LLINF = 1LL << 60; constexpr int MOD = 1000000007; const int dx[4] = {1, 0, -1, 0}; const int dy[4] = {0, 1, 0, -1}; void Case(int i) { cout << "Case #" << i << ": "; } int popcount(int x) { return __builtin_popcount(x); } ll popcount(ll x) { return __builtin_popcountll(x); } #pragma endregion Macros #include <atcoder/modint> using mint = atcoder::modint; #line 1 "/home/siro53/kyo-pro/compro_library/math/combination.hpp" struct Combination { int MAX; vector<mint> fac, finv; Combination(int _MAX) : fac(_MAX), finv(_MAX), MAX(_MAX) { fac[0] = fac[1] = 1; for(int i = 2; i < MAX; i++) { fac[i] = fac[i - 1] * i; } finv[MAX - 1] = fac[MAX - 1].inv(); for(int i = MAX - 1; i >= 1; i--) { finv[i - 1] = finv[i] * i; } } mint COM(int n, int k) { if(n < k || n < 0 || k < 0) { return 0; } return fac[n] * finv[k] * finv[n - k]; } }; #line 75 "main.cpp" mint cayley(int n) { if(n <= 1) return 1; return mint(n).pow(n-2); } int main() { int N, mod; cin >> N >> mod; mint::set_mod(mod); Combination comb(100000); vector dp(N+1, vector<mint>(N+1, 0)); dp[0][0] = 1; REP(i, N) { REP(j, N) { FOR(k, 1, N+1) { if(i + k > N) break; dp[i+k][j+1] += dp[i][j] * comb.COM(N-i-1, k-1) * cayley(k); } } } REP(M, N) { int m = N - M; mint ans = dp[N][m]; cout << ans.val() << endl; } }