結果

問題 No.75 回数の期待値の問題
ユーザー tktk_snsntktk_snsn
提出日時 2020-06-05 18:43:28
言語 Python3
(3.9.3 + numpy 1.14.5 + scipy 1.1.0)
結果
AC  
実行時間 154 ms / 5,000 ms
コード長 371 Byte
コンパイル時間 95 ms
使用メモリ 31,020 KB
最終ジャッジ日時 2020-09-09 20:35:23
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
使用メモリ
testcase_00 AC 145 ms
28,600 KB
testcase_01 AC 144 ms
28,808 KB
testcase_02 AC 146 ms
28,444 KB
testcase_03 AC 149 ms
28,856 KB
testcase_04 AC 149 ms
28,664 KB
testcase_05 AC 151 ms
28,792 KB
testcase_06 AC 151 ms
28,624 KB
testcase_07 AC 149 ms
28,444 KB
testcase_08 AC 147 ms
28,652 KB
testcase_09 AC 145 ms
28,612 KB
testcase_10 AC 146 ms
28,608 KB
testcase_11 AC 146 ms
28,864 KB
testcase_12 AC 153 ms
28,864 KB
testcase_13 AC 151 ms
28,748 KB
testcase_14 AC 147 ms
28,464 KB
testcase_15 AC 147 ms
28,572 KB
testcase_16 AC 147 ms
29,848 KB
testcase_17 AC 148 ms
28,856 KB
testcase_18 AC 154 ms
28,852 KB
testcase_19 AC 151 ms
31,020 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

import numpy as np
K = int(input())
A = np.zeros((K+1, K+1), dtype = np.float64)
for i in range(K + 1):
    A[i, i] = 1.0
    p = 1.0
    for j in range(1, 7):
        if i + j > K:
            break
        A[i, i + j] -= 1.0 / 6.0
        p -= 1.0 / 6.0
    if i < K:
        A[i, 0] -= p

B = np.ones(K+1, np.float64)
B[K] = 0.0

E = np.linalg.solve(A, B)
print(E[0])
0