結果

問題 No.301 サイコロで確率問題 (1)
ユーザー chocoruskchocorusk
提出日時 2020-09-22 01:04:35
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
MLE  
実行時間 -
コード長 728 bytes
コンパイル時間 106 ms
コンパイル使用メモリ 12,800 KB
実行使用メモリ 94,200 KB
最終ジャッジ日時 2024-06-25 00:30:41
合計ジャッジ時間 4,812 ms
ジャッジサーバーID
(参考情報)
judge1 / judge5
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 881 ms
44,228 KB
testcase_01 MLE -
権限があれば一括ダウンロードができます

ソースコード

diff #

import sys
read=sys.stdin.buffer.read
readline=sys.stdin.buffer.readline
readlines=sys.stdin.buffer.readlines

import numpy as np
mat=np.array([[0]*6 for i in range(6)], dtype=np.float64)
for i in range(6):
    mat[0][i]=1.0/6.0
for i in range(5):
    mat[i+1][i]=1.0
t=int(readline())
matps=[np.identity(6) for i in range(60)]
matps[0]=mat
for i in range(1, 60):
    matps[i]=np.dot(matps[i-1], matps[i-1])
for _ in range(t):
    n=int(readline())
    matp=np.identity(6)
    for j in range(60):
        if n&(1<<j):
            matp=np.dot(matp, matps[j])
    #matp=np.linalg.matrix_power(mat, n)
    x=matp[0][0]
    y=0
    for i in range(6):
        y+=2.0/7.0*i*matp[0][i]
    y+=2.0/7.0*n
    print('{:.13f}'.format(y/x))
0