結果

問題 No.155 生放送とBGM
ユーザー te-sh
提出日時 2017-04-28 16:03:29
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 436 ms / 6,000 ms
コード長 992 bytes
コンパイル時間 2,884 ms
コンパイル使用メモリ 157,260 KB
実行使用メモリ 491,884 KB
最終ジャッジ日時 2024-06-12 18:58:15
合計ジャッジ時間 5,800 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 15
権限があれば一括ダウンロードができます

ソースコード

diff # 

import std.algorithm, std.conv, std.range, std.stdio, std.string;
import std.bitmanip;  // BitArray

void main()
{
  auto rd = readln.split.to!(int[]), n = rd[0], l = rd[1] * 60;
  auto si = readln.split.map!(s => s[0..2].to!int * 60 + s[3..5].to!int).array;

  if (l >= si.sum) {
    writeln(n);
    return;
  }

  auto dp1 = new long[][](l, n + 1);
  dp1[0][0] = 1;
  foreach (i; 0..n)
    foreach_reverse (j; si[i]..l)
      dp1[j][1..$] += dp1[j - si[i]][0..$-1];

  auto dp2 = new long[][][](n, l, n + 1);
  foreach (i; 0..n)
    foreach (j; 0..l) {
      dp2[i][j][] = dp1[j][];
      if (j >= si[i])
        dp2[i][j][1..$] -= dp2[i][j-si[i]][0..$-1];
    }

  auto fact = new real[](n + 1);
  fact[0] = 1;
  foreach (i; 1..n+1)
    fact[i] = fact[i-1] * i;

  auto r = real(0);
  foreach (j; 0..n)
    foreach (m; max(0, l-si[j])..l)
      foreach (k; 0..n) {
        auto a = dp2[j][m][k];
        r += a * fact[k] * fact[n-k-1] * (k+1);
      }

  writefln("%.10f", r / fact[n]);
}
0