#include<stdio.h>
#include<string.h>
const int N=18;
double dp[1<<N][N],p[N][N];
double mmax(double x,double y)
{
	return x>y?x:y;
}
void solve()
{
	int n;
	memset(dp,0,sizeof(dp));
	scanf("%d",&n);
	for(int i=0;i<n;++i)
	for(int j=0;j<n;++j)
	scanf("%lf",&p[i][j]);
	dp[1][0]=1.0;
	for(int i=0;i<(1<<n);++i)
	for(int j=0;j<n;++j)
	if((i&(1<<j))) 
	for(int k=0;k<n;++k)
	if((i&(1<<k))&&j!=k) 
	dp[i][j]=mmax(dp[i][j],dp[i^(1<<j)][k]*p[k][j]+dp[i^(1<<k)][j]*p[j][k]);
	double ans=0.0;
	for(int i=0;i<n;++i)
	if(dp[(1<<n)-1][i]>ans)ans=dp[(1<<n)-1][i];
	printf("%.8f\n",ans);
}
int main()
{
	solve();
	return 0;
}