forked from ngthanhtrung23/CompetitiveProgramming
-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathD.cpp
More file actions
63 lines (45 loc) · 1.11 KB
/
D.cpp
File metadata and controls
63 lines (45 loc) · 1.11 KB
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
#include<bits/stdc++.h>
using namespace std;
#define double long double
int n, x;
double P[105];
double Pmiss[105][100005];
double Pmcur, Phcur;
int guess[105];
double E;
int main() {
cin >> n;
for(int i=1; i<=n; i++) {
cin >> x;
P[i] = x/100.0;
}
for(int i=1; i<=n; i++) {
Pmiss[i][1] = 1 - P[i];
for(int j=2; j<=100000; j+=2) {
Pmiss[i][j] = Pmiss[i][j/2] * Pmiss[i][j/2];
Pmiss[i][j+1] = Pmiss[i][j] * Pmiss[i][1];
}
}
E = n * 1.0;
Phcur = 1.0;
for(int i=1; i<=n; i++) {
guess[i] = 1;
Phcur *= P[i];
}
int cur = 0;
while(true) {
E += (1 - Phcur); //Pmisscur
int next = 1;
for(int i=2; i<=n; i++) {
if ((1-Pmiss[i][guess[i]+1]) / (1-Pmiss[i][guess[i]]) > (1-Pmiss[next][guess[next]+1])/(1-Pmiss[next][guess[next]])) next = i;
}
guess[next]++;
if (next >= 100000) break;
cur++;
if (cur > n * 10000) break;
Phcur = 1.0;
for(int i=1; i<=n; i++) Phcur *= (1-Pmiss[i][guess[i]]);
}
for(int i=1; i<=n; i++) E += Pmiss[i][guess[i]] / P[i];
cout << (fixed) << setprecision(9) << E;
}