HDU 4800/zoj 3735 Josephina and RPG 2013 长沙现场赛J题

第一年参加现场赛，比赛的时候就A了这一道，基本全场都A的签到题竟然A不出来，结果题目重现的时候1A，好受打击 ORZ.....

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4800

题目大意：给定C(3,N)支队伍之间对战的获胜概率，再给定一个序列存放队伍编号，每次获胜之后可以选择和当前战胜的对手换队伍。问按给定序列依次挑战全部胜利的最大概率。

解题思路：状压DP

dp[i][j]表示使用队伍i从编号j开始挑战全胜的概率，ai[i]表示i位置的队伍，rate[i][j]表示队伍i战胜队伍j的概率。

状态转移方程：

dp[i][j]=r[i][ai[j]]*max(dp[i][j+1],dp[ai[j]][j+1])

代码如下：

 #include<cmath>
 #include<cstdio>
 #include<cstring>
 #include<iostream>
 #include<algorithm>
 using namespace std;
 #define N 10005

 int m,n;
 int ai[N];
 double dp[][N],rate[][];
 bool vis[][N];

 double search(int team,int lo)
 {
     if(lo>m-||team<||team>n)
         return ;
     if(vis[team][lo]==true)
         return dp[team][lo];
     else
     {
         double t=search(team,lo+);
         if(team!=ai[lo])
         {
             double t2=search(ai[lo],lo+);
             if(t2>t)
                 t=t2;
         }
         vis[team][lo]=true;
         dp[team][lo]=rate[team][ai[lo]]*t;
         return dp[team][lo];
     }
 }

 int main()
 {
     int i,j;
     while(scanf("%d",&n)!=EOF)
     {
         n=n*(n-)*(n-);
         n/=;
         for(i=;i<n;i++){
             for(j=;j<n;j++)
             {
                 scanf("%lf",&rate[i][j]);
             }
         }
         scanf("%d",&m);
         for(i=;i<m;i++)
             scanf("%d",&ai[i]);
         memset(vis,false,sizeof(vis));
         for(i=;i<n;i++)
         {
             dp[i][m-]=rate[i][ai[m-]];
             vis[i][m-]=true;
         }
         double _max=;
         for(i=;i<n;i++)
         {
             _max=max(_max,search(i,));
         }
         printf("%.6lf\n",_max);
     }
     return ;
 }