結果

問題 No.155 生放送とBGM
コンテスト
ユーザー te-sh
提出日時 2020-01-25 17:45:53
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 386 ms / 6,000 ms
コード長 3,379 bytes
コンパイル時間 3,299 ms
コンパイル使用メモリ 164,920 KB
実行使用メモリ 385,416 KB
最終ジャッジ日時 2024-06-22 04:29:55
合計ジャッジ時間 5,732 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 15
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

// URL: https://yukicoder.me/problems/no/155

import std.algorithm, std.array, std.container, std.math, std.range, std.typecons, std.string;
import std.conv;

version(unittest) {} else
void main()
{
  int N, L; io.getV(N, L);
  string[] st; io.getA(N, st);

  L *= 60;
  auto S = st.map!(sti => sti[0..2].to!int*60+sti[3..5].to!int).array;

  if (S.sum < L) { io.put(N); return; }

  auto dp1 = new long[][][](N+1, N+1, L); dp1[0][0][0] = 1;
  foreach (i, Si; S)
    foreach (j; 0..N+1)
      foreach (k; 0..L) {
        dp1[i+1][j][k] = dp1[i][j][k];
        if (j >= 1 && k >= Si)
          dp1[i+1][j][k] += dp1[i][j-1][k-Si];
      }

  auto fact = Fact!(double, false)(N+1);

  auto r = 0.0, dp2 = new long[][](N+1, L);
  foreach (i, Si; S) {
    foreach (j; 0..N+1)
      foreach (k; 0..L) {
        dp2[j][k] = dp1[N][j][k];
        if (j >= 1 && k >= Si)
          dp2[j][k] -= dp2[j-1][k-Si];
      }

    foreach (j; 0..N)
      foreach (k; max(0, L-Si)..L)
        r += dp2[j][k]*fact[j+1]*fact[N-1-j];
  }

  io.put(r/fact[N]);
}

struct Fact(T, bool inv = true)
{
  int n;
  T[] table;
  alias table this;
  static if (inv) T[] invTable;

  this(int n)
  {
    this.n = n;

    table = new T[](n+1); table[0] = 1;
    foreach (i; 1..n+1) table[i] = table[i-1]*i;

    static if (inv) {
      invTable = new T[](n+1); invTable[n] = T(1)/table[n];
      foreach_reverse (i; 1..n+1) invTable[i-1] = invTable[i]*i;
    }
  }

  static if (inv) {
    pure T perm(size_t a, size_t b) in { assert(n >= a && a >= b); } do
    { return table[a]*invTable[a-b]; }

    pure T combi(size_t a, size_t b) in { assert(n >= a && a >= b); } do
    { return table[a]*invTable[b]*invTable[a-b]; }

    pure T homo(size_t a, size_t b) in { assert(n >= a+b-1); } do
    { return combi(a+b-1, b); }
  }
}

auto io = IO!()();
import std.stdio;
struct IO(string floatFormat = "%.10f", string delimiter = " ", alias IN = stdin, alias OUT = stdout)
{
  import std.conv, std.format, std.meta, std.traits;
  alias assignable = hasAssignableElements;

  auto getV(T...)(ref T v) { foreach (ref w; v) get(w); }
  auto getA(T)(size_t n, ref T v) if (assignable!T) { v = new T(n); foreach (ref w; v) get(w); }
  auto getC(T...)(size_t n, ref T v) if (allSatisfy!(assignable, T))
  {
    foreach (ref w; v) w = new typeof(w)(n);
    foreach (i; 0..n) foreach (ref w; v) get(w[i]);
  }
  auto getM(T)(size_t r, size_t c, ref T v) if (assignable!T && assignable!(ElementType!T))
  {
    v = new T(r); foreach (ref w; v) getA(c, w);
  }

  auto put(bool flush = false, T...)(T v)
  {
    foreach (i, w; v) { putA(w); if (i < v.length-1) OUT.write(delimiter); }
    OUT.writeln;
    static if (flush) OUT.flush();
  }
  auto putB(S, T)(bool c, S t, T f) { if (c) put(t); else put(f); }
  auto putRaw(T...)(T v) { OUT.write(v); OUT.writeln; }

  private
  {
    dchar[] buf;
    auto sp = (new dchar[](0)).splitter;
    void nextLine() { IN.readln(buf); sp = buf.splitter; }
    auto get(T)(ref T v) { if (sp.empty) nextLine(); v = sp.front.to!T; sp.popFront(); }

    auto putR(T)(T v)
    {
      auto w = v;
      while (!w.empty) { putA(w.front); w.popFront(); if (!w.empty) OUT.write(delimiter); }
    }
    auto putA(T)(T v)
    {
      static if (isInputRange!T && !isSomeString!T) putR(v);
      else if (isFloatingPoint!T) OUT.write(format(floatFormat, v));
      else OUT.write(v);
    }
  }
}
0