poj 3071 Football (概率 dp)

链接:poj 3071

题意:有2^n个队,相邻的两两打淘汰赛,n轮后必定会决出冠军,

      求最后哪个队夺冠的概率最大

分析:dp[i][j]表示第i轮的时候,第j去支队伍赢的概率

则dp[i][j]的前提就是i-1轮的时候,j是赢的,且第i轮赢了对方

接下来就是找到第i轮的时候,他的可能队手(难点

所有对手的编号为1-2^n,通过二进制可以发现规律,

若第j队和对手k队在第i轮比赛,那么(j-1)和(k-1)二进制第i位刚好相反,

从i+1位开始所有高位是一样的,可以利用位运算判断是否为对手

状态方程为:

dp[i][j]=sigma(dp[i-1][j]*dp[i-1][k]*p[j][k])


#include<stdio.h>
#include<string.h>
int main()
{
    double p[150][150],dp[150][150];
    int n,i,j,k;
    while(scanf("%d",&n)!=EOF){
        if(n==-1)
            break;
        int m=1<<n;
        for(i=1;i<=m;i++)
            for(j=1;j<=m;j++)
                scanf("%lf",&p[i][j]);
        memset(dp,0,sizeof(dp));
        for(i=1;i<=m;i++)
            dp[0][i]=1.0;
        for(i=1;i<=n;i++)
            for(j=1;j<=m;j++)
                for(k=1;k<=m;k++)
                    if(((j-1)>>(i-1)^1)==((k-1)>>(i-1)))  //判断是否为对手
                        dp[i][j]+=dp[i-1][j]*dp[i-1][k]*p[j][k];
        double max_p=0;
        int ans=0;
        for(i=1;i<=m;i++)
            if(dp[n][i]>max_p){
                max_p=dp[n][i];
                ans=i;
            }
        printf("%d\n",ans);
    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值