题目:
403机房最近决定举行一场锦标赛。锦标赛共有N个人参加,共进行N-1轮。第一轮随机挑选两名选手进行决斗,胜者进入下一轮的比赛,第二轮到第N-1轮再每轮随机挑选1名选手与上一轮胜利的选手决斗,最后只剩一轮选手。第i名选手与第j名选手决斗,第i名选手胜利的概率是a[i][j].
作为一号选手的富榄想知道如何安排每轮出场的选手可以使得他获胜的概率最大,并求出这个最大概率。
样例输入:
第一个一个整数N,表示参赛人数。
接下来N行,每行N个小数表示a[i][j]
保证满足a[i][i]=0,a[i][j]+a[j][i]=1;
3
0.0 0.5 0.8
0.5 0.0 0.4
0.2 0.6 0.0
样例输出:
一个小数,表示最大概率。与答案误差不超过10^(-5) 即算正确。
0.6800000
数据范围:
对于30%的数据,N<=10
对于60%的数据,N<=15
对于100%的数据,N<=18
剖解题目:
题目简洁明了,我的语文水平已经无法剖解了。。。
思路:
看到数据范围,N<=18,嗯,顷刻间想到了状压DP,然后看看题目,发现可以往下想。
解法:
确实是状压dp,但是在推的时候推着推着就跪了,因为正推发现下一个状态有win与lost两种情况,然后推出状态会增多,然后就炸了!
事实上,思路是正确的,只需要倒过来推便可以了。
设f[i][s]表示,在第i个人是winner时,剩下人状态为s的最大获胜概率,其中s为一个二进制数,1位表示这个人还没死透,0表示这个人已经死透了。
那么,倒过来推就有:
fi,s=max(fi,s−2j∗ai,j+fj,s−2i∗aj,i)
其中有s and 2j ==1 且s and 2i ==1且i!=j.
时间复杂度: O(2n∗n2) .
代码:
#include<cstdio>
#include<cstdlib>
#include<cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int maxn=19;
int n;
double a[maxn][maxn],f[maxn][1<<maxn],ans;
double max(double x,double y)
{
if (x>y) return x;
else return y;
}
int main()
{
//freopen("T2.in","r",stdin);
scanf("%d",&n);
fo(i,1,n)
fo(j,1,n) scanf("%lf",&a[i][j]);
f[1][1]=1;
fo(s,1,(1<<n)-1){
fo(i,1,n)
if (s&(1<<i-1))
fo(j,1,n)
if ((j!=i) && (s&(1<<j-1)))
f[i][s]=max(f[i][s],f[i][s-(1<<j-1)]*a[i][j]+f[j][s-(1<<i-1)]*a[j][i]);
}
fo(i,1,n) ans=max(ans,f[i][(1<<n)-1]);
printf("%.7lf",ans);
//fclose(stdin);
}