結果
問題 | No.155 生放送とBGM |
ユーザー |
👑 |
提出日時 | 2015-02-17 23:57:27 |
言語 | D (dmd 2.109.1) |
結果 |
AC
|
実行時間 | 843 ms / 6,000 ms |
コード長 | 2,630 bytes |
コンパイル時間 | 3,027 ms |
コンパイル使用メモリ | 209,004 KB |
実行使用メモリ | 18,316 KB |
最終ジャッジ日時 | 2024-06-12 02:22:13 |
合計ジャッジ時間 | 8,745 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 15 |
ソースコード
import core.thread;import std.conv, std.stdio, std.string;import std.algorithm, std.array, std.bigint, std.container, std.math, std.range, std.regex;// Inputclass EOFException : Throwable { this() { super("EOF"); } }string[] tokens;string readToken() { for (; tokens.empty; ) { tokens = readln.split; if (tokens.empty && stdin.eof) throw new EOFException; } auto token = tokens[0];tokens.popFront; return token; }int readInt() { return to!int(readToken); }long readLong() { return to!long(readToken); }real readReal() { return to!real(readToken); }// chmin/chmaxvoid chmin(T)(ref T t, in T f) { if (t > f) t = f; }void chmax(T)(ref T t, in T f) { if (t < f) t = f; }// Pairstruct Pair(S, T) {S x; T y;int opCmp( const Pair p) const { return (x < p.x) ? -1 : (x > p.x) ? +1 : (y < p.y) ? -1 : (y > p.y) ? +1 : 0; }int opCmp(ref const Pair p) const { return (x < p.x) ? -1 : (x > p.x) ? +1 : (y < p.y) ? -1 : (y > p.y) ? +1 : 0; }string toString() const { return "(" ~ to!string(x) ~ ", " ~ to!string(y) ~ ")"; }}auto pair(S, T)(inout(S) x, inout(T) y) { return Pair!(S, T)(x, y); }// Arrayint binarySearch(T)(in T[] as, in bool delegate(T) test) { int low = -1, upp = as.length; for (; low + 1 < upp; ) { int mid = (low + upp) >> 1; (test(as[mid]) ? low : upp) = mid; } return upp; }int lowerBound(T)(in T[] as, in T val) { return as.binarySearch((T a) { return (a < val); }); }int upperBound(T)(in T[] as, in T val) { return as.binarySearch((T a) { return (a <= val); }); }T[] unique(T)(in T[] as) { T[] bs; foreach (a; as) if (bs.empty || bs[$ - 1] != a) bs ~= a; return bs; }int N, L;int[] A;void main(string[] args) {try {for (; ; ) {N = readInt;L = readInt * 60;A = new int[N];foreach (i; 0 .. N) {auto ss = readToken.split(":");A[i] = ss[0].to!int * 60 + ss[1].to!int;}auto fac = new real[N + 1];fac[0] = 1.0;foreach (k; 1 .. N + 1) {fac[k] = fac[k - 1] * k;}foreach (i; 0 .. N) {chmin(A[i], L);}auto dp = new long[][](N + 1, L * 2 + 1);dp[0][0] = 1;foreach (i; 0 .. N) {foreach_reverse (k; 0 .. N) foreach (x; 0 .. L * 2 + 1) {dp[k + 1][min(x + A[i], L * 2)] += dp[k][x];}}real ans = 0.0;foreach (i; 0 .. N) {foreach (k; 0 .. N) foreach (x; A[i] .. L * 2 + 1) {dp[k + 1][x] -= dp[k][x - A[i]];}foreach (k; 0 .. N) foreach (x; 0 .. L) {ans += (fac[k] * fac[N - 1 - k] / fac[N]) * dp[k][x];}foreach_reverse (k; 0 .. N) foreach (x; A[i] .. L * 2 + 1) {dp[k + 1][x] += dp[k][x - A[i]];}}writefln("%.10f", ans);}} catch (EOFException) {}}