#include <cstdio>

#include <iostream>

#include <ctime>

#include <vector>

#include <cmath>

#include <map>

#include <queue>

#include <algorithm>

#include <cstring>

using namespace std;

typedef long long LL;

const int N=5e2+;

const int INF=0x3f3f3f3f;

const int mod=1e9+;

double p[][],dp[<<][];

int n;

double gao(int s,int cur){

  if(dp[s][cur]!=-1.0) return dp[s][cur];

  for(int i=;i<n;++i){

     if(i!=cur&&(s&(<<i))){

       dp[s][cur]=max(dp[s][cur],p[cur][i]*gao(s^(<<i),cur)+p[i][cur]*gao(s^(<<cur),i));

     }

  }

  return dp[s][cur];

}

int main(){

   scanf("%d",&n);

   for(int i=;i<n;++i){

    for(int j=;j<n;++j)

      scanf("%lf",&p[i][j]);

   }

   for(int i=;i<(<<n);++i)

    for(int j=;j<;++j)

       dp[i][j]=-;

   for(int i=;i<n;++i)

    dp[(<<i)][i]=;

   dp[][]=;

   double ans=;

   for(int i=;i<n;++i){

     ans=max(ans,gao((<<n)-,i));

   }

   printf("%.10f\n",ans);

   return ;

}