結果

問題 No.1006 Share an Integer
ユーザー 👑 obakyanobakyan
提出日時 2020-03-20 21:07:25
言語 Lua
(LuaJit 2.1.1696795921)
結果
AC  
実行時間 862 ms / 2,000 ms
コード長 1,361 bytes
コンパイル時間 310 ms
コンパイル使用メモリ 6,944 KB
実行使用メモリ 25,784 KB
最終ジャッジ日時 2024-05-08 18:25:11
合計ジャッジ時間 8,759 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 2 ms
6,944 KB
testcase_02 AC 2 ms
6,940 KB
testcase_03 AC 2 ms
6,944 KB
testcase_04 AC 3 ms
6,940 KB
testcase_05 AC 3 ms
6,940 KB
testcase_06 AC 3 ms
6,940 KB
testcase_07 AC 3 ms
6,948 KB
testcase_08 AC 3 ms
6,940 KB
testcase_09 AC 3 ms
6,944 KB
testcase_10 AC 3 ms
6,944 KB
testcase_11 AC 389 ms
23,092 KB
testcase_12 AC 802 ms
25,380 KB
testcase_13 AC 859 ms
25,748 KB
testcase_14 AC 862 ms
25,720 KB
testcase_15 AC 631 ms
25,784 KB
testcase_16 AC 275 ms
14,348 KB
testcase_17 AC 400 ms
24,180 KB
testcase_18 AC 642 ms
25,036 KB
testcase_19 AC 680 ms
24,864 KB
testcase_20 AC 659 ms
25,024 KB
testcase_21 AC 755 ms
25,700 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