ZOJ 3735 Josephina and RPG

实在无法吐槽了。。。先上个图。。。。。最开始看错了N的范围。。以为就是200左右。。结果就是如此坑爹的提交情况。。。哎。。其实还是一个比较简单的dp，，不过还是纠结了好久。。。没办法毕竟太弱了。。。当时看到这个题的时候就觉得是个dp但是觉得搜索可能也能做。。但是试了几次发现TLE。。，果断不想再写了。就开始推转移方程。。。dp【i】【j】  i表示打完第ai【i】组，j表示自己当前为j组，要特判j=ai【i-1】。。,因为有多种选择情况

dp【i】【j】=dp【i-1】【j】*k【j】【ai【i】】；  //j不等于ai【i-1】；

dp【i】【j】=MAX（dp【i-1】【b】*k【b】【ai【i】】）；//j等于ai【i-1】；

#include<iostream>  
#include<cstdio>  
#include<algorithm>  
#include<cstring>  
#include<string>
using namespace std;  
double dp[10005][200];
double k[200][200];
int ai[10005];
int n,m;
int main()
{
    while(~scanf("%d",&n))
    {
        n=n*(n-1)*(n-2)/6;
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<n;j++)
            {
                scanf("%lf",&k[i][j]);
            }
        }
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<n;j++)
            dp[i][j]=0.0;
        }
        scanf("%d",&m);
        for(int i=1;i<=m;i++)
        scanf("%d",&ai[i]);
        for(int i=0;i<n;i++)
        dp[1][i]=k[i][ai[1]];
        double ans=0.0;
        int flag=0;
        for(int i=2;i<=m;i++)
        {
            for(int j=0;j<n;j++)
            {
                if(j!=ai[i-1])
                dp[i][j]=dp[i-1][j]*k[j][ai[i]];   //j不是 ai【i-1】说明dp【i-1】时的状态只能是j 
                else
                {
                    double s=0.0;
                    for(int b=0;b<n;b++)         //j为ai【i-1】，可以是其它组打败j之后才为j的，状态可以有多组 
                    {
                        s=max(s,dp[i-1][b]*k[j][ai[i]]);
                    }
                    dp[i][j]=s;
                }
            }
        }
        for(int k=0;k<n;k++)
            ans=max(ans,dp[m][k]);
        printf("%lf\n",ans);
    }
    return 0;
}