結果
| 問題 |
No.210 探し物はどこですか?
|
| コンテスト | |
| ユーザー |
👑 |
| 提出日時 | 2021-08-12 21:52:48 |
| 言語 | Lua (LuaJit 2.1.1734355927) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 1,795 bytes |
| コンパイル時間 | 248 ms |
| コンパイル使用メモリ | 7,072 KB |
| 実行使用メモリ | 6,824 KB |
| 最終ジャッジ日時 | 2024-10-02 03:10:34 |
| 合計ジャッジ時間 | 3,820 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 31 WA * 12 |
ソースコード
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 lim = 200000
for i = 1, lim do
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))