【SSL】1010方格取数

【SSL】1010方格取数

Time Limit:50000MS

Memory Limit:65536K

Description

设有N*N的方格图(N<=10,我们将其中的某些方格中填入正整数,而其他的方格中则放入数字0。如下图所示（见样例）：

　　某人从图的左上角的A 点出发，可以向下行走，也可以向右走，直到到达右下角的B点。在走过的路上，他可以取走方格中的数（取走后的方格中将变为数字0）。
　　此人从A点到B 点共走两次，试找出2条这样的路径，使得取得的数之和为最大。

Input

输入的第一行为一个整数N（表示N*N的方格图），接下来的每行有三个整数，前两个表示位置，第三个数为该位置上所放的数。一行单独的0表示输入结束。

Output

只需输出一个整数，表示2条路径上取得的最大的和。

Sample Input

8
2  3  13
2  6   6
3  5   7
4  4  14
5  2  21 
5  6   4
6　3  15
7　2  14
0　0  0

Sample Output

67

思路

实际上本问题完全可以用动态规划解决，只是递推起来更为复杂些而已，前面在考虑只走一次的情况，只需考虑一个人到达某个格子（i，j）的情况，得出f[i][j]=max(sum[i-1][j],sum[i][j-1])+a[i][j]，现在考虑两个人同时从A出发，则需考虑两个人到达任意两个格子(i1，j1)与(i2，j2)的情况，显然要到达这两个格子，其前一状态必为(i1-1，j1)，(i2-1，j2)；(i1-1，j1)，(i2，j2-1)；(i1，j1-1)，(i2-1，j2)；(i1，j1-1)，(i2，j2-1) 四种情况之一，类似地，可以推导出：设p=max(f[i1-1][j1][i2-1][j2],f[i1-1][j1][i2][j2-1],f[i1][j1-1][i2-1][j2],f[i1][j1-1][i2][j2-1])，则
i1=i2&&j1=j2:f[i1][j1][i2][j2]=p+a[i][j]
i1!=i2 || j1!=j2:f[i1][j1][i2][j2]=p+f[i1][j1]+a[i2][j2]
1<=i1,j1,i2,j2<=n

代码

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int n,a[55][55],f[55][55][55][55];
void input()
{
	int x,y,w;
	memset(a,0,sizeof(a));
	scanf("%d",&n);
	for(scanf("%d%d%d",&x,&y,&w);x||y||w;scanf("%d%d%d",&x,&y,&w))
		a[x][y]=w;
	return;
}
void DP()
{
	int i1,j1,i2,j2;
	for(i1=1;i1<=n;i1++)
		for(j1=1;j1<=n;j1++)
			for(i2=1;i2<=n;i2++)
				for(j2=1;j2<=n;j2++)
					if (i1==i2&&j1==j2)//状态转移方程
						f[i1][j1][i2][j2]=max(max(f[i1-1][j1][i2-1][j2],f[i1-1][j1][i2][j2-1]),max(f[i1][j1-1][i2-1][j2],f[i1][j1-1][i2][j2-1]))+a[i1][j1];
					else
						f[i1][j1][i2][j2]=max(max(f[i1-1][j1][i2-1][j2],f[i1-1][j1][i2][j2-1]),max(f[i1][j1-1][i2-1][j2],f[i1][j1-1][i2][j2-1]))+a[i1][j1]+a[i2][j2];
	return;
}
int main()
{
	input();
	DP();
	printf("%d",f[n][n][n][n]);
	return 0;
 }