1023: Pong’s Birds（概率）

1023: Pong’s Birds

时间限制: 1 Sec  内存限制: 128 MB
提交: 94  解决: 21
[提交][状态][讨论版]

题目描述

In order to train his birds, Pong is holding a competition for them. (He have birds, does he?) He have 2n birds,
so he want to hold elimination series. Literally, if n = 2 and he has 4 birds identified as 1,2,3,4, he would first hold competition for 1 and 2, 3 and 4, the one fails would be eliminated and then he holds competition for the winners. 
According to his long observation, for each pair of birds he knows the probability to win for each bird, and he want to know the probability of each bird to be the final winner.

输入

For the first row there is an integer T(T ≤ 100), which is the number of test cases.
For each case , first row is an integer n(1 ≤ n ≤ 7), for the next 2n rows , each row has 2n real numbers as the probability to win for the i-th when competed to j-th bird. It is promised that for each i, j p[i][j] + p[j][i] = 1 and p[i][i] = 0;

输出

For each case you should output a line with 2n real numbers, which is the probability of i-th bird be the final winner. Rounded to 3 decimal places.
Make sure to output a space character after each number, to prevent the Presentation Error.

样例输入

1
1
0 0.5
0.5 0

样例输出

0.500 0.500

提示

求队伍夺冠概率，画图找下规律递推，

 #include <bits/stdc++.h>
 using namespace std;

 const int MAXN = ( << ) + ;
 double p[MAXN][MAXN];
 double dp[][MAXN];//dp[i][j]表示第i轮j队伍出线概率

 int main()
 {
     int t;
     int n;
     int n2;
     int i, j, k;
     int grp;
     for (j = ; j < MAXN; ++j) {
         dp[][j] = ;
     }
     scanf("%d", &t);
     while (t--) {
         scanf("%d", &n);
         n2 =  << n;
         for (i = ; i < n2; ++i) {
             for (j = ; j < n2; ++j) {
                 scanf("%lf", &p[i][j]);
             }
         }

         for (i = ; i <= n; ++i) {
             for (j = ; j < n2; ++j) {
                 grp = j / ( << (i - ));//
                 grp ^= ;//奇偶互换
                 dp[i][j] = ;
                 for (k = grp << (i - ); k < ((grp + ) << (i - )); ++k) {
                     dp[i][j] += dp[i - ][j] * dp[i - ][k] * p[j][k];
                 }
             }
         }

         for (i = ; i < n2 - ; ++i) {
             printf("%.3lf ", dp[n][i]);
         }
         printf("%.3lf\n", dp[n][n2 - ]);
     }
     return ;
 }

记忆化搜索

 #include <bits/stdc++.h>
 using namespace std;

 const int MAXN = ( << ) + ;
 double p[MAXN][MAXN];
 double dp[][MAXN];//dp[i][j]表示第i轮j队伍出线概率

 double dfs(int i, int j)
 {
     if (i == ) {
         return ;
     }
     if (dp[i][j] >= ) {
         //printf("%lf\n", dp[i][j]);
         return dp[i][j];
     }
     int k;
     int grp;
     grp = j / ( << (i - ));//
     grp ^= ;//奇偶互换
     dp[i][j] = ;
     for (k = grp << (i - ); k < ((grp + ) << (i - )); ++k) {
         dp[i][j] += dfs(i - , j) * dfs(i - , k) * p[j][k];
     }
     return dp[i][j];
 }

 int main()
 {
     int t;
     int n;
     int n2;
     int i, j, k;
     int grp;
     scanf("%d", &t);
     while (t--) {
         scanf("%d", &n);
         n2 =  << n;
         for (i = ; i < n2; ++i) {
             for (j = ; j < n2; ++j) {
                 scanf("%lf", &p[i][j]);
             }
         }

         //memset(dp, -1, sizeof(dp));//不能使用
         for (i = ; i <= n; ++i) {
             for (j = ; j < n2; ++j) {
                 dp[i][j] = -;
             }
         }

         for (i = ; i < n2 - ; ++i) {
             printf("%.3lf ", dfs(n, i));
         }
         printf("%.3lf\n", dfs(n, n2 - ));
     }
     return ;
 }