結果

問題 No.155 生放送とBGM
ユーザー vwxyz
提出日時 2021-09-24 21:52:11
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 596 ms / 6,000 ms
コード長 670 bytes
コンパイル時間 274 ms
コンパイル使用メモリ 12,544 KB
実行使用メモリ 64,404 KB
最終ジャッジ日時 2024-07-05 10:27:28
合計ジャッジ時間 8,487 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 15
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
readline=sys.stdin.readline
import numpy as np

N,L=map(int,readline().split())
L*=60
S=[]
for s in readline().split():
    mm,ss=map(int,s.split(":"))
    S.append(mm*60+ss)
ans=0
S_sum=sum(S)
if L>=S_sum:
    ans=N
else:
    fact=np.zeros(51)
    fact[0]=1.0
    for i in range(1,51):
        fact[i]=fact[i-1]*i
    dp=np.zeros((N+1,18001))
    dp[0,0]=1
    for s in S:
        dp[1:,s:]+=dp[:-1,:-s]
    for s in S:
        for i in range(1,N+1):
            dp[i,s:]-=dp[i-1,:-s]
        for i in range(N):
            ans+=fact[i]*dp[i,:L].sum()*fact[N-i-1]
        for i in range(N,0,-1):
            dp[i,s:]+=dp[i-1,:-s]
    ans/=fact[N]
print(ans)
0