結果
問題 | No.2874 Gunegune Tree |
ユーザー | 寝癖 |
提出日時 | 2024-07-14 18:00:07 |
言語 | C++23 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 1,071 ms / 2,000 ms |
コード長 | 1,991 bytes |
コンパイル時間 | 3,549 ms |
コンパイル使用メモリ | 261,456 KB |
実行使用メモリ | 14,156 KB |
最終ジャッジ日時 | 2024-07-14 18:05:48 |
合計ジャッジ時間 | 18,839 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
6,816 KB |
testcase_01 | AC | 2 ms
6,816 KB |
testcase_02 | AC | 2 ms
6,940 KB |
testcase_03 | AC | 75 ms
6,940 KB |
testcase_04 | AC | 434 ms
7,808 KB |
testcase_05 | AC | 157 ms
6,940 KB |
testcase_06 | AC | 3 ms
6,944 KB |
testcase_07 | AC | 1,008 ms
13,492 KB |
testcase_08 | AC | 467 ms
8,068 KB |
testcase_09 | AC | 45 ms
6,940 KB |
testcase_10 | AC | 851 ms
11,992 KB |
testcase_11 | AC | 658 ms
10,076 KB |
testcase_12 | AC | 204 ms
6,944 KB |
testcase_13 | AC | 975 ms
13,444 KB |
testcase_14 | AC | 808 ms
11,648 KB |
testcase_15 | AC | 136 ms
6,944 KB |
testcase_16 | AC | 1,042 ms
13,828 KB |
testcase_17 | AC | 115 ms
6,944 KB |
testcase_18 | AC | 125 ms
6,944 KB |
testcase_19 | AC | 84 ms
6,940 KB |
testcase_20 | AC | 110 ms
6,940 KB |
testcase_21 | AC | 331 ms
6,940 KB |
testcase_22 | AC | 492 ms
8,448 KB |
testcase_23 | AC | 310 ms
6,940 KB |
testcase_24 | AC | 410 ms
7,424 KB |
testcase_25 | AC | 703 ms
10,696 KB |
testcase_26 | AC | 432 ms
7,680 KB |
testcase_27 | AC | 1,023 ms
13,604 KB |
testcase_28 | AC | 229 ms
6,940 KB |
testcase_29 | AC | 585 ms
9,344 KB |
testcase_30 | AC | 463 ms
8,192 KB |
testcase_31 | AC | 343 ms
6,944 KB |
testcase_32 | AC | 1,071 ms
14,156 KB |
ソースコード
#include <bits/stdc++.h> #include <atcoder/modint> using namespace std; using namespace atcoder; using mint = modint998244353; const mint zero = 0, one = 1; // 状態i->jに遷移する確率 const vector<vector<mint>> P = { {zero, one/5, one/5, one/5, one/5, one/5}, {zero, one/2, one/2, zero , zero , zero }, {zero, one/3, one/3, one/3, zero , zero }, {zero, zero , one/3, one/3, one/3, zero }, {zero, zero , zero , one/3, one/3, one/3}, {zero, zero , zero , zero , one/2, one/2} }; // 単位行列 const vector<vector<mint>> E = { {one, zero, zero, zero, zero, zero}, {zero, one, zero, zero, zero, zero}, {zero, zero, one, zero, zero, zero}, {zero, zero, zero, one, zero, zero}, {zero, zero, zero, zero, one, zero}, {zero, zero, zero, zero, zero, one} }; // 行列積 vector<vector<mint>> operator*(const vector<vector<mint>> &A, const vector<vector<mint>> &B) { int N = A.size(); vector<vector<mint>> C(N, vector<mint>(N, 0)); for (int i = 0; i < N; i++) { for (int j = 0; j < N; j++) { for (int k = 0; k < N; k++) { C[i][j] += A[i][k]*B[k][j]; } } } return C; } int main() { int N; cin >> N; // 0: 初期状態, 1: 右, 2: 右上, 3: 上, 4: 左上, 5: 左 vector<vector<mint>> dp(N+1, vector<mint>(6)); // P^i, P^{i+1} vector<vector<mint>> Pi = E, Pi1 = P; for (int i = 0; i < N; i++) { for (int j = 1; j < 6; j++) { for (int k = 0; k < 6; k++) { // (i+1)日目の状態がjであるという条件のもとで、i日目の状態がkである確率 mint pk = P[k][j]*Pi[0][k]/Pi1[0][j]; // 1<j<5ならば高さが1増える dp[i+1][j] += (dp[i][k]+(1<j&&j<5))*pk; } } Pi = Pi1; Pi1 = Pi*P; } mint ans = 0; for (int i = 0; i < 6; i++) ans += dp[N][i]*Pi[0][i]; cout << ans.val() << endl; }