結果
問題 | No.498 ワープクリスタル (給料日編) |
ユーザー | brthyyjp |
提出日時 | 2022-02-16 13:06:10 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,810 ms / 2,000 ms |
コード長 | 1,060 bytes |
コンパイル時間 | 383 ms |
コンパイル使用メモリ | 86,788 KB |
実行使用メモリ | 321,100 KB |
最終ジャッジ日時 | 2023-09-11 17:11:23 |
合計ジャッジ時間 | 14,031 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge15 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 92 ms
71,568 KB |
testcase_01 | AC | 89 ms
71,296 KB |
testcase_02 | AC | 93 ms
71,172 KB |
testcase_03 | AC | 90 ms
71,200 KB |
testcase_04 | AC | 192 ms
79,364 KB |
testcase_05 | AC | 96 ms
75,968 KB |
testcase_06 | AC | 1,681 ms
320,664 KB |
testcase_07 | AC | 652 ms
108,416 KB |
testcase_08 | AC | 1,810 ms
321,100 KB |
testcase_09 | AC | 628 ms
176,824 KB |
testcase_10 | AC | 93 ms
76,112 KB |
testcase_11 | AC | 149 ms
78,236 KB |
testcase_12 | AC | 92 ms
76,084 KB |
testcase_13 | AC | 190 ms
86,052 KB |
testcase_14 | AC | 97 ms
75,852 KB |
testcase_15 | AC | 128 ms
78,132 KB |
testcase_16 | AC | 96 ms
75,952 KB |
testcase_17 | AC | 116 ms
77,756 KB |
testcase_18 | AC | 1,196 ms
286,984 KB |
testcase_19 | AC | 1,440 ms
286,352 KB |
testcase_20 | AC | 1,752 ms
320,744 KB |
testcase_21 | AC | 1,419 ms
286,096 KB |
testcase_22 | AC | 1,068 ms
287,264 KB |
testcase_23 | AC | 1,651 ms
320,640 KB |
testcase_24 | AC | 1,797 ms
321,072 KB |
ソースコード
N = 100 mod = 10**9+7 fac = [1]*(N+1) finv = [1]*(N+1) for i in range(N): fac[i+1] = fac[i] * (i+1) % mod finv[-1] = pow(fac[-1], mod-2, mod) for i in reversed(range(N)): finv[i] = finv[i+1] * (i+1) % mod def cmb1(n, r): if r <0 or r > n: return 0 r = min(r, n-r) return fac[n] * finv[r] * finv[n-r] % mod from collections import defaultdict gx, gy, k = map(int, input().split()) XYN = [] for i in range(k): x,y,n = map(int, input().split()) XYN.append((x, y, n)) dp = defaultdict(lambda: 0) dp[(0, 0, 0)] = 1 for i, (x, y, n) in enumerate(XYN): nx = defaultdict(lambda: 0) for a in range(0, n+1): for (cx, cy, s), v in dp.items(): if abs(cx-gx) > 10**4*15*(k-i) or abs(cy-gy) > 10**4*15*(k-i): continue ncx, ncy = cx+a*x, cy+a*y ns = s+a nx[(ncx,ncy,ns)] += v*cmb1(ns,a) nx[(ncx,ncy,ns)] %= mod dp = nx ans = 0 for (cx, cy, s), v in dp.items(): if cx == gx and cy == gy: ans += v ans %= mod print(ans)