結果

問題 No.2472 Tea time in the grand garden
ユーザー 👑 SPD_9X2SPD_9X2
提出日時 2023-07-19 00:26:43
言語 PyPy3
(7.3.15)
結果
TLE  
(最新)
AC  
(最初)
実行時間 -
コード長 1,667 bytes
コンパイル時間 835 ms
コンパイル使用メモリ 82,588 KB
実行使用メモリ 179,916 KB
最終ジャッジ日時 2024-11-07 03:05:11
合計ジャッジ時間 3,894 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 39 ms
60,300 KB
testcase_01 AC 38 ms
52,832 KB
testcase_02 TLE -
testcase_03 -- -
testcase_04 -- -
testcase_05 -- -
testcase_06 -- -
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 -- -
testcase_30 -- -
testcase_31 -- -
testcase_32 -- -
testcase_33 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

"""

Tea Time (5) 想定解

定数倍高速化Ver

注: O(N+K^3) なので、通りません

"""

import sys
from sys import stdin

mod = 998244353

def modfac(n, MOD):
 
    f = 1
    factorials = [1]
    for m in range(1, n + 1):
        f *= m
        f %= MOD
        factorials.append(f)
    inv = pow(f, MOD - 2, MOD)
    invs = [1] * (n + 1)
    invs[n] = inv
    for m in range(n, 1, -1):
        inv *= m
        inv %= MOD
        invs[m - 1] = inv
    return factorials, invs

def modnCr(n,r):
    if n < 0 or r < 0 or n < r:
        return 0
    return fac[n] * inv[n-r] * inv[r] % mod

fac,inv = modfac(1000,mod)


N,K = map(int,stdin.readline().split())

# 制約チェック
assert 1 <= N <= 400
assert 0 <= K <= 400

if K == 0:
    print (1)
    sys.exit()

# 1つ目の括弧は置いた状態から開始
dp = [[[0],[0]] for l in range(2)]
dp[0][1][0] = 1 #last,s,p の順番

for i in range(2,2*K+1):

    ndp = [[[0] * i for nnum in range(min(K+1,i+1))] for l in range(2)]

    for put in range(2): #置く括弧 0 = ( , 1 = )
        for last in range(2):
            put_neq_last = (1 if last != put else 0)
            for ls in range(len(dp[0])):
                news = ls + (1,-1)[put]
                if 0 <= news < len(ndp[0]):
                    for lp in range(len(dp[0][0])):
                        newp = lp + put_neq_last
                        ndp[put][news][newp] = (ndp[put][news][newp] + dp[last][ls][lp]) % mod
    
    dp = ndp

    #print (dp)
ans = 0

for p in range(len(dp[0][0])):

    now = dp[1][0][p] * modnCr(2*K+N-p,N-p)
    #print (p,dp[0][p][1] , modnCr(2*K+N-p,N-p))
    ans += now
    ans %= mod

print (ans)
0