結果

問題 No.1501 酔歩
ユーザー suisen
提出日時 2021-05-07 23:08:10
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 764 ms / 2,000 ms
コード長 652 bytes
コンパイル時間 464 ms
コンパイル使用メモリ 81,792 KB
実行使用メモリ 114,048 KB
最終ジャッジ日時 2024-09-15 11:14:27
合計ジャッジ時間 27,898 ms
ジャッジサーバーID
(参考情報)
judge4 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 53
権限があれば一括ダウンロードができます

ソースコード

diff #

from fractions import Fraction
N, K = map(int, input().split())
A = tuple(map(int, input().split()))
def solve():
    if K == N:
        print('1/1')
        return
    elif K == 1:
        print('0')
        return
    P = [Fraction(0) if i == 0 else Fraction(A[i + 1], A[i - 1] + A[i + 1]) for i in range(N - 1)]
    p2 = Fraction(1)
    p1 = Fraction(1) / P[1]
    for i in range(3, N):
        p0 = (p1 - (1 - P[i - 1]) * p2) / P[i - 1]
        p2 = p1
        p1 = p0
    ans = [0] * N
    ans[0] = 0
    ans[1] = 1 / p1
    for i in range(2, N):
        ans[i] = (ans[i - 1] - (1 - P[i - 1]) * ans[i - 2]) / P[i - 1]
    print(ans[K - 1])
solve()
0