結果

問題 No.1036 Make One With GCD 2
コンテスト
ユーザー obakyan
提出日時 2020-05-01 00:54:12
言語 Lua
(LuaJit 2.1.1765228720)
結果
TLE  
実行時間 -
コード長 1,859 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 405 ms
コンパイル使用メモリ 6,944 KB
実行使用メモリ 416,348 KB
最終ジャッジ日時 2024-12-21 07:39:55
合計ジャッジ時間 70,806 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 22 TLE * 19
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

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 ~= -x do
    x, y = y % x, x
  end
  return y
end
local function getgcd2(x, y)
  if x < y then
    while x ~= -x do
      if y < x + x then
        x, y = y - x, x
      else
        x, y = y % x, x
      end
    end
    return y
  else
    while y ~= -y do
    -- 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")
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