結果

問題 No.210 探し物はどこですか？
ユーザー 👑 obakyan
提出日時 2021-08-12 21:54:04
言語 Lua
(LuaJit 2.1.1734355927)
結果
AC  
実行時間 1,516 ms / 2,000 ms
コード長 1,808 bytes
コンパイル時間 90 ms
コンパイル使用メモリ 5,120 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-10-02 03:13:04
合計ジャッジ時間 72,389 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 43
権限があれば一括ダウンロードができます

ソースコード

diff # 

local bls, brs = bit.lshift, bit.rshift

local Heapq = {}
Heapq.create = function(self, lt)
  self.lt = lt
  self.cnt = 0
  self.t = {}
end

Heapq.push = function(self, v)
  local hqlt = self.lt
  local hqt = self.t
  local c = self.cnt + 1
  self.cnt = c
  hqt[c] = v
  while 1 < c do
    local p = brs(c, 1)
    if hqlt(hqt[c], hqt[p]) then
      hqt[c], hqt[p] = hqt[p], hqt[c]
      c = p
    else
      break
    end
  end
end

Heapq.pop = function(self)
  local hqlt = self.lt
  local hqt = self.t
  local ret = hqt[1]
  local c = self.cnt
  hqt[1] = hqt[c]
  c = c - 1
  self.cnt = c
  local p = 1
  while true do
    local d1, d2 = p * 2, p * 2 + 1
    if c < d1 then break
    elseif c < d2 then
      if hqlt(hqt[d1], hqt[p]) then
        hqt[d1], hqt[p] = hqt[p], hqt[d1]
      end
      break
    else
      if hqlt(hqt[d1], hqt[d2]) then
        if hqlt(hqt[d1], hqt[p]) then
          hqt[d1], hqt[p] = hqt[p], hqt[d1]
          p = d1
        else break
        end
      else
        if hqlt(hqt[d2], hqt[p]) then
          hqt[d2], hqt[p] = hqt[p], hqt[d2]
          p = d2
        else break
        end
      end
    end
  end
  return ret
end

Heapq.new = function(lt)
  local obj = {}
  setmetatable(obj, {__index = Heapq})
  obj:create(lt)
  return obj
end

local n = io.read("*n")
local p, q = {}, {}
for i = 1, n do
  p[i] = io.read("*n") / 1000.0
end
for i = 1, n do
  q[i] = io.read("*n") / 100.0
end
local chance = {}
for i = 1, n do
  chance[i] = p[i] * q[i]
end
local hq = Heapq.new(function(x, y) return chance[x] > chance[y] end)
for i = 1, n do
  hq:push(i)
end
local ret = 0
local i = 0
while os.clock() < 1.5 do
  i = i + 1
  local idx = hq:pop()
  ret = ret + i * chance[idx]
  chance[idx] = chance[idx] * (1 - q[idx])
  hq:push(idx)
end
print(string.format("%.10f", ret))
0