結果
問題 |
No.1201 お菓子配り-4
|
ユーザー |
👑 |
提出日時 | 2021-07-29 22:24:25 |
言語 | Lua (LuaJit 2.1.1734355927) |
結果 |
AC
|
実行時間 | 2,716 ms / 4,000 ms |
コード長 | 1,058 bytes |
コンパイル時間 | 562 ms |
コンパイル使用メモリ | 5,504 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-09-14 18:49:04 |
合計ジャッジ時間 | 18,969 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 36 |
ソースコード
local mfl, mce = math.floor, math.ceil local mod = 1000000007 local half = 500000004 local z = 1000000000000000 local function bmul(x, y) if x * y < z then return (x * y) % mod end local x1, y1 = mfl(x / 31623), mfl(y / 31623) local x0, y0 = x - x1 * 31623, y - y1 * 31623 return (x1 * y1 * 14122 + (x1 * y0 + x0 * y1) * 31623 + x0 * y0) % mod end local function badd(x, y) return (x + y) % mod end local ret = 0 local function floorSum(n_m, m, a, b) while true do local n = n_m + m if m <= a then local v = mfl(a / m) ret = badd(ret, bmul(v, bmul(n, n - 1))) a = a % m end if m <= b then ret = badd(ret, bmul(2 * n, mfl(b / m))) b = b % m end if a == 0 then break end n_m, m, a, b = mfl((a * n_m + b) / m), a, m, (a * n_m + b) % m end end local n, m = io.read("*n", "*n") local a, b = {}, {} for i = 1, n do a[i] = io.read("*n") end for i = 1, m do b[i] = io.read("*n") end for i = 1, n do for j = 1, m do floorSum(1, b[j], a[i], 0) end end print(ret) -- print(os.clock())