cf 442B Andrey and Problem
B. Andrey and Problem
time limit per test
2 seconds
memory limit per test
256 megabytes
input
standard input
output
standard output
Andrey needs one more problem to conduct a programming contest. He has n friends who are always willing to help. He can ask some of them to come up with a contest problem. Andrey knows one value for each of his fiends -- the probability that this friend will come up with a problem if Andrey asks him.
Help Andrey choose people to ask. As he needs only one problem, Andrey is going to be really upset if no one comes up with a problem or if he gets more than one problem from his friends. You need to choose such a set of people that maximizes the chances of Andrey not getting upset.
Input
The first line contains a single integer n (1 ≤ n ≤ 100) -- the number of Andrey's friends. The second line contains n real numbers p__i(0.0 ≤ p__i ≤ 1.0) -- the probability that the i-th friend can come up with a problem. The probabilities are given with at most 6 digits after decimal point.
Output
Print a single real number -- the probability that Andrey won't get upset at the optimal choice of friends. The answer will be considered valid if it differs from the correct one by at most 10 - 9.
Sample test(s)
input
4<br></br>0.1 0.2 0.3 0.8
output
0.800000000000
input
2<br></br>0.1 0.2
output
0.260000000000
Note
In the first sample the best strategy for Andrey is to ask only one of his friends, the most reliable one.
In the second sample the best strategy for Andrey is to ask all of his friends to come up with a problem. Then the probability that he will get exactly one problem is 0.10.8 + 0.90.2 = 0.26.
可以选1个,可以选2个,选i个....
但是选哪i个,似乎还要枚举...这样计算起来很麻烦
但是再一看,很容易发现,选i个概率最大的情况一定是选本身概率最大的i个的情况
那么我们只需要枚举选的个数,取每种情况的最大值就好了.
1
2
3 /*************************************************************************
4 > File Name: code/2015summer/#3/C.cpp
5 > Author: 111qqz
6 > Email: rkz2013@126.com
7 > Created Time: 2015年07月28日 星期二 12时27分50秒
8 ************************************************************************/
9
10 #include<iostream>
11 #include<iomanip>
12 #include<cstdio>
13 #include<algorithm>
14 #include<cmath>
15 #include<cstring>
16 #include<string>
17 #include<map>
18 #include<set>
19 #include<queue>
20 #include<vector>
21 #include<stack>
22 #define y0 abc111qqz
23 #define y1 hust111qqz
24 #define yn hez111qqz
25 #define j1 cute111qqz
26 #define tm crazy111qqz
27 #define lr dying111qqz
28 using namespace std;
29 #define REP(i, n) for (int i=0;i<int(n);++i)
30 typedef long long LL;
31 typedef unsigned long long ULL;
32 const int N=1E2+7;
33 double p[N];
34 int n;
35
36 bool cmp(double a,double b)
37 {
38 if (a>b) return true;
39 return false;
40 }
41 double solve (int x)
42 {
43 double res=0;
44 double poss;
45 for ( int i = 1 ; i <= x; i++)
46 {
47 poss = 1.0;
48 for ( int j = 1 ; j <= x; j++ )
49 {
50 if (j==i) poss=poss*p[j-1];
51 else poss = poss * (1-p[j-1]);
52 }
53 res = res + poss;
54 }
55 return res;
56 }
57 int main()
58 {
59 cin>>n;
60 for ( int i = 0 ; i < n ; i++ )
61 {
62 cin>>p[i];
63 }
64 sort(p,p+n,cmp);
65 double ans = 0;
66 for ( int i = 0 ; i < n ; i++ )
67 {
68 ans = max(ans,solve(i+1));
69 }
70 cout<<fixed<<setprecision(12)<<ans<<endl;
71
72 return 0;
73 }
74
75