poj 3071 简单概率dp

题意很清晰，就是问最有可能获得冠军的队伍。

需要注意的是每轮比赛中，每个队都是和自己旁边的一个队比赛，采用淘汰赛制，所以需要决定第i轮的时候j队可以和哪些队比赛，然后求概率dp即可。

状态转移方程：

　　dp[i][j] += dp[i - 1][j] * dp[i - 1][k] * p[j][k];

　　dp[i][j]表示第i轮j队获胜的概率 = 第i-1轮j队获胜的概率 * 第i-1轮k队获胜的概率 * j队打败k队的概率（之和）。

　　前提是j和k在第i轮可能遇到。

　　判断方法： ( j >> ( i - 1 ) ) ^ 1 == k >> ( i - 1 ) 

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;

const int N = 8;
const int M = 1 << 7;
double dp[N][M];
double p[M][M];

int main ()
{
    int n, m;
    while ( scanf("%d", &n) != EOF )
    {
        if ( n == -1 ) break;
        m = 1 << n;
        for ( int i = 0; i < m; i++ )
        {
            for ( int j = 0; j < m; j++ )
            {
                scanf("%lf", &p[i][j]);
            }
        }
        for ( int j = 0; j < m; j++ )
        {
            dp[0][j] = 1.0;
        }
        for ( int i = 1; i <= n; i++ )
        {
            for ( int j = 0; j < m; j++ )
            {
                dp[i][j] = 0;
                for ( int k = 0; k < m; k++ )
                {
                    if ( ( ( j >> ( i - 1 ) ) ^ 1 ) == ( k >> ( i - 1 ) ) )
                    {
                        dp[i][j] += dp[i - 1][j] * dp[i - 1][k] * p[j][k];
                    }
                }
            }
        }
        int ans = 0;
        for ( int j = 1; j < m; j++ )
        {
            if ( dp[n][j] > dp[n][ans] )
            {
                ans = j;
            }
        }
        printf("%d\n", ans + 1);
    }
    return 0;
}