結果

問題 No.155 生放送とBGM
コンテスト
ユーザー mamekin
提出日時 2015-02-18 23:41:25
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 1,106 ms / 6,000 ms
コード長 1,837 bytes
コンパイル時間 1,524 ms
コンパイル使用メモリ 108,220 KB
実行使用メモリ 362,336 KB
最終ジャッジ日時 2024-06-23 21:14:05
合計ジャッジ時間 10,078 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 15
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <cstdio>
#include <iostream>
#include <sstream>
#include <fstream>
#include <iomanip>
#include <algorithm>
#include <cmath>
#include <string>
#include <vector>
#include <list>
#include <queue>
#include <stack>
#include <set>
#include <map>
#include <bitset>
#include <numeric>
#include <limits>
#include <climits>
#include <cfloat>
#include <functional>
using namespace std;

void solve(const vector<int>& s, int len, vector<vector<vector<double> > >& dp)
{
    int n = s.size();
    dp.assign(1, vector<vector<double> >(1, vector<double>(len, 0.0)));
    dp[0][0][0] = 1.0;
    for(int i=0; i<n-1; ++i){
        dp.push_back(vector<vector<double> >(i+2, vector<double>(len, 0.0)));
        for(int a=0; a<=i; ++a){
            for(int b=0; b<len; ++b){
                dp[i+1][a][b] += dp[i][a][b];
                if(b + s[i] < len)
                    dp[i+1][a+1][b+s[i]] += dp[i][a][b];
            }
        }
    }
}

int main()
{
    int n, len;
    cin >> n >> len;
    len *= 60;

    vector<int> s(n);
    for(int i=0; i<n; ++i){
        int a, b;
        char c;
        cin >> a >> c >> b;
        s[i] = a * 60 + b;
    }

    vector<double> power(n, 1);
    for(int i=1; i<n; ++i)
        power[i] = power[i-1] * i;

    vector<vector<vector<double> > > dp1, dp2;
    solve(s, len, dp1);
    reverse(s.begin(), s.end());
    solve(s, len, dp2);

    double ret = 0.0;
    for(int i=0; i<n; ++i){
        for(int a=0; a<=i; ++a){
            for(int b=0; b<n-i; ++b){
                double tmp = 0.0;
                for(int j=len-1; j>=0; --j){
                    tmp += dp2[n-1-i][b][len-1-j];
                    ret += dp1[i][a][j] * tmp * power[a+b] * power[n-1-a-b];
                }
            }
        }
    }

    for(int i=1; i<=n; ++i)
        ret /= i;
    printf("%.10f\n", ret);

    return 0;
}
0