結果
問題 | No.425 ジャンケンの必勝法 |
ユーザー | 👑 obakyan |
提出日時 | 2021-08-12 23:19:38 |
言語 | Lua (LuaJit 2.1.1696795921) |
結果 |
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 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 4 ms
5,248 KB |
testcase_01 | AC | 4 ms
5,248 KB |
testcase_02 | AC | 4 ms
5,248 KB |
testcase_03 | AC | 4 ms
5,248 KB |
testcase_04 | AC | 4 ms
5,248 KB |
testcase_05 | AC | 4 ms
5,248 KB |
testcase_06 | AC | 4 ms
5,248 KB |
testcase_07 | AC | 4 ms
5,248 KB |
testcase_08 | AC | 4 ms
5,248 KB |
testcase_09 | AC | 4 ms
5,248 KB |
testcase_10 | AC | 4 ms
5,248 KB |
testcase_11 | AC | 4 ms
5,248 KB |
testcase_12 | AC | 4 ms
5,248 KB |
testcase_13 | AC | 3 ms
5,248 KB |
testcase_14 | AC | 4 ms
5,248 KB |
testcase_15 | AC | 4 ms
5,248 KB |
testcase_16 | AC | 4 ms
5,248 KB |
testcase_17 | AC | 3 ms
5,248 KB |
testcase_18 | AC | 4 ms
5,248 KB |
testcase_19 | AC | 4 ms
5,248 KB |
testcase_20 | AC | 4 ms
5,248 KB |
testcase_21 | AC | 4 ms
5,248 KB |
testcase_22 | AC | 4 ms
5,248 KB |
ソースコード
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)