poj 3071 Football

http://poj.org/problem?id=3071

 记忆化搜索

代码及其注释：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<map>
#include<cmath>
#define LL long long

using namespace std;

const int N=130;

double win[N][N];//第几队 在第几次比赛中获胜的概率
double beat[N][N];// i赢j的概率
double pro(int k,int x,int l,int r)//第k队 第x次比赛 l和r用来确认范围
{
    if(win[k][x]>=0.0)
    return win[k][x];//记忆化
    if(x==0)
    {
        win[k][x]=1.0;//第0次一定获胜 为边界
        return win[k][x];
    }
    int mid=(l+r)>>1;
    double ftemp=0.0;
    if(k<=mid)//判定k 所在的范围 和他的对手范围
    {
        for(int i=mid+1;i<=r;++i)
        {
            ftemp+=(beat[k][i]*pro(i,x-1,mid+1,r));//累加  对手出现的概率和打败概率的乘积
        }
        win[k][x]=pro(k,x-1,l,mid)*ftemp;//最后乘上自己出现的概率
    }else
    {
        for(int i=l;i<=mid;++i)
        {
            ftemp+=(beat[k][i]*pro(i,x-1,l,mid));
        }
        win[k][x]=pro(k,x-1,mid+1,r)*ftemp;
    }
    return win[k][x];
}
int main()
{
    int n;
    while(scanf("%d",&n)!=EOF)
    {
        if(n==-1)
        break;
        int m=int(pow(2,n));
        for(int i=1;i<=m;++i)
        {
            for(int j=1;j<=m;++j)
            {
                scanf("%lf",&beat[i][j]);
            }
        }
        for(int i=1;i<=m;++i)
        {
            for(int j=0;j<=n;++j)
            {
                win[i][j]=-1.0;//初始化
            }
        }
        int ans=1;
        for(int i=1;i<=m;++i)
        {
            if(pro(i,n,1,m)>win[ans][n])//枚举每队的获胜概率 记录最大的队伍
            ans=i;
        }
        printf("%d\n",ans);
    }
    return 0;
}