1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
#include <bits/stdc++.h>
#define debug(x) cout<<#x<<" = "<<x<<"\n"
#define pb push_back
#define ins insert
#define fi first
#define se second
 
using namespace std;
using ull = unsigned long long;
using ll = long long;
using ld = double; //zmiana!
 
template <typename H, typename T> 
ostream& operator<<(ostream& os, pair<H, T> m){
	return os <<"("<< m.fi<<", "<<m.se<<")";
}
 
template <typename H> 
ostream& operator<<(ostream& os, vector<H> V){
	os<<"{";
	for(int i=0; i<V.size(); i++){
		if(i)os<<" ";
		os<<V[i];
	}
	os<<"}";
	return os;
}
 
void solve();
  
int main(){
	ios_base::sync_with_stdio(false);
	cin.tie(NULL);
	cout<<fixed<<setprecision(10);
	int t = 1;
	//cin>>t;
	while(t--) solve();
	return 0;
}
  
const int OFFSET = 1800;
const int SRODEK = 20;
ld rozklad[OFFSET+1][2];

 
void solve(){
	//pamietaj o formatowaniu!!!!!!!
	int n,t;
	int cnt = 0;
	ld wynik = 0;
	cin>>n>>t;
	ld odp = 0.0;
	
	for(int i = 0; i < OFFSET+1; i++) rozklad[i][0] = 0.0;
	rozklad[SRODEK][0] = 1.0;
	vector<ld> vec;
	for(int i = 0; i < n; i++){
		ld x; cin>>x; vec.pb(x);
	}
	sort(vec.begin(),vec.end(),greater<ld>());
	bool  pewna = false;
	
	for(int i = 0; i < vec.size(); i++){
		ld p = vec[i];
		if(p >= 0.5) pewna = true; else pewna = false;
		if(pewna){
			wynik+=1;
			for(int j  = 1;j < OFFSET; j++){ // zostawiam zera na bokach intencjonalnie
				rozklad[j][cnt^1] = rozklad[j][cnt]*p +rozklad[j+1][cnt] *(1.0-p);

			}



		}else{
			wynik-=1;
			for(int j  = 1;j < OFFSET; j++){ // zostawiam zera na bokach intencjonalnie
				rozklad[j][cnt^1] = rozklad[j][cnt]*(1.0-p) +rozklad[j-1][cnt] *(p);
			}



		}
		ld sum = 0.0;
		ld total = 0.0;
		ld total2 = 0.0;
		int mini = -1;
		int maxi = -1;
		for(int j  = 1;j < OFFSET; j++){ // zostawiam zera na bokach intencjonalnie
			if((j-SRODEK)*2+wynik >= t)sum += rozklad[j][cnt^1];
			total2+=rozklad[j][cnt^1];

		}
		for(int j  = 1;j < OFFSET; j++){ // zostawiam zera na bokach intencjonalnie
			total+=rozklad[j][cnt^1];
			if(total > 0.000000000001*total2 && mini == -1) mini = j;  
			if(total > (1-0.000000000001)*total2 && maxi == -1) maxi = j;  
		}

		odp = max(odp,sum);

	//	debug(maxi-mini);
	//	debug(mini);
		int diff = SRODEK - mini;
		wynik -= diff*2;
		for(int j = 1; j <  OFFSET; j++){
			
			if(j - diff > 0 && j-diff < OFFSET) rozklad[j][cnt] = rozklad[j-diff][cnt^1];
			else rozklad[j][cnt] = 0.0;

		}
	//	debug(p);
	//	debug(1.0-total);
	//	debug(rozklad[2][0]);
	//	debug(rozklad[2*OFFSET-1][0]);
	//	if(total < 0.5) assert(0==2);
		//debug(rozklad[OFFSET+50][cnt^1]);



		//cnt^=1;
	//	if(i +1 == vec.size()) debug(total);
	}
	cout<<odp<<"\n";

	
	


}