結果

問題 No.1036 Make One With GCD 2
ユーザー 👑 obakyanobakyan
提出日時 2020-05-01 00:54:12
言語 Lua
(LuaJit 2.1.1696795921)
結果
TLE  
実行時間 -
コード長 1,859 bytes
コンパイル時間 170 ms
コンパイル使用メモリ 5,380 KB
実行使用メモリ 360,424 KB
最終ジャッジ日時 2023-08-23 07:07:19
合計ジャッジ時間 16,404 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1,279 ms
187,608 KB
testcase_01 AC 703 ms
137,676 KB
testcase_02 AC 660 ms
139,100 KB
testcase_03 AC 171 ms
49,668 KB
testcase_04 AC 308 ms
99,780 KB
testcase_05 AC 1 ms
4,380 KB
testcase_06 AC 1 ms
4,376 KB
testcase_07 AC 328 ms
63,132 KB
testcase_08 AC 276 ms
59,696 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 ~= -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