結果

問題 No.2788 4-33 Hard
コンテスト
ユーザー detteiuu
提出日時 2026-07-26 00:32:37
言語 PyPy3
(7.3.17)
コンパイル:
pypy3 -mpy_compile _filename_
実行:
pypy3 _filename_
結果
AC  
実行時間 151 ms / 2,000 ms
+ 577µs
コード長 1,182 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 242 ms
コンパイル使用メモリ 95,852 KB
実行使用メモリ 85,940 KB
最終ジャッジ日時 2026-07-26 00:33:13
合計ジャッジ時間 10,433 ms
ジャッジサーバーID
（参考情報） 		judge1_0 / judge2_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 52
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

MOD = 998244353

def inverse(n, d):
    return n * pow(d, -1, MOD) % MOD

def C(N, K):
    if N < K:
        return 0
    if N < K*2:
        K = N-K
    top = 1
    for i in range(N, N-K, -1):
        top *= i%MOD
        top %= MOD
    bottom = 1
    for i in range(1, K+1):
        bottom *= i%MOD
        bottom %= MOD
    return inverse(top, bottom)

O = [list(map(int, input().split())) for _ in range(5)]
X = [list(map(int, input().split())) for _ in range(5)]

dp = [[[0]*34 for _ in range(5)] for _ in range(9)]
dp[0][0][0] = 1
for i in range(5):
    for j in range(34):
        ndp = [[[0]*34 for _ in range(5)] for _ in range(9)]
        for k in range(9):
            for l in range(5):
                for m in range(34):
                    if dp[k][l][m] == 0: continue
                    for n in range(O[i][j]+1):
                        a, b, c = k+n, l+i*n, m+j*n
                        if 8 < a or 4 < b or 33 < c: break
                        ndp[a][b][c] += dp[k][l][m]*C(O[i][j], n)%MOD
                        ndp[a][b][c] %= MOD
        dp = ndp

ans = 0
for i in range(5):
    j = 4-i
    ans += dp[8][i][33]*(X[j][0]%MOD)%MOD
    ans %= MOD

print(ans)
0