nyoj712 - 探 寻 宝 藏 - 双程动态规划（详解）

探 寻 宝 藏

时间限制： 1000 ms  |  内存限制： 65535 KB

难度： 5

描述

传说HMH大沙漠中有一个M*N迷宫，里面藏有许多宝物。某天，Dr.Kong找到了迷宫的地图，他发现迷宫内处处有宝物，最珍贵的宝物就藏在右下角，迷宫的进出口在左上角。当然，迷宫中的通路不是平坦的，到处都是陷阱。Dr.Kong决定让他的机器人卡多去探险。

但机器人卡多从左上角走到右下角时，只会向下走或者向右走。从右下角往回走到左上角时，只会向上走或者向左走，而且卡多不走回头路。（即：一个点最多经过一次）。当然卡多顺手也拿走沿路的每个宝物。

Dr.Kong希望他的机器人卡多尽量多地带出宝物。请你编写程序，帮助Dr.Kong计算一下，卡多最多能带出多少宝物。

输入

第一行： K 表示有多少组测试数据。 
接下来对每组测试数据：
第1行: M N
第2~M+1行： Ai1 Ai2 ……AiN (i=1,…..,m)


【约束条件】
2≤k≤5 1≤M, N≤50 0≤Aij≤100 (i=1,….,M; j=1,…,N)
所有数据都是整数。 数据之间有一个空格。

输出

对于每组测试数据，输出一行：机器人卡多携带出最多价值的宝物数

样例输入

22 30 10 1010 10 803 30 3 92 8 55 7 100

样例输出

120134

来源

第六届河南省程序设计大赛

上传者

ACM_赵铭浩

分析：

可以看成过两个人同时从起点向终点前进，但要保证路线不能相交，除了在终点，两人速度相同，如果某一时刻在相同的行或者列，则他们一定会在相同的列或行，因此保证他们行不相同，既可以保证列不相交

dp[k][i][j]：走到第k步（在起点时是第二步），A走到i行，B走到j行时所得的最大宝藏数，则A在的列为k-i,B为k-j

比如从(1,1)走到(1,4)走了5步，在第1行，所以在第5-1列

则，一共要走m+n-1-1+2步，最后求的是dp[m+n][n][n]

k>=i&&k>=j&&i!=j

状态转移方程：

当l<k的时候  表示走了l步时，A到i行，B到j行的最大值

也就是四个值的最大值加上a[i][l-i]+a[j][l-j]

dp[l-1][i-1][j]      上一步A在i-1行，B在j行

dp[l-1][i][j-1])     上一步A在i行，B在j-1行

dp[l-1][i-1][j-1]   上一步A在i-1行，B在j-1行

dp[l-1][i][j]         上一步A在i行，B在j行

dp[l][i][j]=max(max(dp[l-1][i-1][j],dp[l-1][i][j-1]),max(dp[l-1][i-1][j-1],dp[l-1][i][j]))+a[i][l-i]+a[j][l-j];

求dp[m+n][n][n]时，只有两种情况，AB分别从上或右方到达终点：

dp[k][row][row]=max(dp[k-1][row][row-1],dp[k-1][row-1][row])+a[row][col];

代码如下：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<string>
#include<queue>

using namespace std;
int dp[110][55][55];
int a[55][55],n,m;
void find_ans(int k,int row,int col){
    for(int l=2;l<k;l++){
        for(int i=1;i<=row;i++){
            for(int j=1;j<=row;j++){
                if(i!=j&&l>=i&&l>=j){
                    dp[l][i][j]=max(max(dp[l-1][i-1][j],dp[l-1][i][j-1]),max(dp[l-1][i-1][j-1],dp[l-1][i][j]))+a[i][l-i]+a[j][l-j];
                }
            }
        }
    }
    dp[k][row][row]=max(dp[k-1][row][row-1],dp[k-1][row-1][row])+a[row][col];
    printf("%d\n",dp[k][row][row]);
}
int main(){
    int k;
    scanf("%d",&k);
    while(k--){
        scanf("%d%d",&n,&m);
        for(int i=1;i<=n;i++){
            for(int j=1;j<=m;j++){
                scanf("%d",&a[i][j]);
            }
        }
        if(n==0||m==0){printf("0\n");continue;}
        memset(dp,0,sizeof(dp));
        find_ans(n+m,n,m);
    }
    return 0;
}