Codeforces Beta Round #16 E. Fish (状压dp)（概率dp）

Codeforces Beta Round #16 (Div. 2 Only) E. Fish

题目链接：## 点击打开链接

题意：##

有 (n) 条鱼，每两条鱼相遇都会有其中一只吃掉对方，现在给你一个 (n * n)的矩阵，表示 (i) 吃掉 (j) 的概率，最后问你每条鱼存活的概率。

题解：##

最多有 (18) 条鱼，吃掉的概率都不一样，可以用状态压缩，设(dp[1<<n])种状态，最多有 (1<<18) 种状态。
$ dp[i]$ 表示当前鱼的状态为 (i) 时的概率。
那么，(dp( i 吃掉 j ) = dp( i 和 j 同时存在) * p( i 战胜 j )的概率 * prob(i 和 j 相遇)的概率)。

f[i - (1<<k)] += f[i] * prob[j][k] * 1.0 / (cnt * (cnt - 1) / 2) ;// k 被吃掉 
f[i - (1<<j)] += f[i] * prob[k][j] * 1.0 / (cnt * (cnt - 1) / 2); // j 被吃掉 

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
double prob[18][18];
double f[1<<18];
int cnt;
int main()
{
    int n;
    cin>>n;
    for(int i=0;i<n;i++)
    {
    	for(int j=0;j<n;j++)
    	{
    		cin>>prob[i][j];
		}
	}
	f[(1<<n)-1] = 1;//所有鱼都在，初始化概率都是 1 
	for(int i=(1<<n)-1;i>0;--i)//遍历所有状态 
	{
		cnt = 0;
		for(int j=0;j<n;j++) //遍历获胜的鱼的编号
		{
			if(i & (1<<j))//j号鱼还活着
			{
				cnt++;
			}
		}
		for(int j=0;j<n;j++)//遍历鱼j要吃掉的k号鱼  
		{
			if((i&(1<<j))==0)continue;//j号鱼不存活 
			
			for(int k=j+1;k<n;k++) 
			{
				if((i&(1<<k))==0)continue;//k号鱼不存活 
				//f( j 吃掉 k ) = f( j 和 k 同时存在) * f( j 战胜k) * f( j 和 k 相遇)
				//1<<k : 第k位的1一定存在 
				f[i - (1<<k)] += f[i] * prob[j][k] * 1.0 / (cnt * (cnt - 1) / 2);//k被吃掉 
				f[i - (1<<j)] += f[i] * prob[k][j] * 1.0 / (cnt * (cnt - 1) / 2);//j被吃掉 
			}
		}
	}
	
	for(int i=0;i<n-1;i++)
	{
		printf("%.6f ",f[1<<i]);//分别输出只有第 i 条鱼存活的概率
	}
	printf("%.6f
",f[1<<(n-1)]);
    return 0;
}