結果
問題 | No.1667 Forest |
ユーザー | milanis48663220 |
提出日時 | 2021-09-03 22:54:55 |
言語 | C++17 (gcc 13.2.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 2,346 ms / 3,000 ms |
コード長 | 2,291 bytes |
コンパイル時間 | 1,259 ms |
コンパイル使用メモリ | 122,204 KB |
実行使用メモリ | 116,752 KB |
最終ジャッジ日時 | 2023-08-22 00:10:27 |
合計ジャッジ時間 | 17,768 ms |
ジャッジサーバーID (参考情報) |
judge11 / judge12 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2,346 ms
116,556 KB |
testcase_01 | AC | 2,321 ms
116,528 KB |
testcase_02 | AC | 2,274 ms
116,588 KB |
testcase_03 | AC | 54 ms
116,752 KB |
testcase_04 | AC | 2,248 ms
116,512 KB |
testcase_05 | AC | 1,327 ms
116,544 KB |
testcase_06 | AC | 806 ms
116,516 KB |
testcase_07 | AC | 498 ms
116,524 KB |
testcase_08 | AC | 252 ms
116,616 KB |
testcase_09 | AC | 160 ms
116,680 KB |
testcase_10 | AC | 96 ms
116,524 KB |
testcase_11 | AC | 44 ms
116,548 KB |
testcase_12 | AC | 39 ms
116,544 KB |
testcase_13 | AC | 39 ms
116,608 KB |
testcase_14 | AC | 39 ms
116,612 KB |
testcase_15 | AC | 39 ms
116,512 KB |
testcase_16 | AC | 38 ms
116,584 KB |
testcase_17 | AC | 38 ms
116,744 KB |
ソースコード
#include <iostream> #include <algorithm> #include <iomanip> #include <vector> #include <queue> #include <set> #include <map> #include <tuple> #include <cmath> #include <numeric> #include <functional> #include <cassert> #include <atcoder/modint> #define debug_value(x) cerr << "line" << __LINE__ << ":<" << __func__ << ">:" << #x << "=" << x << endl; #define debug(x) cerr << "line" << __LINE__ << ":<" << __func__ << ">:" << x << endl; 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; } using namespace std; typedef long long ll; using mint = atcoder::modint; ll MOD = 1; mint pow(mint a, ll n) { assert(n >= 0); mint ans = 1; mint tmp = a; for (int i = 0; i <= 60; i++) { ll m = (ll)1 << i; if (m & n) { ans *= tmp; } tmp *= tmp; } return ans; } #define N_MAX 200002 mint inv[N_MAX],fac[N_MAX],finv[N_MAX]; void init(){ fac[0]=1;fac[1]=1; finv[0]=1;finv[1]=1; inv[1]=1; for(int i=2;i<N_MAX;i++){ inv[i]=mint(MOD)-inv[MOD%i]*(MOD/i); fac[i]=fac[i-1]*(ll) i; finv[i]=finv[i-1]*inv[i]; } } mint comb(ll n, ll r){ if(n < r){ return mint(0); }else{ return fac[n]*finv[r]*finv[n-r]; } } mint dp[305][305][305]; int main(){ ios::sync_with_stdio(false); cin.tie(0); cout << setprecision(10) << fixed; int n; cin >> n >> MOD; atcoder::modint::set_mod(MOD); init(); vector<mint> v(n+1); v[0] = 0; v[1] = 1; for(int i = 2; i <= n; i++) v[i] = pow(mint(i), i-2); dp[0][0][0] = 1; for(int i = 1; i <= n; i++){ for(int j = 0; j <= n; j++){ for(int k = 0; k <= n; k++){ mint cur = 1; mint iv = 1; for(int l = 0; i*l+j <= n; l++){ int nxj = i*l+j; int nxk = k+l; int rem = n-j; dp[i][nxj][nxk] += dp[i-1][j][k]*comb(rem, i*l)*finv[l]*fac[i*l]*iv*cur; cur *= v[i]; iv *= finv[i]; } } } } for(int i = 0; i <= n-1; i++){ int j = n-i; cout << dp[n][n][j].val() << endl; } }