結果
問題 | No.813 ユキちゃんの冒険 |
ユーザー |
|
提出日時 | 2021-07-02 00:41:52 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
WA
|
実行時間 | - |
コード長 | 735 bytes |
コンパイル時間 | 78 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 44,632 KB |
最終ジャッジ日時 | 2024-06-28 13:24:50 |
合計ジャッジ時間 | 15,747 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 21 WA * 5 |
ソースコード
import sysinput = lambda : sys.stdin.readline().rstrip()sys.setrecursionlimit(2*10**5+10)write = lambda x: sys.stdout.write(x+"\n")debug = lambda x: sys.stderr.write(x+"\n")writef = lambda x: print("{:.12f}".format(x))import numpy as npn = int(input())p = float(input())q = float(input())if p==1:x = 1elif q==0:x = pelse:dp0 = [0]*ndp1 = [0]*ndp0[0] = np.array((1,0), dtype=np.float64)dp1[0] = np.array((0,1), dtype=np.float64)for i in range(1,n):dp0[i] = (dp0[i-1]- p*dp1[i-1]) / qdp1[i] = q*dp1[i-1] + p*dp0[i]v = np.linalg.norm(dp0[i])dp0[i], dp1[i] = dp0[i]/v, dp1[i]/va0,b0 = dp0[n-1]a1,b1 = dp1[n-1]x = (b1-b0) / (a0-a1)print(x)