結果

問題 No.155 生放送とBGM
ユーザー shimomireshimomire
提出日時 2015-02-18 01:31:08
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 571 ms / 6,000 ms
コード長 5,709 bytes
コンパイル時間 1,119 ms
コンパイル使用メモリ 119,796 KB
実行使用メモリ 32,116 KB
最終ジャッジ日時 2023-09-06 01:55:20
合計ジャッジ時間 5,884 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 567 ms
32,112 KB
testcase_01 AC 556 ms
32,116 KB
testcase_02 AC 552 ms
31,984 KB
testcase_03 AC 1 ms
4,380 KB
testcase_04 AC 20 ms
4,988 KB
testcase_05 AC 1 ms
4,380 KB
testcase_06 AC 571 ms
32,068 KB
testcase_07 AC 4 ms
4,380 KB
testcase_08 AC 3 ms
4,380 KB
testcase_09 AC 7 ms
5,208 KB
testcase_10 AC 2 ms
4,380 KB
testcase_11 AC 9 ms
6,060 KB
testcase_12 AC 11 ms
6,056 KB
testcase_13 AC 3 ms
4,376 KB
testcase_14 AC 419 ms
26,556 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cassert>// c
#include <iostream>// io
#include <iomanip>
#include <fstream>
#include <sstream>
#include <vector>// container
#include <map>
#include <set>
#include <queue>
#include <bitset>
#include <stack>
#include <algorithm>// other
#include <complex>
#include <numeric>
#include <functional>
#include <random>
#include <regex>
using namespace std;

using i32=int32_t;using i64=int64_t;using ll =i64;using uint=uint32_t;using ull=uint64_t;
template<typename T> using matrix=vector<vector<T> >;

#define ALL(c) (begin(c)),(end(c))
#define REP(i,n) FOR(i,0,n)
#define REPr(i,n) FORr(i,0,n)
#define FOR(i,l,r) for(int i=(int)(l);i<(int)(r);++i)
#define FORr(i,l,r) for(int i=(int)(r)-1;i>=(int)(l);--i)
#define EACH(it,o) for(auto it = (o).begin(); it != (o).end(); ++it)
#define IN(l,v,r) ((l)<=(v) && (v)<(r))
#define UNIQUE(v) v.erase(unique(ALL(v)),v.end())
//debug
#define DUMP(x)  cerr << #x << " = " << (x)
#define LINE()    cerr<< " (L" << __LINE__ << ")"

class range {
private:
	struct Iter{
		int v;
		int operator*(){return v;}
		bool operator!=(Iter& itr) {return v < itr.v;}
		void operator++() {++v;}
	};
	Iter i, n;
public:
	range(int n) : i({0}), n({n}) {}
	range(int i, int n) : i({i}), n({n}) {}
	Iter& begin() {return i;}
	Iter& end() {return n;}
};

//output
template<typename T> ostream& operator << (ostream& os, const vector<T>& as){REP(i,as.size()){if(i!=0)os<<" "; os<<as[i];}return os;}
template<typename T> ostream& operator << (ostream& os, const vector<vector<T>>& as){REP(i,as.size()){if(i!=0)os<<endl; os<<as[i];}return os;}
template<typename T> ostream& operator << (ostream& os, const set<T>& ss){for(auto a:ss){if(a!=ss.begin())os<<" "; os<<a;}return os;}
template<typename T1,typename T2> ostream& operator << (ostream& os, const pair<T1,T2>& p){os<<p.first<<" "<<p.second;return os;}
template<typename K,typename V> ostream& operator << (ostream& os, const map<K,V>& m){bool isF=true;for(auto& p:m){if(!isF)os<<endl;os<<p;isF=false;}return os;}
template<typename T1> ostream& operator << (ostream& os, const tuple<T1>& t){os << get<0>(t);return os;}
template<typename T1,typename T2> ostream& operator << (ostream& os, const tuple<T1,T2>& t){os << get<0>(t)<<" "<<get<1>(t);return os;}
template<typename T1,typename T2,typename T3> ostream& operator << (ostream& os, const tuple<T1,T2,T3>& t){os << get<0>(t)<<" "<<get<1>(t)<<" "<<get<2>(t);return os;}
template<typename T1,typename T2,typename T3,typename T4> ostream& operator << (ostream& os, const tuple<T1,T2,T3,T4>& t){os << get<0>(t)<<" "<<get<1>(t)<<" "<<get<2>(t)<<" "<<get<3>(t);return os;}
template<typename T1,typename T2,typename T3,typename T4,typename T5> ostream& operator << (ostream& os, const tuple<T1,T2,T3,T4,T5>& t){os << get<0>(t)<<" "<<get<1>(t)<<" "<<get<2>(t)<<" "<<get<3>(t)<<" "<<get<4>(t);return os;}
template<typename T1,typename T2,typename T3,typename T4,typename T5,typename T6> ostream& operator << (ostream& os, const tuple<T1,T2,T3,T4,T5,T6>& t){os << get<0>(t)<<" "<<get<1>(t)<<" "<<get<2>(t)<<" "<<get<3>(t)<<" "<<get<4>(t)<<" "<<get<5>(t);return os;}
template<typename T1,typename T2,typename T3,typename T4,typename T5,typename T6,typename T7> ostream& operator << (ostream& os, const tuple<T1,T2,T3,T4,T5,T6,T7>& t){os << get<0>(t)<<" "<<get<1>(t)<<" "<<get<2>(t)<<" "<<get<3>(t)<<" "<<get<4>(t)<<" "<<get<5>(t)<<" "<<get<6>(t);return os;}

//input
char tmp[1000];
#define nextInt(n) scanf("%d",&n)
#define nextLong(n) scanf("%lld",&n) //I64d
#define nextDouble(n) scanf("%lf",&n) 
#define nextChar(n) scanf("%c",&n)
#define nextString(n) scanf("%s",tmp);n=tmp

// values
template<typename T> T INF(){assert(false);};
template<> int INF<int>(){return 1<<28;};
template<> ll INF<ll>(){return 1LL<<58;};
template<> double INF<double>(){return 1e16;};

template<class T> T EPS(){assert(false);};
template<> int EPS<int>(){return 1;};
template<> ll EPS<ll>(){return 1LL;};
template<> double EPS<double>(){return 1e-8;};
template<> long double EPS<long double>(){return 1e-8;};

template<typename T,typename U> T pfmod(T v,U M){return fmod(fmod(v,M)+M,M);}

ll gcd_positive(ll a,ll b) { return b == 0 ? a : gcd_positive(b,a%b); }
ll gcd(ll a,ll b) { return gcd_positive(abs(a), abs(b)); }
ll lcm(ll a,ll b){return a/gcd(a,b)*b;}


class Main{
	public:

	void run(){
		int N,L;scanf("%d %d",&N,&L);
		L*=60;
		vector<int> ms(N),ss(N);
		REP(i,N)scanf("%d:%d",&ms[i],&ss[i]);
		vector<int> ts(N);
		REP(i,N)ts[i]=ms[i]*60+ss[i];

		// ループ
		int T=0;REP(i,N)T+=ts[i];
		if(T<=L){
			cout << N <<endl;return;
		}

		vector<double> fs(N+1);fs[0]=1;for(int i=1;i<=N;i++)fs[i]=i*fs[i-1];

		// DP = f_1...f_n
		// 条件: iの順序に関して可換 f_if_j = f_jf_i
		// (逆元はそのまま逆実行) f_i^(-1)
		// f_i^(-1)(DP) = f_1..f_{i-1}f_{i+1}..f_n

		vector<vector<double>> dp(N+1,vector<double>(2*L+1));
		dp[0][0]=1;
		REP(i,N){
			// f_i
			REPr(c,N+1)if(c-1>=0)REPr(l,2*L+1)if(l-ts[i]>=0)
				dp[c][l] += dp[c-1][l-ts[i]];
		}

		vector<vector<double>> tmp = dp;

		vector<double> ag(N+1);
		double all=0;
		REP(i,N){
			// f_i^(-1)
			REP(c,N+1)if(c-1>=0)REP(l,2*L+1)if(l-ts[i]>=0)
				dp[c][l] -= dp[c-1][l-ts[i]];

			REP(c,N+1)for(int l=0;l<L;l++)if(L<=l+ts[i]){
				ag[c+1] += fs[c] * dp[c][l] * fs[N-(c+1)];
				all += fs[c] * dp[c][l] * fs[N-(c+1)];
			}
			// f_i
			REPr(c,N+1)if(c-1>=0)REPr(l,2*L+1)if(l-ts[i]>=0)
				dp[c][l] += dp[c-1][l-ts[i]];

			REP(c,N+1)REP(l,2*L+1)assert(dp[c][l]==tmp[c][l]);
		}

		double res=0;
		REP(c,N+1){
			double v = c; double p = 1;
			p *= ag[c]/all;
			res+=v * p;
		}

		cout << res <<endl;
	}
};

int main(){
	cout <<fixed<<setprecision(20);
	cin.tie(0);
	ios::sync_with_stdio(false);
	Main().run();
	return 0;
}
0