1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
#include <bits/stdc++.h>
using namespace std;
const int MX=50005;
int n,m;
double a[MX],f[MX*2],best;
int main() {
  scanf("%d%d",&n,&m);
  for (int i=1; i<=n; i++) scanf("%lf",&a[i]);
  sort(a+1,a+n+1);
  reverse(a+1,a+n+1);
  f[MX]=1;
  for (int i=1; i<n; i++) {
    double cur=0;
    for (int j=(i&1); j<=i; j+=2) {
      if (j) f[MX-j]=f[MX-j-1]*a[i]+f[MX-j+1]*(1.-a[i]);
      f[MX+j]=f[MX+j-1]*a[i]+f[MX+j+1]*(1.-a[i]);
      if (j>=m) cur+=f[MX+j];
      // TODO: lft..rgh
    }
    best=max(best,cur);
  }
  printf("%.10lf\n",best);
  return 0;
}