結果

問題 No.76 回数の期待値で練習
コンテスト
ユーザー Mao-beta
提出日時 2024-03-04 23:10:04
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 119 ms / 5,000 ms
コード長 1,334 bytes
コンパイル時間 356 ms
コンパイル使用メモリ 82,364 KB
実行使用メモリ 91,904 KB
最終ジャッジ日時 2024-09-29 17:34:03
合計ジャッジ時間 1,231 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 2
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
import math
import bisect
from heapq import heapify, heappop, heappush
from collections import deque, defaultdict, Counter
from functools import lru_cache
from itertools import accumulate, combinations, permutations, product

sys.setrecursionlimit(1000000)
MOD = 10 ** 9 + 7
MOD99 = 998244353

input = lambda: sys.stdin.readline().strip()
NI = lambda: int(input())
NMI = lambda: map(int, input().split())
NLI = lambda: list(NMI())
SI = lambda: input()
SMI = lambda: input().split()
SLI = lambda: list(SMI())
EI = lambda m: [NLI() for _ in range(m)]


def main():
    T = NI()
    M = 10**6
    E = [0.0] * (M+7)
    E[M-1] = 1.0000000000000000
    E[M-2] = 1.0833333333333333
    E[M-3] = 1.2569444444444444
    E[M-4] = 1.5353009259259260
    E[M-5] = 1.6915991512345676
    E[M-6] = 2.0513639724794235

    P = [0.0] * 7
    for i in range(2, 7):
        # P[i-1]を順に求める
        # E[M-i] = 1 + E[M-i+1]*P[1] + E[M-i+2]*P[2] * ...
        e = E[M-i]-1
        a = E[M-1]
        for x in range(1, i-1):
            e -= E[M-i+x] * P[x]
        P[i-1] = e / a
    P[6] = 1-sum(P)
    # print(P)

    for i in range(M-7, -1, -1):
        E[i] = 1
        for x in range(1, 7):
            E[i] += P[x] * E[i+x]
    for _ in range(T):
        N = NI()
        print(E[M-N])


if __name__ == "__main__":
    main()
0