結果
問題 | 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.ceillocal mod = 1000000007local half = 500000004local z = 1000000000000000local function bmul(x, y)if x * y < z then return (x * y) % mod endlocal 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) % mod endlocal ret = 0local function floorSum(n_m, m, a, b)while true dolocal n = n_m + mif m <= a thenlocal v = mfl(a / m)ret = badd(ret, bmul(v, bmul(n, n - 1)))a = a % mendif m <= b thenret = badd(ret, bmul(2 * n, mfl(b / m)))b = b % mendif a == 0 then break endn_m, m, a, b = mfl((a * n_m + b) / m), a, m, (a * n_m + b) % mendendlocal n, m = io.read("*n", "*n")local a, b = {}, {}for i = 1, n doa[i] = io.read("*n")endfor i = 1, m dob[i] = io.read("*n")endfor i = 1, n dofor j = 1, m dofloorSum(1, b[j], a[i], 0)endendprint(ret)-- print(os.clock())