#include using namespace std; #define DEBUG(x) cerr<<#x<<": "< #define vl vector #define vii vector< vector > #define vll vector< vector > #define vs vector #define pii pair #define pis pair #define psi pair #define pll pair template pair operator+(const pair &s, const pair &t) { return pair(s.first + t.first, s.second + t.second); } template pair operator-(const pair &s, const pair &t) { return pair(s.first - t.first, s.second - t.second); } template ostream& operator<<(ostream& os, pair p) { os << "(" << p.first << ", " << p.second << ")"; return os; } #define X first #define Y second #define rep(i,n) for(int i=0;i<(int)(n);i++) #define rep1(i,n) for(int i=1;i<=(int)(n);i++) #define rrep(i,n) for(int i=(int)(n)-1;i>=0;i--) #define rrep1(i,n) for(int i=(int)(n);i>0;i--) #define REP(i,a,b) for(int i=a;i bool chmax(T &a, const T &b) { if (a bool chmin(T &a, const T &b) { if (a>b) { a = b; return 1; } return 0; } #define UNIQUE(v) v.erase(std::unique(v.begin(), v.end()), v.end()); const ll inf = 1000000001; const ll INF = (ll)1e18 + 1; const long double pi = 3.1415926535897932384626433832795028841971L; #define Sp(p) cout<>> memo(1 << N, vector>(N * 20 + 1, vector(2, -1))); double calc(int state, int score, int turn) { if (memo[state][score][turn] >= 0) return memo[state][score][turn]; int sum = 0; rep (i, n) { if (state & (1 << i)) sum += a[i]; } int you = sum - score; pair ma(-inf, -1); rep (i, n) { if (state & (1 << i)) continue; if (a[i] == 1) { double ex = 1 - calc(state | (1 << i), you, 1 - turn); chmax(ma, {ex, i}); continue; } double X = calc(state | (1 << i), you, 1 - turn); double mi = inf; rep (j, n) { if (state & (1 << j)) continue; double Y = calc(state | (1 << j), score, turn); if (a[j] == 1) { double f = 1. / a[i] * (1 - X) + (1 - 1. / a[i]) * (1 - (1 - Y)); chmin(mi, f); continue; } double r = 1. / a[i]; double f = r * (1 - X) + (1 - r) * (1 - 1. / a[j] * (1 - Y)) - (1 - r) * (1 - 1./a[j]); f /= (1 - (1. - r) * (1. - 1. / a[j])); chmin(mi, f); } chmax(ma, {mi, i}); } // DEBUG(state); // DEBUG(score); // DEBUG(turn); // DEBUG(ma); memo[state][score][turn] = ma.first; return ma.first; } signed main() { fio(); cin >> n; a.resize(n); int sum = 0; rep (i, n) { cin >> a[i]; sum += a[i]; } rep (i, sum + 1) { rep (j, 2) { int me = i, you = sum - i; if (me != you) { if (me > you) memo[(1 << n) - 1][i][j] = 1; else memo[(1 << n) - 1][i][j] = 0; } else { if (j == 0) memo[(1 << n) - 1][i][j] = 0; else memo[(1 << n) - 1][i][j] = 1; } } } Sp(calc(0, 0, 0)); }