#include using namespace std; typedef long long ll; #define pb push_back #define fi first #define se second typedef pair P; using VP = vector

; using VVP = vector; using VI = vector; using VVI = vector; using VVVI = vector; const int inf = 1e9 + 7; const ll INF = 1LL << 61; const ll mod = 1e9 + 7; template inline bool chmax(T &a, T b) { if (a < b) { a = b; return true; } return false; } template inline bool chmin(T &a, T b) { if (a > b) { a = b; return true; } return false; } int n; VI makevec(int x) { VI v(n); for (int j = 0; j < n; j++) { if (x % 3 == 0) v[j] = 0; else if (x % 3 == 1) v[j] = 1; else v[j] = 2; x /= 3; } return v; } int rev(int x) { int j; int y = 0; int k = 1; for (j = 0; j < n; j++) { if (x % 3 == 1) y += k * 2; else if (x % 3 == 2) y += k; x /= 3; k *= 3; } return y; } int main() { ios::sync_with_stdio(false); cin.tie(nullptr); int i, j, k; cin >> n; int a[n]; for (i = 0; i < n; i++) cin >> a[i]; int all = 0; for (i = 0; i < n; i++) all += a[i]; int m = pow(3, n); VI fl(m, -1); for (i = 0; i < m; i++) { VI v = makevec(i); int p = 0, q = 0; bool fin = 1; for (j = 0; j < n; j++) { if (v[j] == 1) p += a[j]; else if (v[j] == 2) q += a[j]; else fin = 0; } if (p > all - p) fl[i] = 1; else if (q >= all - q) fl[i] = 2; } double get[n]; double iget[n]; for (i = 0; i < n; i++) { get[i] = double(1.0) / a[i]; iget[i] = 1.0 - get[i]; } vector dp(m, 0.0); for (i = m - 1; i >= 0; i--) { if (fl[i] == 1) { dp[i] = 1.0; //cout << i << " " << dp[i] << endl; continue; } else if (fl[i] == 2) { dp[i] = 0.0; // cout << i << " " << dp[i] << endl; continue; } VI v = makevec(i); for (j = 0; j < n; j++) { if (v[j] != 0) continue; // 取れた場合 double w = 1.0; double w2 = 1.0; if (dp[i + pow(3, j)] < 0.00000001) w = 0.0; for (k = 0; k < n; k++) { if (j == k) continue; if (v[k] != 0) continue; chmin(w, get[k] * dp[i + pow(3, j) + 2 * pow(3, k)] + iget[k] * dp[i + pow(3, j)]); } // 取れない bool f = 0; for (k = 0; k < n; k++) { if (v[k] != 0) continue; f = 1; double cont = iget[j] * iget[k]; chmin(w2, (get[j] * w + iget[j] * get[k] * dp[i + 2 * pow(3, k)]) / (1.0 - cont)); } if (f) chmax(dp[i], w2); } //cout << i << " " << dp[i] << endl; } cout << fixed << setprecision(15) << dp[0] << endl; }