結果

問題 No.425 ジャンケンの必勝法
ユーザー 👑 obakyanobakyan
提出日時 2021-08-12 23:19:38
言語 Lua
(LuaJit 2.1.1734355927)
結果
AC  
実行時間 4 ms / 2,000 ms
コード長 1,057 bytes
コンパイル時間 156 ms
コンパイル使用メモリ 5,376 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-02 05:40:33
合計ジャッジ時間 1,304 ms
ジャッジサーバーID
(参考情報)
judge3 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 5
other AC * 18
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

local p0, q = io.read("*n", "*n")
local mat = {}
for i = 1, 101 do
mat[i] = {}
for j = 1, 101 do
mat[i][j] = 0
end
end
local vec = {}
--[[
f(p) = p * 1/2 + p * 1/2 * f(p - q) + (1-p) * 1/3 + (1-p) * 1/3 * f(p + q)
f(p) - p/2 * f(p-q) - (1-p)/3 * f(p+q) = 1/3 + p/6
]]
for p = 0, 100 do
local left = p - q
if left < 0 then left = 0 end
local right = p + q
if 100 < right then right = 100 end
local r = p / 100
mat[p + 1][p + 1] = mat[p + 1][p + 1] + 1
mat[p + 1][left + 1] = mat[p + 1][left + 1] + (-r / 2)
mat[p + 1][right + 1] = mat[p + 1][right + 1] + (-(1 - r) / 3)
vec[p + 1] = 1 / 3 + r / 6
end
for i = 1, 101 do
local base = mat[i][i]
for k = i, 101 do
mat[i][k] = mat[i][k] / base
end
vec[i] = vec[i] / base
for j = i + 1, 101 do
local mul = mat[j][i]
for k = i, 101 do
mat[j][k] = mat[j][k] - mul * mat[i][k]
end
vec[j] = vec[j] - mul * vec[i]
end
end
for i = 101, 1, -1 do
for j = i - 1, 1, -1 do
vec[j] = vec[j] - mat[j][i] * vec[i]
end
end
print((1 + vec[p0 + 1]) / 3)
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0