結果

問題 No.1873 Bracket Swapping
コンテスト
ユーザー cologne
提出日時 2022-03-13 14:54:01
言語 C++17(clang)
(17.0.6 + boost 1.87.0)
結果
AC  
実行時間 112 ms / 2,000 ms
コード長 2,330 bytes
コンパイル時間 2,979 ms
コンパイル使用メモリ 126,080 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-18 21:50:39
合計ジャッジ時間 3,627 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 27
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <string>
#include <vector>

#include <atcoder/modint>
using namespace std;

using Mint = atcoder::modint998244353;
using Mat = vector<vector<Mint>>;

Mat mul(Mat a, Mat b)
{
    int N = a.size();
    Mat 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][k] += a[i][j] * b[j][k];
    return c;
}

Mat I(int N)
{
    Mat a(N, vector<Mint>(N, 0));
    for (int i = 0; i < N; i++)
        a[i][i] = 1;
    return a;
}

Mat pow(Mat a, int b)
{
    if (b == 0)
        return I(a.size());
    if (b == 1)
        return a;
    Mat r = pow(a, b / 2);
    r = mul(r, r);
    if (b & 1)
        r = mul(r, a);
    return r;
}

int main()
{
    string S;
    int K;
    cin >> S >> K;

    int HN = S.size() / 2;
    Mat M(HN + 1, vector<Mint>(HN + 1, 0));
    for (int i = 0; i <= HN; i++)
    {
        // curr diff: 2*i
        M[i][i] = HN * (HN - 1) + 2 * i * (HN - i);
        if (i != HN)
            M[i + 1][i] = (HN - i) * (HN - i);
        if (i != 0)
            M[i - 1][i] = i * i;
    }
    M = pow(M, K);

    Mat DP(1, vector<Mint>(1, 1));
    for (int i = 0; i < (int)S.size(); i++)
    {
        Mat N(i + 2, vector<Mint>(i + 2));
        // prefix sum = j, total difference = k
        for (int j = 0; j <= i; j++)
            for (int k = 0; k <= i; k++)
            {
                // choose '('
                {
                    int new_j = j + 1;
                    int new_k = k + (S[i] != '(');
                    if (new_j >= 0 && new_k >= 0)
                        N[new_j][new_k] += DP[j][k];
                }
                // choose ')'
                {
                    int new_j = j - 1;
                    int new_k = k + (S[i] != ')');
                    if (new_j >= 0 && new_k >= 0)
                        N[new_j][new_k] += DP[j][k];
                }
            }
        DP = N;
    }

    Mat C(HN + 1, vector<Mint>(HN + 1));
    for (int i = 0; i <= HN; i++)
    {
        C[i][i] = C[i][0] = 1;
        for (int j = 1; j < i; j++)
            C[i][j] = C[i - 1][j - 1] + C[i - 1][j];
    }

    Mint ans = 0;
    for (int i = 0; i <= HN; i++)
        ans += M[i][0] * DP[0][2 * i] * (C[HN][i] * C[HN][i]).inv();
    cout << ans.val() << endl;
}
0