結果
問題 | No.1164 GCD Products hard |
ユーザー | 👑 obakyan |
提出日時 | 2021-08-10 11:23:24 |
言語 | Lua (LuaJit 2.1.1696795921) |
結果 |
AC
|
実行時間 | 1,324 ms / 2,500 ms |
コード長 | 1,640 bytes |
コンパイル時間 | 374 ms |
コンパイル使用メモリ | 5,336 KB |
実行使用メモリ | 141,800 KB |
最終ジャッジ日時 | 2023-10-23 18:49:37 |
合計ジャッジ時間 | 20,612 ms |
ジャッジサーバーID (参考情報) |
judge14 / judge12 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 819 ms
72,172 KB |
testcase_01 | AC | 1,014 ms
72,172 KB |
testcase_02 | AC | 729 ms
72,172 KB |
testcase_03 | AC | 240 ms
20,940 KB |
testcase_04 | AC | 235 ms
20,940 KB |
testcase_05 | AC | 874 ms
72,172 KB |
testcase_06 | AC | 996 ms
76,264 KB |
testcase_07 | AC | 1,029 ms
141,800 KB |
testcase_08 | AC | 983 ms
141,800 KB |
testcase_09 | AC | 711 ms
72,172 KB |
testcase_10 | AC | 236 ms
19,900 KB |
testcase_11 | AC | 773 ms
72,172 KB |
testcase_12 | AC | 1,090 ms
76,264 KB |
testcase_13 | AC | 643 ms
72,172 KB |
testcase_14 | AC | 598 ms
72,172 KB |
testcase_15 | AC | 1,058 ms
141,800 KB |
testcase_16 | AC | 662 ms
72,172 KB |
testcase_17 | AC | 858 ms
72,172 KB |
testcase_18 | AC | 739 ms
72,172 KB |
testcase_19 | AC | 258 ms
37,324 KB |
testcase_20 | AC | 415 ms
37,324 KB |
testcase_21 | AC | 1,148 ms
141,800 KB |
testcase_22 | AC | 2 ms
4,348 KB |
testcase_23 | AC | 977 ms
141,800 KB |
testcase_24 | AC | 1,324 ms
141,800 KB |
testcase_25 | AC | 1 ms
4,348 KB |
testcase_26 | AC | 1 ms
4,348 KB |
testcase_27 | AC | 1 ms
4,348 KB |
testcase_28 | AC | 1 ms
4,348 KB |
ソースコード
local mce, mfl, msq, mmi, mma, mab = math.ceil, math.floor, math.sqrt, math.min, math.max, math.abs local function getprimes(x) local primes = {} local allnums = {} for i = 1, x do allnums[i] = true end for i = 2, x do if allnums[i] then table.insert(primes, i) local lim = mfl(x / i) for j = 2, lim do allnums[j * i] = false end end end return primes end local mod6 = 1000000006 local function bmul6(x, y) local x1, y1 = mfl(x / 31623), mfl(y / 31623) local x0, y0 = x - x1 * 31623, y - y1 * 31623 return (x1 * y1 * 14123 + (x1 * y0 + x0 * y1) * 31623 + x0 * y0) % mod6 end local function modpow6(src, pow) local res = 1 while 0 < pow do if pow % 2 == 1 then res = bmul6(res, src) pow = pow - 1 end src = bmul6(src, src) pow = mfl(pow / 2) end return res end local mod = 1000000007 local function bmul(x, y) 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 modpow(src, pow) local res = 1 while 0 < pow do if pow % 2 == 1 then res = bmul(res, src) pow = pow - 1 end src = bmul(src, src) pow = mfl(pow / 2) end return res end local a, b, n = io.read("*n", "*n", "*n") local ans = 1 local primes = getprimes(b) for i = 1, #primes do local p = primes[i] local rawp = p while true do local c1 = mfl(b / p) - mfl((a - 1) / p) local v = modpow6(c1, n) ans = bmul(ans, modpow(rawp, v)) p = p * rawp if b < p then break end end end print(ans)