結果

問題 No.301 サイコロで確率問題 (1)
ユーザー chocoruskchocorusk
提出日時 2020-09-22 01:22:39
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
WA  
実行時間 -
コード長 803 bytes
コンパイル時間 259 ms
コンパイル使用メモリ 10,980 KB
実行使用メモリ 29,828 KB
最終ジャッジ日時 2023-09-07 06:35:16
合計ジャッジ時間 2,965 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge13
このコードへのチャレンジ（β）

テストケース

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

ソースコード

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
    print(x, y)
    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