結果

問題 No.1036 Make One With GCD 2
ユーザー 👑 obakyanobakyan
提出日時 2020-05-01 00:56:36
言語 Lua
(LuaJit 2.1.1696795921)
結果
MLE  
実行時間 -
コード長 1,838 bytes
コンパイル時間 201 ms
コンパイル使用メモリ 5,068 KB
実行使用メモリ 534,940 KB
最終ジャッジ日時 2023-08-23 07:11:41
合計ジャッジ時間 14,031 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 MLE -
testcase_01 AC 801 ms
137,640 KB
testcase_02 AC 629 ms
140,164 KB
testcase_03 AC 173 ms
49,836 KB
testcase_04 AC 306 ms
99,772 KB
testcase_05 AC 1 ms
4,380 KB
testcase_06 AC 1 ms
4,376 KB
testcase_07 AC 308 ms
63,088 KB
testcase_08 AC 256 ms
59,612 KB
testcase_09 TLE -
testcase_10 TLE -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
testcase_30 -- -
testcase_31 -- -
testcase_32 -- -
testcase_33 -- -
testcase_34 -- -
testcase_35 -- -
testcase_36 -- -
testcase_37 -- -
testcase_38 -- -
testcase_39 -- -
testcase_40 -- -
testcase_41 -- -
testcase_42 -- -
testcase_43 -- -
testcase_44 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

local mfl, mce = math.floor, math.ceil
local mmi, mma = math.min, math.max
local bls, brs = bit.lshift, bit.rshift
local ffi = require("ffi")
local C = ffi.C
ffi.cdef[[
long long atoll(const char*);
]]
local function getgcd1(x, y)
  while x ~= 0LL do
    x, y = y % x, x
  end
  return y
end
local function getgcd2(x, y)
  if x < y then
    while x ~= 0LL do
      if y < x + x then
        x, y = y - x, x
      else
        x, y = y % x, x
      end
    end
    return y
  else
    while y ~= 0LL do
      if x < y + y then
        y, x = x - y, y
      else
        y, x = x % y, y
      end
    end
    return x
  end
end

local n = io.read("*n", "*l")
local roundcnt = 1
do
  local mul = 1
  while mul < n do
    roundcnt = roundcnt + 1
    mul = mul * 2
  end
end
local len = {}
do
  local tmp, c = 1, 1
  for i = 1, n + 1 do
    if tmp < i then tmp, c = tmp * 2, c + 1 end
    len[i] = c
  end
end

local dp = {}
for i = 1, roundcnt do
  dp[i] = {}
  for j = 1, n do
    dp[i][j] = 1LL
  end
end

local s = io.read()
do
  local i = 1
  for w in s:gmatch("%d+") do
    dp[1][i] = C.atoll(w)
    i = i + 1
  end
end
-- print(os.clock())
for i = 2, roundcnt do
  local jump = bls(1, i - 2)
  local dpi = dp[i]
  local dpim1 = dp[i - 1]
  for j = 1, n - jump do
    dpi[j] = getgcd1(dpim1[j], dpim1[j + jump])
  end
end
-- print(os.clock())

local function solve(l, r)
  if dp[1][l] == 1LL then return l end
  local cur = 0LL
  local tmp = 0LL
  while 1 < r - l do
    local stage = len[1 + r - l] - 1
    tmp = getgcd2(cur, dp[stage][l])
    if tmp == 1LL then
      r = l + bls(1, stage - 1) - 1
    else
      cur = tmp
      l = l + bls(1, stage - 1) - 1
    end
  end
  return l + 1
end

local ret = 0
local z = n + 1
for i = n, 1, -1 do
  local p = solve(i, z)
  ret = ret + mma(0, n + 1 - p)
end
print(ret)
-- print(os.clock())
0