POJ 3071 Football(概率DP)

题目链接

不1Y都对不住看过那么多年的球。dp[i][j]表示i队进入第j轮的概率,此题用0-1<<n表示非常方便。

 #include <cstdio>
 #include <cstring>
 #include <iostream>
 #include <cmath>
 #include <algorithm>
 using namespace std;
 double dp[][];
 double p[][];
 int main()
 {
     int n,i,j,mod,c,k,ans;
     double maxz;
     while(scanf("%d",&n)!=EOF)
     {
         if(n < ) break;
         for(i = ;i < (<<n);i ++)
         {
             for(j = ;j < (<<n);j ++)
             {
                 scanf("%lf",&p[i][j]);
             }
         }
         for(i = ;i < (<<n);i ++)
         {
             dp[][i] = ;
         }
         for(i = ;i <= n;i ++)
         {
             for(j = ;j < (<<n);j ++)
             dp[i][j] = ;
         }
         for(i = ;i <= n;i ++)
         {
             for(j = ;j < (<<n);j ++)
             {
                 mod = j%(<<i);
                 if(mod >= (<<(i-)))
                 {
                     c = (j/(<<i))*(<<i);
                     for(k = ;k < (<<(i-));k ++)
                     {
                         dp[i][j] += dp[i-][j]*dp[i-][c+k]*p[j][c+k];
                     }
                 }
                 else
                 {
                     c = (j/(<<i))*(<<i)+(<<(i-));
                     for(k = ;k < (<<(i-));k ++)
                     {
                         dp[i][j] += dp[i-][j]*dp[i-][c+k]*p[j][c+k];
                     }
                 }
             }
         }
         ans = ;maxz = ;
         for(i = ;i < (<<n);i ++)
         {
             if(maxz < dp[n][i])
             {
                 ans = i;
                 maxz = dp[n][i];
             }
         }
         printf("%d\n",ans+);
     }
     return ;
 }