POJ 3071 Football

题目链接：http://poj.org/problem?id=3071

题意：有2^n个队伍，给出每个队伍i对战队伍j的胜率，求这些队伍中胜率最大的队伍编号。

思路：因为n<=7,所以最多7轮可以决出胜负，令f[i][j]表示第i轮队伍j获胜的概率，搞清楚第i次对局时，队伍j和哪些队伍可能会有对决，然后递推计算一下即可。
比如1~8：
第一回合（两个人） 1-2  3-4  5-6  7-8
第二回合（四个人一组，但是每组已经淘汰了两个人）如果1第一回合赢了，那么1可能和3或者4对局，其他同理。
第三回合（八个人一组，每组已经淘汰了6个人）
所以将一组分为两部分，前半部分只可能和后半部分的人对局，因为前半部分的人已经被打败淘汰了。

#include <cstdio>
#include <cmath>
#include <cstring>
#include <string>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <stack>
#include <map>
#include <set>
#include <vector>
#include <sstream>
#include <queue>
#include <utility>
using namespace std;

#define rep(i,j,k) for (int i=j;i<=k;i++)
#define Rrep(i,j,k) for (int i=j;i>=k;i--)

#define Clean(x,y) memset(x,y,sizeof(x))
#define LL long long
#define ULL unsigned long long
#define inf 0x7fffffff
#define mod %100000007
const int maxn = 1<<7;
int n,N;

double p[maxn+10][maxn+10];
double f[10][maxn+10];

bool init()
{
    scanf("%d",&n);
    if ( n == -1 ) return false;
    N = 1<<n;
    Clean(f,0);
    rep(i,1,N)
        rep(j,1,N) scanf("%lf",&p[i][j]);
    rep(i,1,N) f[0][i] = 1;
    return true;
}

void solve()
{
    rep(i,1,n) //第i个回合
    {
        int part = 1<<i;//多少人在一组内对决
        for( int j = 1; j <=N ; j+=part ) //一组一组的计算
        {
            int half = j + part/2 - 1;
            for( int k = j; k <= half; k++ )// 一组的前部分人
            {
                for( int kk = half + 1; kk < part + j; kk++ )//后

部分人
                {
                    f[i][k] += f[i-1][k]*f[i-1][kk]*p[k][kk];
                    f[i][kk] += f[i-1][kk]*f[i-1][k]*p[kk][k];
                }
            }
        }
    }
    double m = 0;
    int pos = 0;
    rep(i,1,N)
        if ( m < f[n][i] )
        {
            pos = i;
            m = f[n][i];
        }
    printf("%d\n",pos);
}

int main()
{
    while( init() )
    {
        solve();
    }
    return 0;
}