結果

問題 No.155 生放送とBGM
ユーザー LayCurseLayCurse
提出日時 2015-01-23 19:35:44
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 564 ms / 6,000 ms
コード長 3,800 bytes
コンパイル時間 1,514 ms
コンパイル使用メモリ 162,820 KB
実行使用メモリ 23,064 KB
最終ジャッジ日時 2024-06-23 00:57:59
合計ジャッジ時間 5,751 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 544 ms
22,772 KB
testcase_01 AC 530 ms
23,064 KB
testcase_02 AC 516 ms
22,736 KB
testcase_03 AC 469 ms
22,892 KB
testcase_04 AC 17 ms
18,008 KB
testcase_05 AC 2 ms
6,940 KB
testcase_06 AC 564 ms
22,848 KB
testcase_07 AC 3 ms
6,944 KB
testcase_08 AC 2 ms
6,944 KB
testcase_09 AC 4 ms
8,148 KB
testcase_10 AC 2 ms
6,948 KB
testcase_11 AC 4 ms
9,996 KB
testcase_12 AC 4 ms
6,944 KB
testcase_13 AC 2 ms
6,940 KB
testcase_14 AC 409 ms
21,968 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘void reader(double*)’:
main.cpp:15:29: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   15 | void reader(double *x){scanf("%lf",x);}
      |                        ~~~~~^~~~~~~~~

ソースコード

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

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("%.15f",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][24000], ini[53][24000];
double solve(void){
  int i, j, k, m, p, loop, block;
  int st, ed, sum, ini_p, ini_sum;
  double res = 0, pr;

  block = sqrt(N);
  rep(loop,block){
    st = loop * N / block;
    ed = (loop + 1) * N / block;
    if(st==ed) continue;

    p = 0; sum = 0;
    rep(i,N+1) rep(j,L) ini[i][j] = 0;
    ini[0][0] = 1;
    rep(k,N) if(k<st || k>=ed){
      p++;
      for(i=p-1;i>=0;i--){
        pr = (double)(i+1) / (p+1);
        for(j=sum;j>=0;j--){
          ini[i+1][j+S[k]] += ini[i][j] * pr;
          ini[i][j] *= 1 - pr;
        }
      }
      sum = min(sum+S[k], L-1);
    }
    ini_p = p;
    ini_sum = sum;
    
    REP(m,st,ed){
      rep(i,N+1) rep(j,L) dp[i][j] = ini[i][j];
      p = ini_p;
      sum = ini_sum;
      REP(k,st,ed) if(k!=m){
        p++;
        for(i=p-1;i>=0;i--){
          pr = (double)(i+1) / (p+1);
          for(j=sum;j>=0;j--){
            dp[i+1][j+S[k]] += dp[i][j] * pr;
            dp[i][j] *= 1 - pr;
          }
        }
        sum = min(sum+S[k], L-1);
      }
      rep(i,N+1) rep(j,L) res += dp[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