結果

問題 No.1036 Make One With GCD 2
ユーザー 👑 obakyan
提出日時 2020-05-01 00:17:00
言語 Lua
(LuaJit 2.1.1734355927)
結果
TLE  
実行時間 -
コード長 1,977 bytes
コンパイル時間 119 ms
コンパイル使用メモリ 6,952 KB
実行使用メモリ 419,168 KB
最終ジャッジ日時 2024-12-18 00:06:48
合計ジャッジ時間 69,694 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 21 TLE * 20
権限があれば一括ダウンロードができます

ソースコード

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

local function getgcd(x, y)
  local z = 0LL
  if x < y then
    while 0LL < x do
      if y - x < x then
        x, y = y - x, x
      else
        x, y = y % x, x
      end
    end
    return y
  else
    while 0LL < y 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")
ffi.cdef[[
long long atoll(const char*);
]]
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] = 0
  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)
  for j = 1, n do
    local dst = j + jump
    if dst <= n then
      dp[i][j] = getgcd(dp[i - 1][j], dp[i - 1][dst])
    else
      dp[i][j] = 1LL
    end
  end
end
-- print(os.clock())

local function solvesub(l, r)
end
local function solve(l, r)
  if dp[1][l] == 1LL then return l end
  local ret = l
  local cur = 0LL
  local tmp = 0LL
  while true do
    -- print(cur, l, r)
    local stage = len[1 + r - l + 1] - 1
    tmp = getgcd(cur, dp[stage][l])
    if tmp == 1LL then
      if l == r then break end
      r = l + bls(1, stage - 1) - 2
    else
      cur = tmp
      ret = l + bls(1, stage - 1) - 1
      l = l + bls(1, stage - 1)
      if r < l then break end
    end
  end
  return ret + 1
end

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