JZOJ.4648 锦标赛

Problem

Description

403机房最近决定举行一场锦标赛。锦标赛共有N个人参加，共进行N-1轮。第一轮随机挑选两名选手进行决斗，胜者进入下一轮的比赛，第二轮到第N-1轮再每轮随机挑选1名选手与上一轮胜利的选手决斗，最后只剩一轮选手。第i名选手与第j名选手决斗，第i名选手胜利的概率是a[i][j].
作为一号选手的富榄想知道如何安排每轮出场的选手可以使得他获胜的概率最大，并求出这个最大概率。

Input

第一个一个整数N，表示参赛人数。
接下来N行，每行N个小数表示a[i][j]
保证满足a[i][i]=0，a[i][j]+a[j][i]=1；

Output

一个小数，表示最大概率。与答案误差不超过10^(-5) 即算正确。

Sample Input

3
0.0 0.5 0.8
0.5 0.0 0.4
0.2 0.6 0.0

Sample Output

0.6800000

Hint

第一轮2与3决斗，胜者与1决斗
如果2获胜，那么最后1获胜的概率为a[2][3]*a[1][2]=0.2
如果3获胜，那么最后1获胜的概率为a[3][2]*a[1][3]=0.48
所以总概率为0.68

Data Constraint

对于30%的数据，N<=10
对于60%的数据，N<=15
对于100%的数据，N<=18

Solution

我们看到N<=18，所以显然是状压DP。设F[s,i]为剩下s状态的人，（s为2进制数），当前擂主是i的最大概率。那么如果我们从前往后转，那么一个状态会变成两个状态，不好转，于是我们可以从后往前转，方程显然。

$$f[s][i]=max(f[s][i],f[s-2^j][i]*a[i][j]+f[s-2^i][j]*a[j][i])$$
初始值 $F_{1,1}=1$，最终答案为 $Max(F[全部人都在,i(i∈1-n)])$
最后不要忘了保留7位小数！！！

Code

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#define N 20
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
double f[262150][N],a[N][N],ans;
int _2[20],i,j,n,s;
int main()
{
    scanf("%d",&n);
    _2[1]=1;
    fo(i,2,19) _2[i]=_2[i-1]*2;
    fo(i,1,n) fo(j,1,n) scanf("%lf",&a[i][j]);
    f[1][1]=1;
    fo(s,0,_2[n+1]-1)
        fo(i,1,n)
            fo(j,1,n)
                if (j!=i && (s&_2[j]) && (s&_2[i]))
                    f[s][i]=max(f[s][i],f[s-_2[j]][i]*a[i][j]+f[s-_2[i]][j]*a[j][i]);
    fo(i,1,n) 
        if (f[_2[n+1]-1][i]>ans) ans=f[_2[n+1]-1][i];
    printf("%.7f",ans);
}

——2016.7.17