poj3071 概率dp

由题可知，本题为概率dp。d[i][j] 表示第i个人从第j轮晋级的概率。
则dp[i][j] = ∑〖dp[k][j-1] * p[i][k]〗* dp[i][j-1],其中k为在当前轮次所有可能与i交手的对手。上式表示本轮中于每一个可能的对手的出现概率与胜率乘积再乘上本人晋级本轮的概率。最后答案只要找到最大值即可。
代码如下：

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cmath>
using namespace std;
const int maxn = 130;
double dp[maxn][8];
double rate[maxn][maxn];
int n,m;
int main()
{
    while(~scanf("%d",&m) && m > 0){
        n = (1 << m);
        for(int i = 1;i <= n;i++){
            for(int j = 1;j <= n;j++){
                scanf("%lf",&rate[i][j]);
            }
            dp[i][0] = 1;
        }
        for(int j = 1;j <= m;j++){
            int p = (1 << (j - 1)),q = (1 << j);
            for(int i = 1;i <= n;i++){
                double tmp = 0;
                for(int k = 1;k <= n;k++){
                    if(i == k) continue;
                    if((i - 1) / p != (k - 1) / p && (i - 1) / q == (k - 1) / q){
                        tmp += dp[k][j-1] * rate[i][k];

                    }       
                }
                dp[i][j] = dp[i][j-1] * tmp;
            }
        }
        int ans = 0;
        double maxv  = 0;
        for(int i = 1;i <= n;i++){
            if(dp[i][m] > maxv){maxv = dp[i][m];ans = i;}
        }
        printf("%d\n",ans);
    }
    return 0;
}