Problem NYOJ-712 探寻宝藏 （双向DP）

探寻宝藏 （双向DP）from第六届河南省程序设计大赛

NYOJ-712 ZZULI oj-1484
题目链接
最近南阳有点卡，提交的话后者还行，回到主题 ，中文易懂，但是有想法实现不了，最后才知道 得get （双向DP）这个技能。
首先 想 dp[i][j]这个状态表示啥，依题意 ，求最多的宝藏 ，就相当于之前做过的从 某位置 到 另一位置 的 二维数表里的获得的最大价值（每一个数都有它的价值 —— a[i][j] ）：

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

开始是想对了，然后后面想着判断条件，怎么走回去……其实，并不然，在这条路的基础上 再来一条路：

dp[p][q] = max( dp[p+1][q] , dp[p][q+1] ) + a[p][q]

把 i，j，p，q 联系组合在一起就是：
在（i+j = p+q）的条件下，
dp[i][j][p][q] = max( dp[i-1][j][p-1][q] , dp[i][j-1][p][q-1]
, dp[i-1][j][p][q-1] , dp[i][j-1][p-1][q])
+ a[i][j] + a[p][q];

备注：a[0][0] 入口 一般开始为0不考虑，最后别忘记 加上 a[m][n]

具体的代码如下：

#include<cstdio>
#include<cstring>
using namespace std;
const int M =55;
int a[M][M],dp[M][M][M][M];
int mmax(int x,int y){
    return x>y?x:y;
}
int main(){
    int k,m,n,i,j,q,p,maxx;
    scanf("%d",&k);
    while(k--){
        memset(a,0,sizeof(a));

        /*dp[i][j] = max( dp[i+1][j] , dp[i][j+1] ) + a[i][j]   
                i+j = p+q
        dp[p][q] = max( dp[p+1][q] , dp[p][q+1] ) + a[p][q]     */

        scanf("%d %d",&m,&n);
        for(i=1; i<=m; i++){
            for(j=1; j<=n; j++){
                scanf("%d",&a[i][j]);
            }
        }
        memset(dp,0,sizeof(dp));
        int sum = 0;
        for(i=1;i<=m;i++){
            for(j=1;j<=n;j++){
                for(p=i+1;p<=m;p++){
                    q=i+j-p;
                    if(q<=0)  continue;
                    dp[i][j][p][q] = mmax( mmax(dp[i-1][j][p-1][q],dp[i][j-1][p][q-1])
                                            , mmax(dp[i-1][j][p][q-1],dp[i][j-1][p-1][q]) ) + a[i][j] + a[p][q];

                }
            }   
        }
     sum = mmax(mmax(dp[m-1][n][m-1][n],dp[m-1][n][m][n-1]),  
                     mmax(dp[m][n-1][m-1][n],dp[m][n-1][m][n-1]));  

        printf("%d\n",sum + a[m][n]);
    }
    return 0;
}/*
2
2 3
0 10 10
10 10 80
3 3
0 3 9
2 8 5
5 7 100
*/