結果
問題 | No.498 ワープクリスタル (給料日編) |
ユーザー | brthyyjp |
提出日時 | 2022-02-16 12:49:56 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,765 ms / 2,000 ms |
コード長 | 946 bytes |
コンパイル時間 | 555 ms |
コンパイル使用メモリ | 86,792 KB |
実行使用メモリ | 333,208 KB |
最終ジャッジ日時 | 2023-09-11 17:10:18 |
合計ジャッジ時間 | 22,070 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge15 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 115 ms
92,296 KB |
testcase_01 | AC | 116 ms
92,164 KB |
testcase_02 | AC | 116 ms
92,192 KB |
testcase_03 | AC | 116 ms
92,220 KB |
testcase_04 | AC | 212 ms
95,500 KB |
testcase_05 | AC | 118 ms
92,296 KB |
testcase_06 | AC | 1,735 ms
332,868 KB |
testcase_07 | AC | 634 ms
122,436 KB |
testcase_08 | AC | 1,752 ms
333,132 KB |
testcase_09 | AC | 1,765 ms
332,940 KB |
testcase_10 | AC | 125 ms
92,264 KB |
testcase_11 | AC | 177 ms
95,152 KB |
testcase_12 | AC | 120 ms
92,220 KB |
testcase_13 | AC | 189 ms
101,416 KB |
testcase_14 | AC | 121 ms
92,456 KB |
testcase_15 | AC | 142 ms
93,840 KB |
testcase_16 | AC | 122 ms
92,328 KB |
testcase_17 | AC | 137 ms
93,752 KB |
testcase_18 | AC | 1,732 ms
333,096 KB |
testcase_19 | AC | 1,756 ms
333,208 KB |
testcase_20 | AC | 1,731 ms
332,876 KB |
testcase_21 | AC | 1,746 ms
332,904 KB |
testcase_22 | AC | 1,737 ms
332,900 KB |
testcase_23 | AC | 1,748 ms
332,952 KB |
testcase_24 | AC | 1,746 ms
332,916 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 x, y, n in XYN: nx = defaultdict(lambda: 0) for a in range(0, n+1): for (cx, cy, s), v in dp.items(): 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)