codeforce678E

#include <iostream>#include <algorithm>#include <iomanip>#include <cstdio>using namespace std;#define MAX 18double d[1<<MAX][MAX];//1<<MAX太大，不能再main函数里面声明//这里1<<MAX是指一个活着的人的集合，0为死，1为生，之后利用位运算进行相关操作//d[s][i]这个数组的值代表：当前状态为s（即现在的活人的集合）时主角夺冠的概率；//                          i先不考虑，最后会算出主角夺冠的n个概率，表示不同的n个人做初始擂主时主角夺冠概率//可得:d[s][i]=max(d[s-j][i]*p[i][j]+d[s-i][j]*p[j][i])//其中，s-j用s^(1<<j)这种位运算表示//比如：d[{1,2,3}][2]表示活人有123，擂主是2时，主角夺冠的概率//下一局，可能下一个选手比如1杀了2，变成d[{1,3}][1]表示活人有13，主角夺冠的概率//    也可能是2杀了下一个选手比如1，变成d[{2,3}][2]表示活人有23，主角夺冠的概率，俩种概率相加int main(){    int n;cin>>n;    double p[MAX][MAX];    for(int i=0;i<n;i++)        for(int j=0;j<n;j++)        cin>>p[i][j];    d[1][0]=1;//1的二进制为00...001（n-1个0）表示最后只有主角存活，从后往前推    for(int s=3;s<(1<<n);s+=2)        for(int i=0;i<n;i++)        if(s&(1<<i))//定位活人的序号    {        for(int j=0;j<n;j++)            if(s&(1<<j) && j!=i)            d[s][i]=max(d[s][i],p[i][j]*d[s^(1<<j)][i]+p[j][i]*d[s^(1<<i)][j]);    }    double ans;    for(int i=0;i<n;i++)        ans=max(ans,d[(1<<n)-1][i]);    printf("%.15lf",ans);    return 0;}