結果

問題 No.301 サイコロで確率問題 (1)
ユーザー chocoruskchocorusk
提出日時 2020-09-22 01:04:35
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
WA  
実行時間 -
コード長 728 bytes
コンパイル時間 151 ms
コンパイル使用メモリ 10,752 KB
実行使用メモリ 58,960 KB
最終ジャッジ日時 2023-09-07 06:05:07
合計ジャッジ時間 4,757 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

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

ソースコード

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