Problem: 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 pi (0.0 ≤ pi ≤ 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.

Examples

Input

4
0.1 0.2 0.3 0.8


Output

0.800000000000


Input

2
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.1·0.8 + 0.9·0.2 = 0.26.

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
template <typename TA, typename TB>
ostream &operator<<(ostream &ostm, const pair<TA, TB> &p) {
  ostm << "(" << p.first << "," << p.second << ")";
  return ostm;
}
template <typename T>
ostream &operator<<(ostream &ostm, const vector<T> &v) {
  ostm << "[ ";
  for (auto i : v) ostm << i << " ";
  ostm << "]";
  return ostm;
}
template <typename TA, typename TB>
ostream &operator<<(ostream &ostm, const map<TA, TB> &mp) {
  ostm << "[ ";
  for (auto &it : mp) ostm << it << " ";
  ostm << "]";
  return ostm;
}
template <typename T>
ostream &operator<<(ostream &ostm, const set<T> &s) {
  ostm << "[ ";
  for (auto &it : s) ostm << it << " ";
  ostm << "]";
  return ostm;
}
template <typename T>
ostream &operator<<(ostream &ostm, const stack<T> &inp) {
  stack<T> st = inp;
  ostm << "[ ";
  while (!st.empty()) {
    ostm << st.top() << " ";
    st.pop();
  }
  ostm << "]";
  return ostm;
}
template <typename T>
ostream &operator<<(ostream &ostm, const queue<T> &inp) {
  queue<T> q = inp;
  ostm << "[ ";
  while (!q.empty()) {
    ostm << q.front() << " ";
    q.pop();
  }
  ostm << "]";
  return ostm;
}
template <typename T>
ostream &operator<<(ostream &ostm, const priority_queue<T> &inp) {
  priority_queue<T> pq = inp;
  ostm << "[ ";
  while (!pq.empty()) {
    ostm << pq.top() << " ";
    pq.pop();
  }
  ostm << "]";
  return ostm;
}
template <typename T>
ostream &operator<<(ostream &ostm, const deque<T> &inp) {
  deque<T> dq = inp;
  ostm << "[ ";
  while (!dq.empty()) {
    ostm << dq.front() << " ";
    dq.pop_front();
  }
  ostm << "]";
  return ostm;
}
inline int gtx() {
  const int N = 1048576;
  static char __buffer[N];
  static char *__p = __buffer, *__end = __buffer;
  if (__p == __end) {
    if ((__end = __buffer + fread(__buffer, 1, N, stdin)) == __buffer)
      return EOF;
    __p = __buffer;
  }
  return *__p++;
}
template <typename T>
inline bool rit(T &x) {
  char c = 0;
  bool fg = 0;
  while (c = gtx(), (c < '0' && c != '-') || c > '9')
    if (c == EOF) return false;
  c == '-' ? (fg = 1, x = 0) : (x = c - '0');
  while (c = gtx(), c >= '0' && c <= '9') x = x * 10 + c - '0';
  if (fg) x = -x;
  return true;
}
template <typename T, typename... Args>
inline bool rit(T &x, Args &...args) {
  return rit(x) && rit(args...);
}
inline void pit(int x) { printf("%d", x); }
inline void pln(long long x) { printf("%I64d", x); }
template <typename... Args>
inline void pit(int x, Args... args) {
  printf("%d ", x);
  pit(args...);
}
template <typename... Args>
inline void pln(long long x, Args... args) {
  printf("%I64d ", x);
  pit(args...);
}
void JIZZ() {
  cout << "";
  exit(0);
}
const long double PI = 3.14159265358979323846264338327950288;
const long double eps = 1e-8;
const long long mod = 1e9 + 7;
long double a[109];
int main() {
  int n;
  cin >> n;
  for (int i = 0; i < n; ++i) cin >> a[i];
  sort(a, a + n);
  reverse(a, a + n);
  long double p = a[0], allnt = 1 - a[0];
  for (int i = 1; i < n; ++i) {
    if (p < allnt * a[i] + p * (1 - a[i])) {
      p = allnt * a[i] + p * (1 - a[i]);
      allnt *= (1 - a[i]);
    }
  }
  cout << fixed << setprecision(10) << p << '\n';
}