#include <iostream>
#include <algorithm>
#include <iomanip>
#include <cstdio>
using namespace std;
#define MAX 18
double 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;
}
codeforce678E
最新推荐文章于 2020-02-05 16:14:57 发布