HDU 2686 Matrix多线程DP

题意：给定一个n*n的矩阵，每个格子中有一个数字，要求从最左上角走到最右下角（只能往下和往右走），再返回最左上角（只能往上或往左走），求经过的数字之和的最大值。要求来回不能重复某个格子。

解析：由于矩阵大小为2<n<30，所以可以使用多线程dp来做。将来回两次变为一次同时两条路径来做。用dp[k][x1][y1][x2][y2]来表示到目前为止走了k步，且第一条路径到（x1, y1），第二条路径到了（x2, y2），此时的最大值。那么转移方程可以表示为：

dp[k][x1][y1][x2][y2] = max(dp[k-1][x1-1][y1][x2-1][y2],dp[k-1][x1-1][y1][x2][y2-1],dp[k-1][x1][y1-1][x2-1][y2],dp[k-1][x1][y1-1][x2][y2-1]) + a[x1][y1] + a[x2][y2];

（这道题n的范围比较小，可以用多线程dp来做，如果n很大的话，就不行了。看大神做法是可以用费用流做的。）

代码：

#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;

const int N = 31;
int a[N][N], dp[N<<1][N][N];
int main(){
	int n;
	while(scanf("%d", &n)!=EOF){
		for(int i = 1; i <= n; i++)
			for(int j = 1; j <= n; j++)
				scanf("%d", &a[i][j]);
		memset(dp, 0, sizeof(dp));
		dp[0][1][1] = a[1][1];
		for(int k = 1; k <= 2*n-2; k++){
			for(int i = 1; i <= n && i <= k + 1; i++)
				for(int j = 1; j <= n && j <= k + 1; j++){
					if(i == j && k != 2*n-2) continue;
					dp[k][i][j] = max(dp[k][i][j], dp[k-1][i][j]);
					dp[k][i][j] = max(dp[k][i][j], dp[k-1][i][j-1]);
					dp[k][i][j] = max(dp[k][i][j], dp[k-1][i-1][j]);
					dp[k][i][j] = max(dp[k][i][j], dp[k-1][i-1][j-1]);
					dp[k][i][j] += a[i][k+2-i] + a[j][k+2-j];
				
				}
		}
		printf("%d\n",dp[2*n-2][n][n]-a[n][n]);
	}
	return 0;
}