codeforce678E

#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;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值