P1719 最大加权矩形

题目

题目描述

为了更好的备战 NOIP2013，电脑组的几个女孩子 LYQ,ZSC,ZHQ 认为，我们不光需要机房，我们还需要运动，于是就决定找校长申请一块电脑组的课余运动场地，听说她们都是电脑组的高手，校长没有马上答应他们，而是先给她们出了一道数学题，并且告诉她们：你们能获得的运动场地的面积就是你们能找到的这个最大的数字。

校长先给他们一个n×n 矩阵。要求矩阵中最大加权矩形，即矩阵的每一个元素都有一权值，权值定义在整数集上。从中找一矩形，矩形大小无限制，是其中包含的所有元素的和最大 。矩阵的每个元素属于 [−127,127] ,例如

在左下角 

和为 15。

几个女孩子有点犯难了，于是就找到了电脑组精打细算的 HZH，TZY 小朋友帮忙计算，但是遗憾的是他们的答案都不一样，涉及土地的事情我们可不能含糊，你能帮忙计算出校长所给的矩形中加权和最大的矩形吗？

输入格式

第一行：n，接下来是 n 行 n 列的矩阵。

输出格式

最大矩形（子矩阵）的和。

输入输出样例

 

代码

#define _CRT_SECURE_NO_WARNINGS 1
#include<stdio.h>
int a[125][125];
int max(int x, int y)
{
	return x > y ? x: y;
}
int main()
{
	int n = 0;
	scanf("%d", &n);
	int ans = -1e9;
//此题目类似于最大子段和，只不过是二维的形式，先用前缀和将二维数组进行压缩，
//使得每一行的每个元素是其列的所有前面的总和，这样在最后获取某列的n个时就可以用相减的形式获得，
//可以减少for循环的嵌套，减少时间复杂度，
//s[j] = a[i][j] - a[i - k][j];此代码就是获得往上包含k行的方程。

//从开始设置最下面的一行为i行，可以往上包含几行，j代表列数，
//只需要考虑所选行数从最第一列到最后一列所包含的，从左往右扫，就相当于从左上角一直辐射到右角，
//将整个二维数组全部过一遍，最后ans就是最大值。
	for (int i = 1; i <= n; i++)
	{
		for (int j = 1; j <= n; j++)
		{
			scanf("%d",&a[i][j]);
			a[i][j] += a[i - 1][j];//前缀和压缩数组
		}
	}
	for (int i = 1; i <= n; i++)//以第i行为最下面的一行
	{
		for (int k = 1; k <= i; k++)//往上包含几行
		{
			int dp[1000];
			int s[1000];
			for (int j = 1; j <= n; j++)//压缩的第几列
			{
				s[j] = a[i][j] - a[i - k][j];
				dp[j] = max(dp[j - 1] + s[j], s[j]);
				ans = max(ans, dp[j]);
			}

		}
	}
	printf("%d", ans);
	return 0;
}

思路

此题目类似于最大子段和，只不过是二维的形式，先用前缀和将二维数组进行压缩，使得每一行的每个元素是其列的所有前面的总和，这样在最后获取某列的n个时就可以用相减的形式获得，可以减少for循环的嵌套，减少时间复杂度，s[j] = a[i][j] - a[i - k][j];此代码就是获得往上包含k行的方程。

从开始设置最下面的一行为i行，可以往上包含几行，j代表列数，只需要考虑所选行数从最第一列到最后一列所包含的，从左往右扫，就相当于从左上角一直辐射到右下角，将整个二维数组全部过一遍，最后ans就是最大值。

写法二

#define _CRT_SECURE_NO_WARNINGS 1
#include<stdio.h>
int a[125][125];
int max(int x, int y)
{
	return x > y ? x: y;
}
int main()
{
	int n = 0;
	scanf("%d", &n);
	int ans = -1e9;
	for (int i = 1; i <= n; i++)
	{
		for (int j = 1; j <= n; j++)
		{
			scanf("%d",&a[i][j]);
			a[i][j] += a[i - 1][j] + a[i][j - 1] - a[i - 1][j - 1];
				//前缀和压缩数组
		}
	}
	for (int i = 1; i <= n; i++)
	{
		for (int j = 1; j <= n; j++)
		{
			for (int x = i; x <= n; x++)
			{
				for (int y = j; y <= n; y++)
				{
					if (a[x][y] - a[x][j - 1] - a[i - 1][y] + a[i - 1][j - 1] > ans)
					{
						ans = a[x][y] - a[x][j - 1] - a[i - 1][y] + a[i - 1][j - 1];
					}
				}
			}
		}
	}
	//for (int i = 1; i <= n; i++)//以第i行为最下面的一行
	//{
	//	for (int k = 1; k <= i; k++)//往上包含几行
	//	{
	//		int dp[1000];
	//		int s[1000];
	//		for (int j = 1; j <= n; j++)//压缩的第几列
	//		{
	//			s[j] = a[i][j] - a[i - k][j];
	//			dp[j] = max(dp[j - 1] + s[j], s[j]);
	//			ans = max(ans, dp[j]);
	//		}

	//	}
	//}
	printf("%d", ans);
	return 0;
}