結果

問題 No.498 ワープクリスタル (給料日編)
ユーザー qibqib
提出日時 2023-02-21 16:01:05
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 367 ms / 2,000 ms
コード長 825 bytes
コンパイル時間 320 ms
コンパイル使用メモリ 87,160 KB
実行使用メモリ 77,372 KB
最終ジャッジ日時 2023-09-29 08:19:45
合計ジャッジ時間 7,351 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 76 ms
71,080 KB
testcase_01 AC 75 ms
71,284 KB
testcase_02 AC 76 ms
71,008 KB
testcase_03 AC 76 ms
71,192 KB
testcase_04 AC 350 ms
76,640 KB
testcase_05 AC 78 ms
70,964 KB
testcase_06 AC 367 ms
76,688 KB
testcase_07 AC 364 ms
76,904 KB
testcase_08 AC 356 ms
76,828 KB
testcase_09 AC 359 ms
76,480 KB
testcase_10 AC 75 ms
70,828 KB
testcase_11 AC 120 ms
76,932 KB
testcase_12 AC 75 ms
71,316 KB
testcase_13 AC 91 ms
76,748 KB
testcase_14 AC 76 ms
71,164 KB
testcase_15 AC 90 ms
76,200 KB
testcase_16 AC 78 ms
75,764 KB
testcase_17 AC 85 ms
76,644 KB
testcase_18 AC 365 ms
76,896 KB
testcase_19 AC 358 ms
77,104 KB
testcase_20 AC 354 ms
77,232 KB
testcase_21 AC 357 ms
77,372 KB
testcase_22 AC 357 ms
76,684 KB
testcase_23 AC 360 ms
76,716 KB
testcase_24 AC 354 ms
76,848 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