HDU 2686 Matrix多线程DP

版权声明:转载请注明文章出处 https://blog.csdn.net/zhj_fly/article/details/79174949

题意:给定一个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;
}

猜你喜欢

转载自blog.csdn.net/zhj_fly/article/details/79174949