結果

問題 No.425 ジャンケンの必勝法
コンテスト
ユーザー rpy3cpp
提出日時 2016-09-23 00:02:39
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 31 ms / 2,000 ms
コード長 526 bytes
コンパイル時間 96 ms
コンパイル使用メモリ 12,672 KB
実行使用メモリ 10,752 KB
最終ジャッジ日時 2024-11-17 19:43:23
合計ジャッジ時間 1,410 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 5
other AC * 18
権限があれば一括ダウンロードができます

ソースコード

diff # 

def solve(p, q):
    dp = [0.0] * 101
    p_win = 1.0/3.0
    dp[p] = 1.0/3.0
    eps = 10**(-9)
    while sum(dp) > eps:
        dp, p_win_delta = update_dp(dp, q)
        p_win += p_win_delta
    return p_win

def update_dp(dp, q):
    new_dp = [0.0] * 101
    pdelta = 0.0
    for i, pi in enumerate(dp):
        pdelta += (i/200 + (100 - i)/300) * pi
        new_dp[max(0, i - q)] += pi*i/200
        new_dp[min(100, i + q)] += pi*(100-i)/300
    return new_dp, pdelta

p, q = map(int, input().split())

print(solve(p, q))
0