题意
给你一个 n*n(1<=n<=300) 的正方形矩阵,现在问你从 (1, 1) 走到 (n, n),再从 (n, n) 走到 (1, 1),路径上所有值的和最大是多少?矩阵中每个元素最多取一次
走法限制:
1、从 (1, 1) 到 (n, n) :只能向下或者向右走
2、从 (n, n) 到 (1, 1) :只能向上或者向左走
做法分析
经典动态规划的变形,其实变的简单了。。。人家原来的是要求两条路径不能重合,现在变成了可以重合,当然,只是考虑的东西比较少了而已
先来看一组样例,n=5:
我们可以把问题转化为两个人同时在 (1, 1) ,向着 (n, n) 走,走的过程中如果走到了相同的地方,那个地方的值就只取一次,这样,我们不难得出上面的图来:每一步所走到的位置就那么几个,阶段的划分就是现在走到第几步了,状态就可以表示成: f[i][a][b] :走到了第 i 步,第一个人处的位置是 a, 第二个人处的位置是 b,其中, a、b是按照上面第一幅图中箭头的方向编号,也就是如第二幅图中所展示的那样,a、b 取的就是1,2,,3,4.。。。
对比上面的图,状态之间的转移方程就不难写出来了:
1、 当前步数小于等于 n 的时候:
每个 a 是由上一步的 a 或者 a-1 得到
每个 a 是由上一步的 a 或者 a-1 得到
每个 b 是由上一步的 b 或者 b-1 得到
当然,得判断上一步中 a 和 a-1 、b 和 b-1 是否合法(超出矩阵范围)
组合起来就是 四 个状态
2、当前步数大于 n 的时候:
每个 a 是由上一步的 a 或者 a+1 得到
每个 b 是由上一步的 b 或者 b+1 得到
组合起来就是 四 个状态
分别求解就行了,由于 步数 是 2*n-1 步,每一步最多有 n 个,而 n 的最大值是 300,于是便需要 600*300*300*4 的内存空间,显然是不够的,所以用滚动数组,将内存变成 300*300*4
AC通道
参考代码
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
const int N=305;
const int INT_INF=0x3fffffff;
int val[N][N], n, f[2][2*N][2*N];
int Max(int i)
{
if(i<=n) return i;
else return 2*n-i;
}
int get_val(int i, int a)
{
if(i<=n) return val[i-a+1][a];
else return val[n-a+1][i-n+a];
}
int main()
{
while(~scanf("%d", &n))
{
for(int i=1; i<=n; i++)
for(int j=1; j<=n; j++)
scanf("%d", &val[i][j]);
memset(f, 0, sizeof(f));
f[1][1][1]=val[1][1];
for(int i=2; i<=n; i++)
{
int len=Max(i);
for(int a=1; a<=len; a++)
for(int b=1; b<=len; b++)
{
int cnt=-INT_INF;
if(a!=1 && b!=1) cnt=max(cnt, f[(i+1)%2][a-1][b-1]);
if(a!=1 && b!=len) cnt=max(cnt, f[(i+1)%2][a-1][b]);
if(a!=len && b!=1) cnt=max(cnt, f[(i+1)%2][a][b-1]);
if(a!=len && b!=len) cnt=max(cnt, f[(i+1)%2][a][b]);
if(a!=b) f[i%2][a][b]=cnt+get_val(i, a)+get_val(i, b);
else f[i%2][a][b]=cnt+get_val(i, a);
}
}
for(int i=n+1; i<=2*n-1; i++)
{
int len=Max(i);
for(int a=1; a<=len; a++)
for(int b=1; b<=len; b++)
{
int cnt=-INT_INF;
cnt=max(cnt, f[(i+1)%2][a+1][b]);
cnt=max(cnt, f[(i+1)%2][a+1][b+1]);
cnt=max(cnt, f[(i+1)%2][a][b]);
cnt=max(cnt, f[(i+1)%2][a][b+1]);
if(a!=b) f[i%2][a][b]=cnt+get_val(i, a)+get_val(i, b);
else f[i%2][a][b]=cnt+get_val(i, a);
}
}
printf("%d\n", f[1][1][1]);
}
return 0;
}