結果

問題 No.2964 Obstruction Bingo
ユーザー srjywrdnprktsrjywrdnprkt
提出日時 2024-11-22 10:41:09
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,126 ms / 2,468 ms
コード長 1,567 bytes
コンパイル時間 2,341 ms
コンパイル使用メモリ 210,300 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-11-22 10:41:46
合計ジャッジ時間 36,244 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,248 KB
testcase_02 AC 8 ms
5,248 KB
testcase_03 AC 2 ms
5,248 KB
testcase_04 AC 4 ms
5,248 KB
testcase_05 AC 213 ms
5,248 KB
testcase_06 AC 520 ms
5,248 KB
testcase_07 AC 55 ms
5,248 KB
testcase_08 AC 508 ms
5,248 KB
testcase_09 AC 268 ms
5,248 KB
testcase_10 AC 646 ms
5,248 KB
testcase_11 AC 134 ms
5,248 KB
testcase_12 AC 782 ms
5,248 KB
testcase_13 AC 324 ms
5,248 KB
testcase_14 AC 567 ms
5,248 KB
testcase_15 AC 79 ms
5,248 KB
testcase_16 AC 154 ms
5,248 KB
testcase_17 AC 98 ms
5,248 KB
testcase_18 AC 245 ms
5,248 KB
testcase_19 AC 159 ms
5,248 KB
testcase_20 AC 542 ms
5,248 KB
testcase_21 AC 572 ms
5,248 KB
testcase_22 AC 94 ms
5,248 KB
testcase_23 AC 25 ms
5,248 KB
testcase_24 AC 14 ms
5,248 KB
testcase_25 AC 1,093 ms
5,248 KB
testcase_26 AC 1,094 ms
5,248 KB
testcase_27 AC 1,116 ms
5,248 KB
testcase_28 AC 1,124 ms
5,248 KB
testcase_29 AC 1,121 ms
5,248 KB
testcase_30 AC 1,117 ms
5,248 KB
testcase_31 AC 1,117 ms
5,248 KB
testcase_32 AC 1,113 ms
5,248 KB
testcase_33 AC 1,120 ms
5,248 KB
testcase_34 AC 1,119 ms
5,248 KB
testcase_35 AC 1,088 ms
5,248 KB
testcase_36 AC 1,120 ms
5,248 KB
testcase_37 AC 1,121 ms
5,248 KB
testcase_38 AC 1,092 ms
5,248 KB
testcase_39 AC 1,115 ms
5,248 KB
testcase_40 AC 1,118 ms
5,248 KB
testcase_41 AC 1,120 ms
5,248 KB
testcase_42 AC 1,115 ms
5,248 KB
testcase_43 AC 1,126 ms
5,248 KB
testcase_44 AC 1,086 ms
5,248 KB
testcase_45 AC 450 ms
5,248 KB
testcase_46 AC 347 ms
5,248 KB
testcase_47 AC 340 ms
5,248 KB
testcase_48 AC 574 ms
5,248 KB
testcase_49 AC 550 ms
5,248 KB
testcase_50 AC 776 ms
5,248 KB
testcase_51 AC 1,052 ms
5,248 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
#include <atcoder/modint>

using namespace std;
using namespace atcoder;
using ll = long long;
using mint = modint998244353;

int main(){
    cin.tie(nullptr);
    ios_base::sync_with_stdio(false);

    /*
       dp(i,j,k) = i回目でP_N=j, P_M=k(|k-j|<=L)である確率
    */
    ll L, K, x, c, d;
    mint ans1=0, ans2=0, sm=0;
    string S, T;
    cin >> L >> K >> S >> T;
    vector<mint> p(26);
    for (int i=0; i<26; i++){
        cin >> x;
        p[i] = x;
        sm += p[i];
    }
    sm = sm.inv();
    for (int i=0; i<26; i++) p[i] *= sm;

    vector dp(K+1, vector<mint>(K+1));
    dp[0][0] = 1;
    for (int i=0; i<K; i++){
        vector pd(K+1, vector<mint>(K+1));
        for (int j=0; j<=K; j++){
            for (int k=0; k<=K; k++){
                if (abs(j-k)>=L) continue;
                c = S[j%L]-'a';
                d = T[k%L]-'a';
                if (c == d){
                    if (j+1<=K && k+1<=K) pd[j+1][k+1] += dp[j][k] * p[c];
                    pd[j][k] += dp[j][k] * (-p[c]+1);
                }
                else{
                    if (j+1<=K) pd[j+1][k] += dp[j][k] * p[c];
                    if (k+1<=K) pd[j][k+1] += dp[j][k] * p[d];
                    pd[j][k] += dp[j][k] * (-p[c]-p[d]+1);
                }
            }
        }
        swap(dp, pd);
        for (int j=0; j<=K; j++){
            if (j-L>=0){
                ans1 += dp[j][j-L];
                ans2 += dp[j-L][j];
            }
        }
    }
    cout << ans1.val() << " " << ans2.val() << endl;

    return 0;
}
0