結果
問題 | No.425 ジャンケンの必勝法 |
ユーザー |
![]() |
提出日時 | 2020-09-26 20:06:56 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 514 ms / 2,000 ms |
コード長 | 461 bytes |
コンパイル時間 | 82 ms |
コンパイル使用メモリ | 12,544 KB |
実行使用メモリ | 45,084 KB |
最終ジャッジ日時 | 2024-06-29 17:41:02 |
合計ジャッジ時間 | 13,523 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 5 |
other | AC * 18 |
ソースコード
import sys read=sys.stdin.buffer.read readline=sys.stdin.buffer.readline readlines=sys.stdin.buffer.readlines import numpy as np p0, q=map(int, readline().split()) mat=np.zeros((101, 101), dtype=np.float64) b=np.zeros(101, dtype=np.float64) for p in range(101): p1=max(0, p-q) p2=min(100, p+q) mat[p][p1]-=p/200 mat[p][p]+=1 mat[p][p2]-=(100-p)/300 b[p]=p/200+(100-p)/300 x=np.linalg.solve(mat, b) print('{:.7f}'.format(1/3+1/3*x[p0]))