3.20记录 坐标DP

（原题 洛谷P1006）
这题是一个矩阵中的DP，可以用区间DP解决。
区间DP其实就是多个点汇入一个点的操作，整体思路如下：
dp0 = max（dp1，dp2，dp3，…）+ 决策
这题最简单粗暴的办法就是，既然我们一个点(i,j)不能过两次，那干脆开一个四维的dp[i][j][k][l]，只要i=k && j=l就直接认为这个四维的点不可通过就行了。（当然了，多次经过(m,n)是不可避免的，所以这地方要特判）
代码如下：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int dp[51][51][51][51],a[51][51];
int main(){
	int i,j,k,l,n,m;
	scanf("%d %d",&n,&m);
	for(i = 1;i <= n;i++){
		for(j = 1;j <= m;j++) scanf("%d",&a[i][j]);
	}
	for(i = 1;i <= n;i++){
		for(j = 1;j <= m;j++){
			for(k = 1;k <= n;k++){
				for(l = 1;l <= m;l++){
					dp[i][j][k][l] = max(dp[i - 1][j][k - 1][l],max(dp[i - 1][j][k][l - 1],max(dp[i][j - 1][k - 1][l],dp[i][j - 1][k][l - 1]))) + a[i][j] + a[k][l];
					if(i == k && j == l && (i != n || j != m)) dp[i][j][k][l] = -2147483647;
				}
			}
		}
	}
	printf("%d",dp[n][m][n][m]);
	return 0;
}

（这里我按自己习惯写的n行m列，对结果没影响）

还有一种优化思路：矩形有一个性质，就是横纵坐标的增加等同于步数的增加，而且所有在同一矩形中的点都是如此。因此对于x1,x2,x3,…，他们的纵坐标都可以用k-x1,k-x2,k-x3,…来表示。因此只要先统一步数，横坐标就能代替对应的纵坐标（换句话说，步数的循环优先级最大），这可以省去一维。
代码如下：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int dp[101][51][51],a[51][51];
int main(){
	int i,j,k,l,n,m;
	scanf("%d %d",&n,&m);
	for(i = 1;i <= n;i++){
		for(j = 1;j <= m;j++) scanf("%d",&a[i][j]);
	}
	for(i = 1;i < m + n - 1;i++){
		for(j = 1;j <= n;j++){
			for(k = 1;k <= n;k++){
				if(i + 2 - j < 1 || i + 2 - k < 1 || i + 2 - j > m || i + 2 - k > m) continue;
				dp[i][j][k] = max(dp[i - 1][j][k],max(dp[i - 1][j][k - 1],max(dp[i - 1][j - 1][k],dp[i - 1][j - 1][k - 1]))) + a[j][i + 2 - j] + a[k][i + 2 - k];
				if(j == k && (j != n || i + 2 - j != m)) dp[i][j][k] = -2147483647;
			}
		}
	}
	printf("%d",dp[m + n - 2][n][n]);
	return 0;
}

（注意：由于起点我们标注的是(1,1)，故对于任意点(x,y)，其与步数z的关系应该为x+y=z-2）