51nod 1084：矩阵取数问题 V2

1084 矩阵取数问题 V2

基准时间限制：2 秒 空间限制：131072 KB 分值: 80  难度：5级算法题

 收藏

 关注

一个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来做的话，自己一开始数塔取数，然后将走过的路径置0，再取一次，发现这样做不对。举例：

1 1 3

3 3 3

3 1 1

然后自己反应过来，这个问题实际上就是POJ3422http://blog.csdn.net/u010885899/article/details/50330303的子问题，这里的K=2了，用网络流跑一遍就得到结果。

实际上这里用DP[k][m][n]表示两个人在走取到的数的最大值。已经走了k步，第一个人走到第m行，第二个人走到第n行，通过走的k步，可知其列号。然后两个人走就有四个方向了，下下、下右、右下、右右。取一遍最大值。

对于位置相等的时候，判定一下，只取一个值。

代码：

#pragma warning(disable:4996)  
#include <iostream>  
#include <algorithm>  
#include <cmath>  
#include <vector>  
#include <string>  
#include <cstring>
#include <queue>
#include <map>
using namespace std;
typedef long long ll;

int res;
int m, n;

int A[205][205];
int dp[405][205][205];

void input()
{
	int i, j;
	scanf("%d%d", &m, &n);

	for (i = 1; i <= n; i++)
	{
		for (j = 1; j <= m; j++)
		{
			scanf("%d", &A[i][j]);
		}
	}
}


void solve()
{
	int i, j, k;

	for (i = 2; i <= n + m + 1; i++)
	{
		for (j = 1; j <= n&&i - j >= 0; j++)
		{
			for (k = 1; k <= n&&i - k >= 0; k++)
			{
				if (j == k)
				{
					dp[i][j][k] = max(dp[i][j][k], dp[i - 1][j - 1][k - 1] + A[j][i - j]);//都是往下走才到达的这里
					dp[i][j][k] = max(dp[i][j][k], dp[i - 1][j][k - 1] + A[j][i - j]);
					dp[i][j][k] = max(dp[i][j][k], dp[i - 1][j - 1][k] + A[j][i - j]);
					dp[i][j][k] = max(dp[i][j][k], dp[i - 1][j][k] + A[j][i - j]);
				}
				else
				{
					dp[i][j][k] = max(dp[i][j][k], dp[i - 1][j - 1][k - 1] + A[j][i - j] + A[k][i - k]);// 都是往下走才到达的这里
					dp[i][j][k] = max(dp[i][j][k], dp[i - 1][j][k - 1] + A[j][i - j] + A[k][i - k]);
					dp[i][j][k] = max(dp[i][j][k], dp[i - 1][j - 1][k] + A[j][i - j] + A[k][i - k]);
					dp[i][j][k] = max(dp[i][j][k], dp[i - 1][j][k] + A[j][i - j] + A[k][i - k]);
				}
			}
		}
	}
	cout << dp[n + m + 1][n][n] << endl;
}

int main()
{
	//freopen("i.txt", "r", stdin);
	//freopen("o.txt", "w", stdout);

	input();
	solve();

	return 0;
}