NBOJv2——Problem 1002: 蛤玮的财宝（多线程DP）

Problem 1002: 蛤玮的财宝

Time Limits: 1000 MS Memory Limits: 65536 KB

64-bit interger IO format: %lld Java class name: Main

Description

蛤玮和他的妹子出海游玩,不小心遭遇了海难,他们醒来之后发现自己到了一座金银岛.岛主非常好心的告诉他们在岛的另一边有船可以送他们回家.

这座岛可以看成n*m的矩阵,蛤玮他们在位置(1,1),而船在位置(n,m).蛤玮发现金银岛遍地都是金子,每个格子里有价值a[i,j]的金子,他和妹子打算在回去的路上带一些走.如果他们路过了位置(i,j),就可以假装系鞋带捡走地上的金子.为了不引起怀疑,他们在走的时候只能往接近码头的方向走,即如果蛤玮现在在(i,j),他只能移动到(i+1,j)或者(i,j+1).为了能拿走更多的金子,蛤玮和妹子决定装作互相不认识,这样他们就可以分开走,从而拿到更多的金子.

蛤玮和他妹子想知道他们最多能拿走多少金子.

注意如果蛤玮和他妹子经过了相同的地方,只能得到一次金子,因为地上的捡完就没有了.

Input

T(1<=T<=10),表示数据组数.

每组数据第一行n,m(1<=n,m<=100),接下来n行,每行m个数,第i行第j列的值a[i,j](1<=a[i,j]<=1000)表示位置(i,j)的金子的价值.

Output

每组数据输出一行,蛤玮和他妹子能拿到的金子总价的最大值.

Sample Input

Output for Sample Input

一个人的情况下还是比较好写的，两个人就完全懵比了。看了别人的讲解才知道怎么写，两个人就用dp[i][j][ii][jj]来表示两个人共同状态，然后转移方程就是

dp[i][j][ii][jj]=max(dp[i-1][j][ii-1][jj]，dp[i-1][j][ii][jj-1]，dp[i][j-1][ii-1][jj]，dp[i][j-1][ii][jj-1])+value[i][j]+(i==ii&&j==jj)?0:value[ii][jj]；

但是这道题显然n、m四维来表示会爆内存，然后就有一种优化像二维DP优化到一维一样，把四维优化到三维。

就有

dp[k][ii][jj]=max(dp[k-1] [ii-1][jj]，dp[k-1] [ii][jj-1]，dp[k-1] [ii-1][jj-1]，dp[k-1] [ii][jj])+value[k-i+1][i]+

(ii==jj)?0:value[ii][jj]；

然后就不会RE了。中间循环的时候i与j变量范围注意下，不然下标可能就变成了负的了，RE两次……

代码：

#include<iostream>

#include<algorithm>

#include<cstdlib>

#include<sstream>

#include<cstring>

#include<cstdio>

#include<string>

#include<deque>

#include<stack>

#include<cmath>

#include<queue>

#include<set>

#include<map>

#define INF 0x3f3f3f3f

#define MM(x) memset(x,0,sizeof(x))

using namespace std;

typedef long long LL;

int pos[202][102];

int dp[202][102][102];

int main(void)

{

	int n,m,i,j,tcase,k;

	scanf("%d",&tcase);

	while (tcase--)

	{

		scanf("%d%d",&n,&m);

		for(i=1; i<=n; i++)

		{

			for (j=1; j<=m; j++)

			{

				scanf("%d",&pos[i][j]);

			}

		}

		for (k=1; k<=n+m; k++)

		{

			for (i=1; i<=k; i++)

			{

				for (j=1; j<=k; j++)

				{

					int temp=max(max(dp[k-1][i-1][j-1],dp[k-1][i-1][j]),max(dp[k-1][i][j-1],dp[k-1][i][j]));

					if(i==j)

						dp[k][i][j]=temp+pos[k-i+1][i];

					else

						dp[k][i][j]=temp+pos[k-i+1][i]+pos[k-j+1][j];

				}

			}

		}

		printf("%d\n",dp[n+m][n][m]);

		memset(dp,0,sizeof(dp));

	}

	return 0;

}