結果
問題 | No.498 ワープクリスタル (給料日編) |
ユーザー |
👑 |
提出日時 | 2021-06-19 17:40:08 |
言語 | Lua (LuaJit 2.1.1734355927) |
結果 |
AC
|
実行時間 | 64 ms / 2,000 ms |
コード長 | 1,533 bytes |
コンパイル時間 | 207 ms |
コンパイル使用メモリ | 5,120 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-06-22 22:12:24 |
合計ジャッジ時間 | 1,875 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 21 |
ソースコード
local mod = 1000000007local mfl = math.floorlocal function bmul(x, y)local x1, y1 = mfl(x / 31623), mfl(y / 31623)local x0, y0 = x - x1 * 31623, y - y1 * 31623return (x1 * y1 * 14122 + (x1 * y0 + x0 * y1) * 31623 + x0 * y0) % modendlocal function badd(x, y)return (x + y) % modendlocal function bsub(x, y)return x < y and x - y + mod or x - yendlocal function modpow(src, pow)local res = 1while 0 < pow doif pow % 2 == 1 thenres = bmul(res, src)pow = pow - 1endsrc = bmul(src, src)pow = mfl(pow / 2)endreturn resendlocal function modinv(src)return modpow(src, mod - 2)endlocal function getComb(n, k)local ret = 1local inv = 1for i = 1, k doret = bmul(ret, n + 1 - i)inv = bmul(inv, i)endreturn bmul(ret, modinv(inv))endlocal gx, gy, n = io.read("*n", "*n", "*n")local x, y, c = {}, {}, {}for i = 1, n dox[i], y[i], c[i] = io.read("*n", "*n", "*n")endlocal tot = 1for i = 1, n dotot = tot * (c[i] + 1)endlocal ret = 0local box = {}for i = 0, tot - 1 dolocal ti = ifor j = 1, n dobox[j] = ti % (c[j] + 1)ti = mfl(ti / (c[j] + 1))endlocal vx, vy = 0, 0for j = 1, n dovx = vx + x[j] * box[j]vy = vy + y[j] * box[j]endif vx == gx and vy == gy thenlocal r = 1local cnt = 0for j = 1, n docnt = cnt + box[j]endfor j = 1, n dor = bmul(r, getComb(cnt, box[j]))cnt = cnt - box[j]endret = badd(ret, r)endendprint(ret)