結果

問題 No.1006 Share an Integer
ユーザー 👑 obakyanobakyan
提出日時 2020-03-20 21:07:25
言語 Lua
(LuaJit 2.1.1696795921)
結果
AC  
実行時間 860 ms / 2,000 ms
コード長 1,361 bytes
コンパイル時間 419 ms
コンパイル使用メモリ 7,204 KB
実行使用メモリ 25,804 KB
最終ジャッジ日時 2024-12-15 00:40:31
合計ジャッジ時間 8,816 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 2 ms
6,816 KB
testcase_02 AC 3 ms
6,816 KB
testcase_03 AC 2 ms
6,820 KB
testcase_04 AC 3 ms
6,820 KB
testcase_05 AC 3 ms
6,820 KB
testcase_06 AC 4 ms
6,816 KB
testcase_07 AC 3 ms
6,816 KB
testcase_08 AC 3 ms
6,816 KB
testcase_09 AC 3 ms
6,816 KB
testcase_10 AC 3 ms
6,820 KB
testcase_11 AC 386 ms
23,304 KB
testcase_12 AC 734 ms
25,424 KB
testcase_13 AC 860 ms
25,660 KB
testcase_14 AC 760 ms
25,676 KB
testcase_15 AC 718 ms
25,692 KB
testcase_16 AC 283 ms
14,304 KB
testcase_17 AC 414 ms
24,028 KB
testcase_18 AC 695 ms
24,824 KB
testcase_19 AC 609 ms
24,792 KB
testcase_20 AC 659 ms
25,000 KB
testcase_21 AC 756 ms
25,804 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

local mfl, mce = math.floor, math.ceil
local msq = math.sqrt
local mab = 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 function solve(x, primes)
  local ret = 1
  local prime_num = #primes
  local lim = mce(msq(x))
  local primepos = 1
  local dv = primes[primepos]
  while primepos <= prime_num and dv <= lim do
    if x % dv == 0 then
      local cnt = 1
      x = x / dv
      while x % dv == 0 do
        x = x / dv
        cnt = cnt + 1
      end
      lim = mce(msq(x))
      ret = ret * (cnt + 1)
    end
    if primepos == prime_num then break end
    primepos = primepos + 1
    dv = primes[primepos]
  end
  if x ~= 1 then
    ret = ret * 2
  end
  return ret
end

local x = io.read("*n")
local primes = getprimes(mce(msq(2000000)))
local t = {}
for i = 1, x - 1 do
  t[i] = i - solve(i, primes)
end
local minval, mins = mab(t[1] - t[x - 1]), {1}
for i = 2, x - 1 do
  local v = mab(t[i] - t[x - i])
  if v < minval then
    minval, mins = v, {i}
  elseif v == minval then
    table.insert(mins, i)
  end
end
for i = 1, #mins do
  print(mins[i] .. " " .. x - mins[i])
end
0