結果

問題 No.301 サイコロで確率問題 (1)
コンテスト
ユーザー chocorusk
提出日時 2020-09-22 01:22:39
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
WA  
実行時間 -
コード長 803 bytes
コンパイル時間 89 ms
コンパイル使用メモリ 12,800 KB
実行使用メモリ 92,676 KB
最終ジャッジ日時 2024-06-25 00:57:41
合計ジャッジ時間 4,003 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other WA * 1 MLE * 1
権限があれば一括ダウンロードができます

ソースコード

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