結果
問題 | No.425 ジャンケンの必勝法 |
ユーザー |
👑 |
提出日時 | 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 |
ソースコード
local p0, q = io.read("*n", "*n")local mat = {}for i = 1, 101 domat[i] = {}for j = 1, 101 domat[i][j] = 0endendlocal 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 dolocal left = p - qif left < 0 then left = 0 endlocal right = p + qif 100 < right then right = 100 endlocal r = p / 100mat[p + 1][p + 1] = mat[p + 1][p + 1] + 1mat[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 / 6endfor i = 1, 101 dolocal base = mat[i][i]for k = i, 101 domat[i][k] = mat[i][k] / baseendvec[i] = vec[i] / basefor j = i + 1, 101 dolocal mul = mat[j][i]for k = i, 101 domat[j][k] = mat[j][k] - mul * mat[i][k]endvec[j] = vec[j] - mul * vec[i]endendfor i = 101, 1, -1 dofor j = i - 1, 1, -1 dovec[j] = vec[j] - mat[j][i] * vec[i]endendprint((1 + vec[p0 + 1]) / 3)