結果
問題 | No.498 ワープクリスタル (給料日編) |
ユーザー | brthyyjp |
提出日時 | 2022-02-16 13:02:51 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,869 ms / 2,000 ms |
コード長 | 1,062 bytes |
コンパイル時間 | 395 ms |
コンパイル使用メモリ | 86,904 KB |
実行使用メモリ | 333,648 KB |
最終ジャッジ日時 | 2023-09-11 17:11:05 |
合計ジャッジ時間 | 19,311 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge11 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 122 ms
92,116 KB |
testcase_01 | AC | 127 ms
92,192 KB |
testcase_02 | AC | 122 ms
92,292 KB |
testcase_03 | AC | 121 ms
92,368 KB |
testcase_04 | AC | 220 ms
95,444 KB |
testcase_05 | AC | 123 ms
92,288 KB |
testcase_06 | AC | 1,733 ms
325,468 KB |
testcase_07 | AC | 696 ms
123,832 KB |
testcase_08 | AC | 1,869 ms
333,648 KB |
testcase_09 | AC | 676 ms
193,508 KB |
testcase_10 | AC | 123 ms
92,480 KB |
testcase_11 | AC | 186 ms
94,924 KB |
testcase_12 | AC | 122 ms
92,408 KB |
testcase_13 | AC | 228 ms
103,016 KB |
testcase_14 | AC | 123 ms
92,468 KB |
testcase_15 | AC | 158 ms
94,656 KB |
testcase_16 | AC | 127 ms
92,352 KB |
testcase_17 | AC | 149 ms
93,896 KB |
testcase_18 | AC | 1,254 ms
278,456 KB |
testcase_19 | AC | 1,515 ms
314,672 KB |
testcase_20 | AC | 1,821 ms
329,740 KB |
testcase_21 | AC | 1,505 ms
314,684 KB |
testcase_22 | AC | 1,107 ms
279,180 KB |
testcase_23 | AC | 1,670 ms
314,880 KB |
testcase_24 | AC | 1,825 ms
331,820 KB |
ソースコード
N = 10**6 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)