結果

問題 No.301 サイコロで確率問題 (1)
ユーザー chocoruskchocorusk
提出日時 2020-09-22 01:23:40
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
MLE  
(最新)
AC  
(最初)
実行時間 -
コード長 787 bytes
コンパイル時間 88 ms
コンパイル使用メモリ 12,928 KB
実行使用メモリ 92,528 KB
最終ジャッジ日時 2024-06-25 00:59:07
合計ジャッジ時間 2,882 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 858 ms
43,472 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())
    if n>=10**13:
        print(n)
        continue
    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*float(i)*matp[0][i]
    y+=2.0/7.0*n
    print('{:.13f}'.format(y/x))
0