P1004 [NOIP2000 提高组] 方格取数（三维线性dp 数字三角形模型）

设有 N×N 的方格图，我们在其中的某些方格中填入正整数，而其它的方格中则放入数字0。如下图所示：

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

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

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

输入格式
第一行为一个整数N，表示 N×N 的方格图。

接下来的每行有三个整数，第一个为行号数，第二个为列号数，第三个为在该行、该列上所放的数。

行和列编号从 1 开始。

一行“0 0 0”表示结束。

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

数据范围
N≤10

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

题意：

题目给定一个 n × n 的矩阵，矩阵中的每个格子上有一个价值为 w 的物品，给定起点 (1, 1)，终点 (n, n)，每次规定只能向右或向下走。

我们可以从起点出发两次，但同一个格子在两次路径中都经过的话，物品价值只会被累加 1 次。

问：两次从 A 到 B，途径格子的物品的总价值最大是多少？

本题依旧是由数字三角形 衍生而来。

思路：

两个路线同时进行

如果我们的思路是分两次走：

分开两次走（贪心）：第一次走到(n, n)求出最大值并记录路径令路径上点收益为0后再走一次。
第一次走为局部最优并且也对第二次走造成了影响，第二次走是在第一次影响下所能走的局部最优，不具备“无后效性”，因此分开两次走并不是全局最优解

三维dp（由四维dp优化而来）

一、状态表示：dp[i1][j1][i2][j2]（四维）

①集合：所有从(1, 1)，(1, 1)分别走到(i1, j1)，(i2, j2)的路径

②属性：两条路径和的最大值

二、优化：四维化三维

令k = i1 + j1 = i2 + j2: 两人同时走，每个人走的步数是一样的。(两人同时走才可以获得全局最优解)推出三维的状态转移方程：

dp[i1,k−i1,i2,k−i2]推出dp[k,i1,i2]:两人同时走k步,从(1, 1)，(1, 1)走到(i1, j1)，(i2, j2)能获得的最大花生数目。

题目要求格子只能够取一次，当两人在同一个格子时有i1 = i2 && j1 = j2，所以当i1 = i2时两人走到同一格。

dp[i1−1,j1,i2−1,j2]推出dp[k−1,i1−1,i2−1]：代表两人都走了k−1步，1号最后一步要从(i1−1, j1)到(i1, j1)，2号要从(i2−1, j2)到(i2, j2)。

所以需要判断(i1, j1)，(i2, j2)是否是同一个格子,如果是,需要加上一个权重,否则两个都需要。
（这里就保证了同一点只能走一次：“取走后的方格中将变为数字0”）

（注意：两次取同一个格子的时候满足条件:k = i1 + j1 = i2 + j2，i1 = i2因此 j1 = k − i1 = j2 = k − i2，w[i1][j2] = w[i2][j2]）

下方 i 即 i1，j 即 i2

三、状态计算：集合划分
（由于有两个人走两条路径，以最后一步划分，集合内部有四种情况）

情况：

一：①路最后一步是往下走到(i1,j1)，②路最后一步是往下走到(i2,j2)
二：①路最后一步是往下走到(i1,j1)，②路最后一步是往右走到(i2,j2)
三：①路最后一步是往右走到(i1,j1)，②路最后一步是往下走到(i2,j2)
四：①路最后一步是往右走到(i1,j1)，②路最后一步是往右走到(i2,j2)

因此由上我们可以推出状态转移方程：
dp[k][i1][i2] = max(dp[k][i1][i2], dp[k-1][i1-1][i2-1] + t, dp[k-1][i1][i2-1] + t, dp[k-1][i1-1][i2] + t, dp[k-1][i1][i2] + t)
（t表示两条路最后到达的点的权值，如果相同取其一，如果不相同两个都要相加，具体见代码）

下方 i 即 i1，j 即 i2

#include<bits/stdc++.h>

using namespace std;

int n;
int r, c, a;
const int N = 15;
int w[N][N];
int dp[N<<1][N][N];//第一维大小为2*N

int main()
{
	ios_base::sync_with_stdio(false), cin.tie(0);
	cin >> n;
	while (cin >> r >> c >> a)
	{
		if (r == 0 && c == 0 && a == 0) break;
		w[r][c] = a;
	}
	for (int k = 2; k <= n + n; ++k)
	{
		for (int i1 = 1; i1 <= n; ++i1)
		{
			for (int i2 = 1; i2 <= n; ++i2)
			{
			    int j1 = k - i1, j2 = k - i2;
				if (j1 >= 1 && j1 <= n && j2 >= 1 && j2 <= n)//不要超过边界
				{
					int t = w[i1][j1];
					if (i1 != i2) t += w[i2][j2];
					int& x = dp[k][i1][i2];
					x = max(x, dp[k - 1][i1 - 1][i2 - 1] + t);
					x = max(x, dp[k - 1][i1][i2 - 1] + t);
					x = max(x, dp[k - 1][i1 - 1][i2] + t);
					x = max(x, dp[k - 1][i1][i2] + t);
				}
			}
		}
	}
	cout << dp[2 * n][n][n] << endl;
	return 0;
}