P1006 [NOIP2008 提高组] 传纸条

题目描述

小渊和小轩是好朋友也是同班同学，他们在一起总有谈不完的话题。一次素质拓展活动中，班上同学安排坐成一个 mm 行 nn 列的矩阵，而小渊和小轩被安排在矩阵对角线的两端，因此，他们就无法直接交谈了。幸运的是，他们可以通过传纸条来进行交流。纸条要经由许多同学传到对方手里，小渊坐在矩阵的左上角，坐标 (1,1)(1,1)，小轩坐在矩阵的右下角，坐标 (m,n)(m,n)。从小渊传到小轩的纸条只可以向下或者向右传递，从小轩传给小渊的纸条只可以向上或者向左传递。

在活动进行中，小渊希望给小轩传递一张纸条，同时希望小轩给他回复。班里每个同学都可以帮他们传递，但只会帮他们一次，也就是说如果此人在小渊递给小轩纸条的时候帮忙，那么在小轩递给小渊的时候就不会再帮忙。反之亦然。

还有一件事情需要注意，全班每个同学愿意帮忙的好感度有高有低（注意：小渊和小轩的好心程度没有定义，输入时用 00 表示），可以用一个 [0,100][0,100] 内的自然数来表示，数越大表示越好心。小渊和小轩希望尽可能找好心程度高的同学来帮忙传纸条，即找到来回两条传递路径，使得这两条路径上同学的好心程度之和最大。现在，请你帮助小渊和小轩找到这样的两条路径。

输入格式

第一行有两个用空格隔开的整数 mm 和 nn，表示班里有 mm 行 nn 列。

接下来的 mm 行是一个 m \times nm×n 的矩阵，矩阵中第 ii 行 jj 列的整数表示坐在第 ii 行 jj 列的学生的好心程度。每行的 nn 个整数之间用空格隔开。

输出格式

输出文件共一行一个整数，表示来回两条路上参与传递纸条的学生的好心程度之和的最大值。

输入输出样例

输入 #1复制

3 3
0 3 9
2 8 5
5 7 0

输出 #1复制

34

 

说明/提示

【限制】

对于 30\%30% 的数据，1 \le m,n \le 101≤m,n≤10；
对于 100\%100% 的数据满足：1 \le m,n \le 501≤m,n≤50。

/*
传纸条问题==四维动态规划问题dp
相当于一个纸条传两编
f[i][j][p][q]为第一张纸条传到(i,j),第二张纸条传到(p,q)所累计下来的好心程度和

f[i][j][p][q]=max(f[i-1][j][p][q-1],f[i-1][j][p-1][q],f[i][j-1][p-1][q],f[i][j-1][p][q-1])+a[i][j]+a[p][q]

*/

可以先做P1004 [NOIP2000 提高组] 方格取数https://www.luogu.com.cn/problem/P1004

同为四维动态规划，不过这个题目讲解比较明白

代码：

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include <iostream>
#define maxnum 55 
using namespace std;

/*
传纸条问题==四维动态规划问题dp
相当于一个纸条传两编
f[i][j][p][q]为第一张纸条传到(i,j),第二张纸条传到(p,q)所累计下来的好心程度和

f[i][j][p][q]=max(f[i-1][j][p][q-1],f[i-1][j][p-1][q],f[i][j-1][p-1][q],f[i][j-1][p][q-1])+a[i][j]+a[p][q]

*/
int Max(int a, int b, int c, int d)
{
	int  res = 0;
	res = max(a, b);
	res = max(res, max(c, d));
	return res;
}
int main()
{
	int m,n;
	cin >> m >> n;
	/*
	int** num = new int* [50+ 2];
	for (int i = 0; i < 50 + 2; i++)
	{
		num[i] = new int[n + 2];
		memset(num[i], 0, sizeof(num[i]));
	}
	*/
	int num[52][52] = { 0 };
	for (int i = 1; i <= m; i++)
	{
		for (int j = 1; j <= n; j++)
		{
			cin >> num[i][j];
		}
	}
	//int dp[52][52][52][52] = { 0 };
	int**** dp = new int*** [maxnum];
	for (int i = 0; i < maxnum; i++)
	{
		dp[i] = new int**[maxnum];
		for (int j = 0; j < maxnum; j++)
		{
			dp[i][j] = new int* [maxnum];
			for (int k = 0; k < maxnum; k++)
			{
				dp[i][j][k] = new int[maxnum];
				memset(dp[i][j][k], 0, sizeof(dp[i][j][k]));

			}
		}
	}
	dp[1][1][1][1] = 0;
	for (int i = 1; i <= m; i++)
	{
		for (int j = 1; j <= n; j++)
		{
			for (int p = 1; p <= m; p++)
			{
				for (int q = j + 1; q <= n; q++)
				{
					dp[i][j][p][q] = Max(dp[i - 1][j][p][q - 1], dp[i - 1][j][p - 1][q], dp[i][j - 1][p - 1][q], dp[i][j - 1][p][q - 1]) + num[i][j] + num[p][q];
				}
			}
		}
	}
	cout << dp[m][n - 1][m - 1][n] << endl;
	return 0;

}

 遇到问题：

在vs中直接使用：int dp[50][50][50][50] 会出现如下错误，导致程序运行失败

C6262

• 项目
• 2022/03/24
• 2 个参与者

警告 C6262：函数使用 stack 的 constant_1 字节：超过了/analyze： stacksize constant_2。 考虑将某些数据移到堆

解决办法：

 1、将int dp[50][50][50][50]定义在main函数外部 ，就可以正常运行。

2、使用动态申请数组空间。

如何申请：

//申请四维
int**** dp = new int*** [maxnum];
	for (int i = 0; i < maxnum; i++)
	{
		dp[i] = new int**[maxnum];
		for (int j = 0; j < maxnum; j++)
		{
			dp[i][j] = new int* [maxnum];
			for (int k = 0; k < maxnum; k++)
			{
				dp[i][j][k] = new int[maxnum];
				memset(dp[i][j][k], 0, sizeof(dp[i][j][k]));

			}
		}
	}

申请二维；

/*
	int** map = new int* [n + 2];
	for (int i = 0; i < n + 2; i++)
	{
		map[i] = new int[n + 2];
		memset(map[i],0, sizeof(map[i]));
	}
	*/