結果

問題 No.2964 Obstruction Bingo
ユーザー eve__fuyukieve__fuyuki
提出日時 2024-11-17 16:34:30
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 54 ms / 2,468 ms
コード長 1,374 bytes
コンパイル時間 2,658 ms
コンパイル使用メモリ 213,132 KB
実行使用メモリ 14,208 KB
最終ジャッジ日時 2024-11-17 16:34:36
合計ジャッジ時間 5,925 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,248 KB
testcase_02 AC 2 ms
5,248 KB
testcase_03 AC 2 ms
5,248 KB
testcase_04 AC 2 ms
5,248 KB
testcase_05 AC 2 ms
5,248 KB
testcase_06 AC 25 ms
8,192 KB
testcase_07 AC 5 ms
5,248 KB
testcase_08 AC 19 ms
6,912 KB
testcase_09 AC 3 ms
5,248 KB
testcase_10 AC 13 ms
5,760 KB
testcase_11 AC 20 ms
7,040 KB
testcase_12 AC 39 ms
11,264 KB
testcase_13 AC 7 ms
5,248 KB
testcase_14 AC 6 ms
5,248 KB
testcase_15 AC 3 ms
5,248 KB
testcase_16 AC 13 ms
5,504 KB
testcase_17 AC 9 ms
5,248 KB
testcase_18 AC 9 ms
5,248 KB
testcase_19 AC 18 ms
6,656 KB
testcase_20 AC 21 ms
7,424 KB
testcase_21 AC 39 ms
11,392 KB
testcase_22 AC 14 ms
6,016 KB
testcase_23 AC 8 ms
5,248 KB
testcase_24 AC 2 ms
5,248 KB
testcase_25 AC 54 ms
14,080 KB
testcase_26 AC 54 ms
14,080 KB
testcase_27 AC 53 ms
13,952 KB
testcase_28 AC 53 ms
14,208 KB
testcase_29 AC 53 ms
14,080 KB
testcase_30 AC 53 ms
14,080 KB
testcase_31 AC 52 ms
14,208 KB
testcase_32 AC 52 ms
14,080 KB
testcase_33 AC 54 ms
13,952 KB
testcase_34 AC 54 ms
14,080 KB
testcase_35 AC 52 ms
14,080 KB
testcase_36 AC 54 ms
13,952 KB
testcase_37 AC 53 ms
14,080 KB
testcase_38 AC 53 ms
14,080 KB
testcase_39 AC 53 ms
13,952 KB
testcase_40 AC 53 ms
13,952 KB
testcase_41 AC 53 ms
14,080 KB
testcase_42 AC 53 ms
14,080 KB
testcase_43 AC 53 ms
14,080 KB
testcase_44 AC 53 ms
14,080 KB
testcase_45 AC 5 ms
5,248 KB
testcase_46 AC 2 ms
5,248 KB
testcase_47 AC 2 ms
5,248 KB
testcase_48 AC 7 ms
5,248 KB
testcase_49 AC 7 ms
5,248 KB
testcase_50 AC 17 ms
6,820 KB
testcase_51 AC 48 ms
13,184 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

void fast_io() {
	ios_base::sync_with_stdio(false);
	cin.tie(nullptr);
}
#include <atcoder/modint>
using mint = atcoder::modint998244353;
int main() {
	fast_io();
	int l, k;
	cin >> l >> k;
	string s, t;
	cin >> s >> t;
	vector<int> a(26);
	long long a_sum = 0;
	for (int i = 0; i < 26; i++) {
		cin >> a[i];
		a_sum += a[i];
	}
	vector<mint> p(26);
	{
		mint asum_inv = mint(a_sum).inv();
		for (int i = 0; i < 26; i++) {
			p[i] = mint(a[i]) * asum_inv;
		}
	}
	// j: pn % l, x: pm - pn + l
	vector dp(k + 1, vector(l, vector<mint>(2 * l + 1, 0)));
	dp[0][0][l] = 1;
	for (int i = 0; i < k; i++) {
		for (int j = 0; j < l; j++) {
			dp[i + 1][j][0] += dp[i][j][0];
			dp[i + 1][j][2 * l] += dp[i][j][2 * l];
			for (int x = 1; x < 2 * l; x++) {
				int pn = j, pm = (j + x) % l;
				int pnn = s[pn] - 'a', pmm = t[pm] - 'a';
				if (pnn == pmm) {
					dp[i + 1][(j + 1) % l][x] += dp[i][j][x] * p[pnn];
					dp[i + 1][j][x] += dp[i][j][x] * (1 - p[pnn]);
				} else {
					dp[i + 1][(j + 1) % l][x - 1] += dp[i][j][x] * p[pnn];
					dp[i + 1][j][x + 1] += dp[i][j][x] * p[pmm];
					dp[i + 1][j][x] += dp[i][j][x] * (1 - p[pnn] - p[pmm]);
				}
			}
		}
	}
	mint ans_n = 0, ans_m = 0;
	for (int j = 0; j < l; j++) {
		ans_n += dp[k][j][0];
		ans_m += dp[k][j][2 * l];
	}
	cout << ans_n.val() << " " << ans_m.val() << endl;
}
0