AcWing 1015. 摘花生（线性dp 数字三角形模型）

Hello Kitty想摘点花生送给她喜欢的米老鼠。

她来到一片有网格状道路的矩形花生地(如下图)，从西北角进去，东南角出来。

地里每个道路的交叉点上都有种着一株花生苗，上面有若干颗花生，经过一株花生苗就能摘走该它上面所有的花生。

Hello Kitty只能向东或向南走，不能向西或向北走。

问Hello Kitty最多能够摘到多少颗花生。

输入格式
第一行是一个整数T，代表一共有多少组数据。

接下来是T组数据。

每组数据的第一行是两个整数，分别代表花生苗的行数R和列数 C。

每组数据的接下来R行数据，从北向南依次描述每行花生苗的情况。每行数据有C个整数，按从西向东的顺序描述了该行每株花生苗上的花生数目M。

输出格式
对每组输入数据，输出一行，内容为Hello Kitty能摘到得最多的花生颗数。

数据范围
1≤T≤100,
1≤R,C≤100,
0≤M≤1000

输入样例：
2
2 2
1 1
3 4
2 3
2 3 4
1 6 5

输出样例：
8
16

本题的思考方式和数字三角形很相似，可以看做是 数字三角形 的衍生题目。

对于这一类路线问题，我们状态表示一般可以用点的坐标来表示，我们使用二维数组。

dp[i][j]的状态表示：

①集合：所有从(1, 1)走到(i, j)的路径

②属性（一般有 min, max, count）：max

综合一下，dp[i][j]表示所有从(1, 1)走到(i, j)所有路线收集到花生数量的max值

dp[i][j]的状态计算（对应集合的划分）：

对于集合划分，有一个很重要的划分依据，我们常常以“最后一步”进行划分

我们在草稿纸上可以用一个椭圆来帮助理解集合划分：

椭圆左边部分所有方案都是先从左上角沿着某些路径走，先走到[i-1][j]，然后往下方走，走到[i][j]，右边部分所有集合也都是像这样的一个形式。

对于左边部分最大值max，可以这样分析：

显然每一种这样形式的路线都可以分为两部分：

第一部分是当前[i][j]这个点（所有路线都包含这个点）。

第二部分是从左上角走到[i-1][j]这个点。

因此我们求max时，只需使变化的部分取max即可，即第二个部分（第二部分最大值含义：从左上角走到[i-1][j]所有路线的最大值）。

“第二部分最大值”根据定义恰好是dp[i-1][j]状态表示（从底部走到[i-1][j]所有路线max）的结果。So，椭圆左边部分max即为dp[i-1][j]+w[i][j]

同理可求椭圆右边部分max：dp[i][j-1]+w[i][j]

上面的思考过程和 数字三角形 几乎一样 。

当椭圆左右两边max都有了，我们可以得出最终结果（也就是状态转移方程）：max(dp[i-1][j]+w[i][j],dp[i][j-1]+w[i][j])

这样就知道了每个状态是由那些状态计算转移而来的，如dp[i][j]由dp[i-1][j]、dp[i][j-1]二者转移而来，同样dp[i-1][j]也可以由某些状态转移而来，注意，可能有某个状态可以同时作为两个状态的转移源头。这些过程可以画成一个拓扑图的形式，图可能会有交叉。

我们的目标是求出dp[n][m]，即从左上方走到(n, m)的所有路线可以摘得花生的最大值。

dp的状态转移本质和顺序问题：所有的状态可以看做一个拓扑图中的各个节点，我们需要按照一个拓扑序来计算，对于线性dp我们可以按照某一个合理的顺序来保证按照拓扑序来循环，对于这道题我们只需按照i从小到大循环，且j从小到大循环就可以保证按照拓扑序来循环，这样一来，算每个状态的时候，它所依赖的状态就可以保证被提前计算过。

时间复杂度 ： O(n^2)

#include<bits/stdc++.h>

using namespace std;

int t;
int r, c;
const int N = 100 + 10;
int a[N][N];
int dp[N][N];

int main()
{
	cin >> t;
	while (t--)
	{
		cin >> r >> c;
		for (int i = 1; i <= r; ++i)
		{
			for (int j = 1; j <= c; ++j)
			{
				cin >> a[i][j];
			}
		}
		dp[1][1] = a[1][1];
		for (int i = 1; i <= r; ++i)
		{
			for (int j = 1; j <= c; ++j)
			{
				dp[i][j] = max(dp[i - 1][j], dp[i][j - 1]) + a[i][j];
			}
		}
		cout << dp[r][c] << endl;
	}
	return 0;
}