結果

問題 No.2964 Obstruction Bingo
ユーザー eve__fuyukieve__fuyuki
提出日時 2024-11-17 16:24:51
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 1,305 bytes
コンパイル時間 2,382 ms
コンパイル使用メモリ 214,404 KB
実行使用メモリ 14,208 KB
最終ジャッジ日時 2024-11-17 16:24:57
合計ジャッジ時間 4,643 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 2 ms
6,816 KB
testcase_02 AC 1 ms
6,816 KB
testcase_03 AC 1 ms
6,816 KB
testcase_04 AC 1 ms
6,816 KB
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 WA -
testcase_20 WA -
testcase_21 WA -
testcase_22 WA -
testcase_23 WA -
testcase_24 WA -
testcase_25 WA -
testcase_26 WA -
testcase_27 WA -
testcase_28 WA -
testcase_29 WA -
testcase_30 WA -
testcase_31 WA -
testcase_32 WA -
testcase_33 WA -
testcase_34 WA -
testcase_35 WA -
testcase_36 WA -
testcase_37 WA -
testcase_38 WA -
testcase_39 WA -
testcase_40 WA -
testcase_41 WA -
testcase_42 WA -
testcase_43 WA -
testcase_44 WA -
testcase_45 WA -
testcase_46 AC 2 ms
6,816 KB
testcase_47 AC 2 ms
6,816 KB
testcase_48 AC 6 ms
6,816 KB
testcase_49 WA -
testcase_50 WA -
testcase_51 WA -
権限があれば一括ダウンロードができます

ソースコード

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);
	int 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++) {
			for (int x = 1; x < 2 * l; x++) {
				int pn = j, pm = (j + x) % l;
				if (s[pn] == t[pm]) {
					dp[i + 1][(j + 1) % l][x] += dp[i][j][x] * p[s[pn] - 'a'];
					dp[i][j][x] += dp[i][j][x] * (1 - p[s[pn] - 'a']);
				} else {
					dp[i + 1][(j + 1) % l][x - 1] +=
						dp[i][j][x] * p[s[pn] - 'a'];
					dp[i + 1][j][x + 1] += dp[i][j][x] * p[t[pm] - 'a'];
					dp[i + 1][j][x] +=
						dp[i][j][x] * (1 - p[s[pn] - 'a'] - p[t[pm] - 'a']);
				}
			}
		}
	}
	mint ans_n = 0, ans_m = 0;
	for (int i = 0; i <= k; i++) {
		ans_n += dp[i][0][0];
		ans_m += dp[i][0][2 * l];
	}
	cout << ans_n.val() << " " << ans_m.val() << endl;
}
0