結果

問題 No.1036 Make One With GCD 2
コンテスト
ユーザー 👑 obakyan
提出日時 2020-05-01 00:33:03
言語 Lua
(LuaJit 2.1.1734355927)
結果
MLE  
実行時間 -
コード長 1,794 bytes
コンパイル時間 260 ms
コンパイル使用メモリ 6,820 KB
実行使用メモリ 537,060 KB
最終ジャッジ日時 2024-12-18 00:43:41
合計ジャッジ時間 68,943 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 21 TLE * 19 MLE * 1
権限があれば一括ダウンロードができます

ソースコード

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 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 = getgcd(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