結果

問題 No.1873 Bracket Swapping
ユーザー minimumminimum
提出日時 2023-06-09 19:41:10
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 1,744 bytes
コンパイル時間 154 ms
コンパイル使用メモリ 82,176 KB
実行使用メモリ 85,684 KB
最終ジャッジ日時 2024-06-10 11:53:04
合計ジャッジ時間 5,661 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 37 ms
58,368 KB
testcase_01 AC 36 ms
52,864 KB
testcase_02 AC 58 ms
65,920 KB
testcase_03 AC 905 ms
77,028 KB
testcase_04 AC 398 ms
76,964 KB
testcase_05 AC 384 ms
76,828 KB
testcase_06 TLE -
testcase_07 -- -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

class Matrix():

    def __init__(self, MOD=-1):
        self.MOD = MOD

    def mul(self, a, b):
        L, M, N = len(a), len(b), len(b[0])
        assert len(a[0]) == M
        c = [[0] * N for _ in range(L)]
        for i in range(L):
            for j in range(N):
                for k in range(M):
                    c[i][j] += a[i][k] * b[k][j]
                    if self.MOD != -1:
                        c[i][j] %= self.MOD
        return c
    
    def pow(self, x, n):
        y = [[0] * len(x) for _ in range(len(x))]
        for i in range(len(x)):
            y[i][i] = 1
        while n > 0:
            if n & 1:
                y = self.mul(x, y)
            x = self.mul(x, x)
            n >>= 1
        return y

MOD = 998244353
mat = Matrix(MOD)

S = input()
N = len(S) // 2
K = int(input())

dp = [[0] * (2 * N + 1) for _ in range(2 * N + 1)]
dp[0][0] = 1
for i in range(2 * N):
    ndp = [[0] * (2 * N + 1) for _ in range(2 * N + 1)]
    for j in range(2 * N + 1):
        for k in range(2 * N + 1):
            if S[i] == '(':
                if j > 0: ndp[j][k] += dp[j - 1][k]
                if k > 0 and j < 2 * N: ndp[j][k] += dp[j + 1][k - 1]
            else:
                if j > 0 and k > 0: ndp[j][k] += dp[j - 1][k - 1]
                if j < 2 * N: ndp[j][k] += dp[j + 1][k]
            ndp[j][k] %= MOD
    dp = ndp
res1 = dp[0]

A = [[0] * (2 * N + 1) for _ in range(2 * N + 1)]
for i in range(N + 1):
    A[2 * i][2 * i] = N * (2 * N - 1) - i * i - (N - i) * (N - i)
    if i > 0:
        A[2 * i - 2][2 * i] = i * i
    if i < N:
        A[2 * i + 2][2 * i] = (N - i) * (N - i)
B = mat.pow(A, K)
res2 = B[0]

ans = 0
for i in range(2 * N + 1):
    ans += res1[i] * res2[i]
    ans %= MOD
print(ans)
0