【更难矩阵取数问题】51nod 1084 矩阵取数问题 V2

Problem Description

一个M*N矩阵中有不同的正整数，经过这个格子，就能获得相应价值的奖励，先从左上走到右下，再从右下走到左上。第1遍时只能向下和向右走，第2遍时只能向上和向左走。两次如果经过同一个格子，则该格子的奖励只计算一次，求能够获得的最大价值。

例如：3 * 3的方格。

1 3 3
2 1 3
2 2 1

能够获得的最大价值为：17。1 -> 3 -> 3 -> 3 -> 1 -> 2 -> 2 -> 2 -> 1。其中起点和终点的奖励只计算1次。

Input

第1行：2个数M N，中间用空格分隔，为矩阵的大小。(2 <= M, N <= 200)
第2 - N + 1行：每行M个数，中间用空格隔开，对应格子中奖励的价值。(1 <= A[i,j] <= 10000)

Output

输出能够获得的最大价值。

Input示例

3 3
1 3 3
2 1 3
2 2 1

Output示例

17

思路：

没招了么？其实我们可以“两个人一起”dp（让两个人同时走）。

用dp[x1][y1][x2][y2]表示第一个人在(x1,y1) 并且第二个人在(x2,y2)时的最大值。

我们有初值dp[1][1][1][1] = a[1][1]， 求的是dp[m][n][m][n]。

问题来了： 每个人走一步，状态转移是什么？

dp[x1][y1][x2][y2] = max{dp[x1’][y1’][x2’][y2’]} + a[x1][y1] + a[x2][y2]
其中(x1’,y1’)是(x1,y1)的邻居，(x2’,y2’)是(x2,y2)的邻居。

事实上，因为我们有这个等式提示我们其实只要用3维就可以表示这个矩阵，因为 y2 = x1 + y1 – x2所以那一维可以用走多少步表示出来。

dp[step + 1][x1][x2] = max{dp[step][x1’][x2’]} + a[x1][y1] + a[x2][y2]
然而这个dp并没有体现出走到相同格子，数字仅计算一次的要求。那么我们加上这个条件：如果x1 = x2，dp[step + 1][x1][x2] = max{dp[step][x1’][x2’]} + a[x1][y1]。

#include<bits/stdc++.h>
using namespace std;
int a[205][205];
int dp[405][205][205];//下标分别代表步数，第一个人x1坐标，第二人的x2坐标，y1 坐标可以由 步数 - x1 + 1弄出来，最开始就算一步。
int main()
{
    int m, n, i, j, x1, x2;
    while(~scanf("%d %d", &m, &n))
    {
        for(i = 1; i <= n; i++)
        {
            for(j = 1; j <= m; j++)
            {
                scanf("%d", &a[i][j]);
            }
        }
        memset(dp, 0, sizeof(dp));//初始化
        for(i = 1; i < m + n; i++)//步数
        {
            for(x1 = 1; x1 <= n && i - x1 >= 0; x1++)//i 如果小于x1就不成立了，这时候y1就是0或者负数，
            {
                for(x2 = 1; x2 <= n && i - x2 >= 0; x2++)
                {
                    int t;
                    if(x1 == x2) t = 0;//如果相等只需要加一次
                    else t = a[x2][i - x2 + 1];//否则两个点的位置都加
                    //两个横坐标，有四种情况
                    dp[i][x1][x2] = max(dp[i][x1][x2], dp[i - 1][x1 - 1][x2 - 1] + a[x1][i - x1 + 1] + t);
                    dp[i][x1][x2] = max(dp[i][x1][x2], dp[i - 1][x1][x2 - 1] + a[x1][i - x1 + 1] + t);
                    dp[i][x1][x2] = max(dp[i][x1][x2], dp[i - 1][x1 - 1][x2] + a[x1][i - x1 + 1] + t);
                    dp[i][x1][x2] = max(dp[i][x1][x2], dp[i - 1][x1][x2] + a[x1][i - x1 + 1] + t);
                }
            }
        }
        printf("%d\n", dp[m + n - 1][n][n]);
    }
}