結果

問題 No.2788 4-33 Hard
ユーザー mattu34mattu34
提出日時 2024-06-21 07:00:22
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 257 ms / 2,000 ms
コード長 3,330 bytes
コンパイル時間 316 ms
コンパイル使用メモリ 82,848 KB
実行使用メモリ 91,148 KB
最終ジャッジ日時 2024-06-21 07:00:41
合計ジャッジ時間 15,041 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 218 ms
90,588 KB
testcase_01 AC 165 ms
90,624 KB
testcase_02 AC 193 ms
90,380 KB
testcase_03 AC 222 ms
90,632 KB
testcase_04 AC 193 ms
91,140 KB
testcase_05 AC 219 ms
90,800 KB
testcase_06 AC 191 ms
90,492 KB
testcase_07 AC 162 ms
90,288 KB
testcase_08 AC 183 ms
91,016 KB
testcase_09 AC 200 ms
91,148 KB
testcase_10 AC 221 ms
90,344 KB
testcase_11 AC 223 ms
90,716 KB
testcase_12 AC 235 ms
90,628 KB
testcase_13 AC 251 ms
90,996 KB
testcase_14 AC 223 ms
90,504 KB
testcase_15 AC 218 ms
90,628 KB
testcase_16 AC 215 ms
90,500 KB
testcase_17 AC 218 ms
90,632 KB
testcase_18 AC 252 ms
90,636 KB
testcase_19 AC 242 ms
90,500 KB
testcase_20 AC 241 ms
90,636 KB
testcase_21 AC 216 ms
90,504 KB
testcase_22 AC 221 ms
90,636 KB
testcase_23 AC 238 ms
90,636 KB
testcase_24 AC 221 ms
90,508 KB
testcase_25 AC 218 ms
90,764 KB
testcase_26 AC 221 ms
90,496 KB
testcase_27 AC 218 ms
90,616 KB
testcase_28 AC 257 ms
90,504 KB
testcase_29 AC 222 ms
90,628 KB
testcase_30 AC 221 ms
90,640 KB
testcase_31 AC 226 ms
90,472 KB
testcase_32 AC 234 ms
90,608 KB
testcase_33 AC 248 ms
90,600 KB
testcase_34 AC 216 ms
90,760 KB
testcase_35 AC 212 ms
90,632 KB
testcase_36 AC 213 ms
90,376 KB
testcase_37 AC 216 ms
90,372 KB
testcase_38 AC 215 ms
90,504 KB
testcase_39 AC 251 ms
90,520 KB
testcase_40 AC 223 ms
90,624 KB
testcase_41 AC 217 ms
90,500 KB
testcase_42 AC 215 ms
90,376 KB
testcase_43 AC 215 ms
90,584 KB
testcase_44 AC 238 ms
90,756 KB
testcase_45 AC 216 ms
90,504 KB
testcase_46 AC 214 ms
90,484 KB
testcase_47 AC 223 ms
90,612 KB
testcase_48 AC 217 ms
90,612 KB
testcase_49 AC 248 ms
90,496 KB
testcase_50 AC 217 ms
90,336 KB
testcase_51 AC 224 ms
90,468 KB
testcase_52 AC 219 ms
90,472 KB
testcase_53 AC 223 ms
90,468 KB
testcase_54 AC 220 ms
90,600 KB
testcase_55 AC 249 ms
90,472 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

from collections import *
import sys
import heapq
import bisect
import itertools
from functools import lru_cache
from types import GeneratorType
from fractions import Fraction
import math
import copy
import random

# sys.setrecursionlimit(int(1e7))
# @lru_cache(maxsize=None) # CPython特化
# @bootstrap # PyPy特化(こっちのほうが速い) yield dfs(), yield Noneを忘れずに


def bootstrap(f, stack=[]):  # yield
    def wrappedfunc(*args, **kwargs):
        if stack:
            return f(*args, **kwargs)
        else:
            to = f(*args, **kwargs)
            while True:
                if type(to) is GeneratorType:
                    stack.append(to)
                    to = next(to)
                else:
                    stack.pop()
                    if not stack:
                        break
                    to = stack[-1].send(to)
            return to

    return wrappedfunc


dxdy1 = ((0, 1), (0, -1), (1, 0), (-1, 0))  # 上下右左
dxdy2 = (
    (0, 1),
    (0, -1),
    (1, 0),
    (-1, 0),
    (1, 1),
    (-1, -1),
    (1, -1),
    (-1, 1),
)  # 8方向すべて
dxdy3 = ((0, 1), (1, 0))  # 右 or 下
dxdy4 = ((1, 1), (1, -1), (-1, 1), (-1, -1))  # 斜め
INF = float("inf")
_INF = 1 << 60
MOD = 998244353
mod = 998244353
MOD2 = 10**9 + 7
mod2 = 10**9 + 7
# memo : len([a,b,...,z])==26
# memo : 2^20 >= 10^6
# 小数の計算を避ける : x/y -> (x*big)//y  ex:big=10**9
# @:小さい文字, ~:大きい文字,None: 空の文字列
# ユークリッドの互除法：gcd(x,y)=gcd(x,y-x)
# memo : d 桁以下の p 進表記を用いると p^d-1 以下のすべての
#        非負整数を表現することができる
# memo : (X,Y) -> (X+Y,X−Y) <=> 点を原点を中心に45度回転し、√2倍に拡大
# memo : (x,y)のx正から見た偏角をラジアンで(-πからπ]：　math.atan2(y, x)
# memo : a < bのとき ⌊a⌋ ≦ ⌊b⌋

input = lambda: sys.stdin.readline().rstrip()
mi = lambda: map(int, input().split())
li = lambda: list(mi())
ii = lambda: int(input())
py = lambda: print("Yes")
pn = lambda: print("No")
pf = lambda: print("First")
ps = lambda: print("Second")

# 階乗 & 逆元計算
factorial = [1]
inverse = [1]
for i in range(1, 10):
    factorial.append(factorial[-1] * i % MOD)
    inverse.append(pow(factorial[-1], MOD - 2, MOD))


# 組み合わせ計算
def nCr(n, r):
    if n < r or r < 0:
        return 0
    elif r == 0:
        return 1
    tmp = 1
    for x in range(n, n - r, -1):
        tmp *= x
        tmp %= MOD
    return tmp * inverse[r] % MOD


O = [li() for _ in range(5)]
X = [li() for _ in range(5)]

dp = [[[0] * 34 for _ in range(5)] for _ in range(9)]
dp[0][0][0] = 1

for a in range(5):
    for b in range(34):
        for j in range(7, -1, -1):
            for x in range(4, -1, -1):
                for y in range(33, -1, -1):
                    for r in range(1, min(9, O[a][b] + 1)):
                        if x + r * a <= 4 and y + r * b <= 33 and j + r <= 8:
                            dp[j + r][x + r * a][y + r * b] += dp[j][x][y] * nCr(
                                O[a][b], r
                            )
                            dp[j + r][x + r * a][y + r * b] %= MOD
ans = 0
for a in range(5):
    for b in range(1):
        ans += dp[8][4 - a][33] * X[a][b]
        ans %= MOD
print(ans)
0