#include <bits/stdc++.h>
using namespace std;
long double pr[18][18];
long double dp[(1 << 18)][18];
int n;
long double rec(int mask, int last) {
  if (dp[mask][last]) {
    return dp[mask][last] - 1;
  }
  if (mask == (1 << n) - 1 && last == 0) {
    dp[mask][last] = 2;
    return 1;
  } else if (mask == (1 << n) - 1 && last != 0) {
    dp[mask][last] = 1;
    return 0;
  }
  int i;
  long double mx = -1;
  for (i = 0; i < n; i++) {
    if (!((mask >> i) & 1)) {
      mx = max(mx, pr[last][i] * rec(mask | (1 << i), last) +
                       pr[i][last] * rec(mask | (1 << i), i));
    }
  }
  dp[mask][last] = mx + 1;
  return mx;
}
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(NULL);
  int i, j;
  long double mx = -1;
  cin >> n;
  if (n == 1) {
    cout << "1";
    return 0;
  }
  for (i = 0; i < n; i++) {
    for (j = 0; j < n; j++) {
      cin >> pr[i][j];
    }
  }
  for (i = 0; i < n; i++) {
    for (j = i + 1; j < n; j++) {
      mx = max(mx, rec((1 << i) | (1 << j), i) * pr[i][j] +
                       rec((1 << i) | (1 << j), j) * pr[j][i]);
    }
  }
  cout << setprecision(15) << mx;
  return 0;
}
