結果

問題 No.155 生放送とBGM
ユーザー LayCurseLayCurse
提出日時 2015-01-23 19:51:08
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 560 ms / 6,000 ms
コード長 3,398 bytes
コンパイル時間 1,437 ms
コンパイル使用メモリ 148,128 KB
実行使用メモリ 52,372 KB
最終ジャッジ日時 2023-09-05 04:20:16
合計ジャッジ時間 6,156 ms
ジャッジサーバーID
(参考情報)
judge12 / judge13
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 560 ms
52,184 KB
testcase_01 AC 549 ms
52,184 KB
testcase_02 AC 560 ms
52,144 KB
testcase_03 AC 551 ms
52,372 KB
testcase_04 AC 41 ms
32,796 KB
testcase_05 AC 2 ms
5,772 KB
testcase_06 AC 556 ms
52,132 KB
testcase_07 AC 6 ms
10,536 KB
testcase_08 AC 2 ms
5,944 KB
testcase_09 AC 8 ms
10,628 KB
testcase_10 AC 3 ms
6,044 KB
testcase_11 AC 10 ms
10,876 KB
testcase_12 AC 11 ms
10,940 KB
testcase_13 AC 5 ms
10,444 KB
testcase_14 AC 431 ms
49,956 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include<bits/stdc++.h>
using namespace std;

#define REP(i,a,b) for(i=a;i<b;i++)
#define rep(i,n) REP(i,0,n)

#define mygc(c) (c)=getchar_unlocked()
#define mypc(c) putchar_unlocked(c)

#define ll long long
#define ull unsigned ll

#define double long double

void reader(int *x){int k,m=0;*x=0;for(;;){mygc(k);if(k=='-'){m=1;break;}if('0'<=k&&k<='9'){*x=k-'0';break;}}for(;;){mygc(k);if(k<'0'||k>'9')break;*x=(*x)*10+k-'0';}if(m)(*x)=-(*x);}
void reader(ll *x){int k,m=0;*x=0;for(;;){mygc(k);if(k=='-'){m=1;break;}if('0'<=k&&k<='9'){*x=k-'0';break;}}for(;;){mygc(k);if(k<'0'||k>'9')break;*x=(*x)*10+k-'0';}if(m)(*x)=-(*x);}
void reader(double *x){scanf("%Lf",x);}
int reader(char c[]){int i,s=0;for(;;){mygc(i);if(i!=' '&&i!='\n'&&i!='\r'&&i!='\t'&&i!=EOF) break;}c[s++]=i;for(;;){mygc(i);if(i==' '||i=='\n'||i=='\r'||i=='\t'||i==EOF) break;c[s++]=i;}c[s]='\0';return s;}
template <class T, class S> void reader(T *x, S *y){reader(x);reader(y);}
template <class T, class S, class U> void reader(T *x, S *y, U *z){reader(x);reader(y);reader(z);}
template <class T, class S, class U, class V> void reader(T *x, S *y, U *z, V *w){reader(x);reader(y);reader(z);reader(w);}

void writer(int x, char c){int s=0,m=0;char f[10];if(x<0)m=1,x=-x;while(x)f[s++]=x%10,x/=10;if(!s)f[s++]=0;if(m)mypc('-');while(s--)mypc(f[s]+'0');mypc(c);}
void writer(ll x, char c){int s=0,m=0;char f[20];if(x<0)m=1,x=-x;while(x)f[s++]=x%10,x/=10;if(!s)f[s++]=0;if(m)mypc('-');while(s--)mypc(f[s]+'0');mypc(c);}
void writer(double x, char c){printf("%.15Lf",x);mypc(c);}
void writer(const char c[]){int i;for(i=0;c[i]!='\0';i++)mypc(c[i]);}
void writer(const char x[], char c){int i;for(i=0;x[i]!='\0';i++)mypc(x[i]);mypc(c);}
template<class T> void writerLn(T x){writer(x,'\n');}
template<class T, class S> void writerLn(T x, S y){writer(x,' ');writer(y,'\n');}
template<class T, class S, class U> void writerLn(T x, S y, U z){writer(x,' ');writer(y,' ');writer(z,'\n');}
template<class T> void writerArr(T x[], int n){int i;if(!n){mypc('\n');return;}rep(i,n-1)writer(x[i],' ');writer(x[n-1],'\n');}

int N, L;
int S[120];

double brute(void){
  int i, j, k, ind[120], rem;
  double res = 0, all = 0;

  rep(i,N) ind[i] = i;
  do{
    rem = L;
    rep(i,N){
      res++;
      rem -= S[ind[i]];
      if(rem <= 0) break;
    }
    all++;
  }while(next_permutation(ind,ind+N));

  return res / all;
}


double dp[53][30000], ini[53][30000];
double solve(void){
  int i, j, k, m, p, mx = 61 * 60;
  double res = 0, pr, inv;

  rep(i,N+1) rep(j,L+mx) dp[i][j] = 0;
  dp[0][0] = 1;
  p = 0;
  rep(k,N){
    p++;
    for(i=p-1;i>=0;i--){
      pr = (double)(i+1) / (p+1);
      for(j=L+mx-1;j>=0;j--){
        dp[i+1][j+S[k]] += dp[i][j] * pr;
        dp[i][j] *= 1 - pr;
      }
    }
  }

  rep(m,N){
    rep(i,N+1) rep(j,L+mx) ini[i][j] = dp[i][j];
    rep(i,N){
      pr = (double)(i+1) / (N+1);
      inv = 1.0 / (1-pr);
      rep(j,L){
        ini[i][j] *= inv;
        ini[i+1][j+S[m]] -= ini[i][j] * pr;
      }
    }
    rep(i,N+1) rep(j,L) res += ini[i][j];
  }
  
  return res;
}

int main(){
  int i, j, k;
  double res;

  reader(&N,&L);
  assert(1 <= N && N <= 50);
  assert(1 <= L && L <= 300);
  L *= 60;
  rep(i,N){
    reader(&j,&k);
    S[i] = j * 60 + k;
    assert(0 <= j && j <= 59);
    assert(0 <= k && k <= 59);
  }

//  res = brute();
//  writerLn(res);
  res = solve();
  writerLn(res);

  return 0;
}
0