POJ-3071 Football解题报告（简单概率dp，水题）

题目描述

题目：https://vjudge.net/problem/POJ-3071#author=goodlife2017

思路分析

这是一道简单的求概率的题目。和算期望不同，算概率的时候是从开始向后推，推出结果的概率（求期望是从后开始向前推）。因此，这道题的思路就是从第一轮开始算。但是算的过程比较麻烦的一点：我们需要确定我们这一轮一共有多少种可能（有几支队伍可能和目前在算概率的这支队伍对上），因此，我们这里有（抄了）一个特殊的判别方法：

if(((k-1)>>i)==(j-1)>>i&&!((k-1)>>(i-1)==(j-1)>>(i-1)))

这里的意思就是说：i为轮数，而每轮过后，我们比较过的（包括前几轮比较过）的队伍数其实是2^i，比如，第1轮，对于1，我们比较过的有1、2两支队伍，第2轮，我们比较过的有1、2、3、4四支队伍……以此类推，因此，若这支队伍这一轮结束的时候需要和我们队进行比较，那么，就需要满足：(k-1)>>i==(j-1)，k为枚举需要比较的队伍，j为目前正在计算的队伍（这里-1是为了防止类似2、4这样的边缘的数不满足）（其实等式左右两边都是为0的）。另外，我们还需要防止之前已经比较过的再被比较一次，所以，需要再满足后面的条件：

!((k-1)>>(i-1)==(j-1)>>(i-1))

也就是说，这支队伍上一轮是没有被比较过的，从这一轮开始才需要比较的。
上面就是这道题的难点。

完整代码

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
const int maxn=1<<7+10;
double mp[200][200]={0};
double dp[8][200]={0};
int n;

int main()
{
	while(scanf("%d",&n)&&n!=-1)
	{
		memset(mp,0,sizeof(mp));
		memset(dp,0,sizeof(dp));
		for(int i=1;i<=(1<<n);i++)
		{
			for(int j=1;j<=(1<<n);j++)
			{
				scanf("%lf",&mp[i][j]);
			}
		}
		for(int i=1;i<=(1<<n);i++)
		{
			dp[0][i]=1.0;
		}
		for(int i=1;i<=n;i++)//轮数 
		{
			for(int j=1;j<=1<<n;j++)//队伍 
			{
				for(int k=1;k<=1<<n;k++)//枚举需要比较的队伍
				{
					if(((k-1)>>i)==(j-1)>>i&&!((k-1)>>(i-1)==(j-1)>>(i-1)))
					{
						dp[i][j]+=dp[i-1][j]*dp[i-1][k]*mp[j][k];
					}
				}
			}
		}
		int ans=1;
		for(int i=2;i<=1<<n;i++)
		{
			if(dp[n][ans]<dp[n][i])
			{
				ans=i;
			}
		}
		cout<<ans<<endl;
		
	}
	return 0;
 }