結果

問題 No.498 ワープクリスタル (給料日編)
ユーザー qibqib
提出日時 2023-02-21 16:01:05
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 335 ms / 2,000 ms
コード長 825 bytes
コンパイル時間 165 ms
コンパイル使用メモリ 82,208 KB
実行使用メモリ 76,380 KB
最終ジャッジ日時 2024-07-22 02:56:16
合計ジャッジ時間 5,829 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 41 ms
52,772 KB
testcase_01 AC 41 ms
53,088 KB
testcase_02 AC 40 ms
53,548 KB
testcase_03 AC 40 ms
53,020 KB
testcase_04 AC 314 ms
75,564 KB
testcase_05 AC 42 ms
52,680 KB
testcase_06 AC 335 ms
76,380 KB
testcase_07 AC 329 ms
75,964 KB
testcase_08 AC 320 ms
76,272 KB
testcase_09 AC 326 ms
75,832 KB
testcase_10 AC 41 ms
52,596 KB
testcase_11 AC 90 ms
75,628 KB
testcase_12 AC 41 ms
52,648 KB
testcase_13 AC 60 ms
68,408 KB
testcase_14 AC 42 ms
54,372 KB
testcase_15 AC 55 ms
63,644 KB
testcase_16 AC 47 ms
59,500 KB
testcase_17 AC 52 ms
62,820 KB
testcase_18 AC 331 ms
75,944 KB
testcase_19 AC 328 ms
76,100 KB
testcase_20 AC 321 ms
76,316 KB
testcase_21 AC 328 ms
75,712 KB
testcase_22 AC 319 ms
75,944 KB
testcase_23 AC 324 ms
76,260 KB
testcase_24 AC 322 ms
76,272 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

MOD = 10 ** 9 + 7

gx, gy, k = map(int, input().split())

m = 75
inv = [1 for _ in range(m + 1)]
fact = [1 for _ in range(m + 1)]
invfact = [1 for _ in range(m + 1)]
for x in range(2, m + 1):
  inv[x] = (- (MOD // x) * inv[MOD % x]) % MOD
  fact[x] = (fact[x - 1] * x) % MOD
  invfact[x] = (invfact[x - 1] * inv[x]) % MOD

items = []
cnt = [1]
for _ in range(k):
  xi, yi, ni = map(int, input().split())
  items.append((xi, yi, ni))
  cnt.append(cnt[-1] * (ni + 1))

ans = 0
for pat in range(cnt[-1]):
  x = 0
  y = 0
  q = []
  s = 0
  for i in range(k):
    xi, yi, _ = items[i]
    ci = (pat % cnt[i + 1]) // cnt[i]
    q.append(ci)
    s += ci
    x += (ci * xi)
    y += (ci * yi)

  if x == gx and y == gy:
    c = fact[s]
    for i in range(k):
      c = (c * invfact[q[i]]) % MOD
    ans = (ans + c) % MOD

print(ans)
0