算法训练（VIP）——方格取数

题目描述
设有 N * N 的方格图, 我们将其中的某些方格中填入正整数,而其他的方格中则放入数字0。

某人从图的左上角的 A 点 (1, 1) 出发，可以向下行走，也可以向右走，直到到达右下角的 B 点 (N, N)。

在走过的路上，他可以取走方格中的数（取走后的方格中将变为数字0）。

此人从 A 点到 B 点共走两次，试找出2条这样的路径，使得取得的数之和为最大。

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

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

样例输入
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

样例输出
67

数据范围
0 < N ≤ 10

---

题解一
记忆化搜索：

#include <iostream>
#include <cstring>
using namespace std;

const int N = 15;

int n, ans;
int w[N][N];
int f[N][N][N][N];

int dfs(int x1, int y1, int x2, int y2)
{	
	if(x1 > n || y1 > n || x2 > n || y2 > n) return 0;
	if(f[x1][y1][x2][y2] != -1) return f[x1][y1][x2][y2];
	
	int t = w[x1][y1];
	if(x1 != x2 || y1 != y2) t += w[x2][y2];                    // 走的不是同一个点
	
	int &x = f[x1][y1][x2][y2];                                 // 用 x 代替 f[x1][y1][x2][y2]，简化代码 
	x = max(x, dfs(x1 + 1, y1, x2 + 1, y2) + t);                
	x = max(x, dfs(x1 + 1, y1, x2, y2 + 1) + t);                
	x = max(x, dfs(x1, y1 + 1, x2 + 1, y2) + t);               
	x = max(x, dfs(x1, y1 + 1, x2, y2 + 1) + t);                
	
	return x;
}

int main()
{
	cin >> n;
	
	int a, b, c;
	while(cin >> a >> b >> c, a || b || c) w[a][b] = c;
	
	memset(f, -1, sizeof f);
	
	cout << dfs(1, 1, 1, 1) << endl;
	return 0;
}

题解二
动态规划：

题目要求从 A 点到 B 点走两次，使路径之和最大。我们可以将它转换为：两个人同时从 A 点出发，然后分别走到 B 点。

一、f[x1][y1][x2][y2]：

• 集合：所有从（1, 1）开始，分别走到 （x1, y1），（x2, y2）的方案的集合
• 属性：最大值

二、f[x1][y1][x2][y2]的上一步：

1. 上上：第一个人从（x1 - 1, y1）走过来，第二个人从（x2 - 1, y2）走过来。
2. 上左：…
3. 左上：…
4. 左左：…

#include <iostream>
using namespace std;

const int N = 15;

int w[N][N];
int f[N][N][N][N];

int main()
{
	int n;
	cin >> n;
	
	int a, b, c;
	while(cin >> a >> b >> c, a || b || c) w[a][b] = c;
	
	for (int x1 = 1; x1 <= n; x1 ++)
		for (int y1 = 1; y1 <= n; y1 ++)
			for (int x2 = 1; x2 <= n; x2 ++)
				for (int y2 = 1; y2 <= n; y2 ++)
				{
					if(x1 + y1 != x2 + y2) continue;
					
					int t = w[x1][y1];
					if(x1 != x2 || y1 != y2) t += w[x2][y2];
					
					int &x = f[x1][y1][x2][y2];
					x = max(x, f[x1 - 1][y1][x2 - 1][y2] + t);
					x = max(x, f[x1 - 1][y1][x2][y2 - 1] + t);
					x = max(x, f[x1][y1 - 1][x2 - 1][y2] + t);
					x = max(x, f[x1][y1 - 1][x2][y2 - 1] + t);
				}
	
	cout << f[n][n][n][n] << endl;
	return 0;
}

题解三
动态规划（优化）：
f[k][x1][x2]

• 集合：所有从（1， 1）开始，分别走到（x1， k - x1），（x2，k - x2）的方案的集合
• 属性：最大值

k：当前的横纵坐标之和

#include <iostream>
using namespace std;

const int N = 15;

int w[N][N];
int f[N + N][N][N];

int main()
{
	int n;
	cin >> n;
	
	int a, b, c;
	while(cin >> a >> b >> c, a || b || c) w[a][b] = c;
	
	for (int k = 2; k <= n + n; k ++)
		for (int x1 = 1; x1 < k; x1 ++)
			for (int x2 = 1; x2 < k; x2 ++)
			{
				int y1 = k - x1, y2 = k - x2;
				if(y1 < 1 || y1 > n || y2 < 1 || y2 > n) continue;
				
				int t = w[x1][y1];
				if(x1 != x2) t += w[x2][y2];
				
				int &x = f[k][x1][x2];
				x = max(x, f[k - 1][x1 - 1][x2 - 1] + t);
				x = max(x, f[k - 1][x1 - 1][x2] + t);
				x = max(x, f[k - 1][x1][x2 - 1] + t);
				x = max(x, f[k - 1][x1][x2] + t);
			}
	
	cout << f[n + n][n][n] << endl;
	return 0;
}