結果
問題 |
No.2472 Tea time in the grand garden
|
ユーザー |
👑 ![]() |
提出日時 | 2023-07-19 00:14:53 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,550 bytes |
コンパイル時間 | 330 ms |
コンパイル使用メモリ | 81,920 KB |
実行使用メモリ | 297,472 KB |
最終ジャッジ日時 | 2024-11-07 03:04:45 |
合計ジャッジ時間 | 3,884 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 WA * 1 TLE * 1 -- * 1 |
other | -- * 30 |
ソースコード
""" Tea Time 想定解 定数倍無視ver 注: ミスっています WAです。 """ 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 # 1つ目の括弧は置いた状態から開始 dp = [ [ [0,0] ] , [ [0,0] ] ] dp[1][0][0] = 1 for i in range(2,2*K+1): ndp = [[[0,0] for pnum in range(i)] for nnum in range(min(K+1,i+1))] for ls in range(len(dp)): for lp in range(len(dp[0])): for last in range(2): for put in range(2): #置く括弧 0 = ( , 1 = ) news = ls + (1,-1)[put] newp = lp + (1 if last != put else 0) if 0 <= news < len(ndp): ndp[news][newp][put] += dp[ls][lp][last] ndp[news][newp][put] %= mod dp = ndp ans = 0 for p in range(len(dp[0])): now = dp[0][p][1] * 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)