【枚举】金字塔

3、金字塔

 

Jaguar 国王在一场战役大胜后决定建造一个金字塔，一方面作为纪念战争胜利的纪念碑，同时亦用作埋葬在战斗中阵亡将士们的墓地。该金字塔将建在战场的所在地，有一个a列b行的矩形底部。在金字塔的底层内将有一个较小的c列d行的矩形墓室，用来存放阵亡将士们的遗体以及他们所用过的武器。

国王的建筑师将该战场分为m列n行小方格的网格。每个小方格的高度用一个整数表示。

金字塔及其内部的墓室均应覆盖网格的完整的小方格，而且其边必须与战场的边平行。金字塔内墓室小方格的高度必须保持原有的高度不变，但金字塔底部的其余地面则需要平整，平整的方法是将高度较高的小方格内的砂土移到较低的小方格内，这样一来，金字塔底部的最后高度将是底部所覆盖的所有小方格（墓室小方格除外）高度的平均值。建筑师有权选择墓室在金字塔内的位置，但墓室四周要有至少一小方格厚度的墙围绕。

建筑师希望在战场上找到一块最佳的位置来建造该金字塔及其内部的墓室，使得对于给定大小的金字塔，其底部的高度尽可能高。

图中所示是一个战场的例子，每个小方格内的整数表示该格所处区域内的高度。灰色部份表示金字塔的底部，灰色部份所包围的白色部份则表示墓室的位置。该图给出了一个最佳位置。

 

任务

 

对于给定尺寸的战场、金字塔和墓室大小以及战场上所有小方格的高度，请编写一程序找出建造金字塔及其墓室的最佳位置，使得金字塔底部的高度为最高。

 

限制条件

 

3 ≤ m ≤ 1000

3 ≤ n ≤ 1000

3 ≤ a ≤ m

3 ≤ b ≤ n

1 ≤c≤ a – 2

1 ≤ d ≤ b – 2

所有的高度值都是1-100之间的整数。

 

输入

 

你的程序必须从文件 pyramid.in 中读入所需数据，其格式如下：

pyramid.in	描述
8 5 5 3 2 1 1 5 10 3 7 1 2 5 6 12 4 4 3 3 1 5 2 4 3 1 6 6 19 8 1 1 1 3 4 2 4 5 6 6 3 3 3 2 2 2	第一行：共有6个以空格分开的正整数，它们顺次分别为：m, n, a, b, c 及 d。 随后n行：每行有m个用空格分开的整数，每行表示一横行小方格的高度。其中，第一行表示战场小方格最上面的一行（第一行），而最后一行表示最下面的一行（第n行）。每行的m个整数依次表示该行上从第一列开始的小方格的高度。

输出

 

你的程序必须将结果数据输出到pyramid.out的文件中。

pyramid.out	描述
4 1 6 2	第一行：必须输出两个由空格分隔的整数，它们表示金字塔的底部的左上角位置，其中第一个整数代表列的坐标，而第二个整数则代表行 的坐标。 第二行：必须输出两个由空格分隔的整数，它们表示金字塔内部的墓室的左上角位置，其中第一个整数表示列 的坐标，而第二个数字则表示行 的坐标。

 

 

注：如果存在多个最优位置，那么你输出其中任意一个均被视为正确。

 

评分

 

对部份总值为30分的测试数据而言，这些测试将满足以下的条件：

3 ≤ m ≤10

3 ≤ n ≤ 10

只做到了70分。其他的超时

写的是O(n^4)的枚举

#include <cstdio>
#include <cstring>
#include <string>

long m,n,a,b,c,d;
long map[1100][1100];
long sum[1100][1100];

long getint()
{
	long rs=0;char tmp;bool sgn=1;
	do tmp = getchar();
	while (!isdigit(tmp)&&tmp-'-');
	if (tmp=='-'){sgn=0;tmp=getchar();}
	do rs=(rs<<3)+(rs<<1)+tmp-'0';
	while (isdigit(tmp=getchar()));
	return sgn?rs:-rs;
}
int main()
{
	freopen("pyramid.in","r",stdin);
	freopen("pyramid.out","w",stdout);

	m = getint();
	n = getint();
	b = getint();
	a = getint();
	d = getint();
	c = getint();

	for (long i=1;i<n+1;i++)
	{
		for (long j=1;j<m+1;j++)
		{
			map[i][j] = getint();
			sum[i][j] = sum[i][j-1]+map[i][j];
		}
	}
	for (long i=1;i<n+1;i++)
	{
		for (long j=1;j<m+1;j++)
		{
			sum[i][j] += sum[i-1][j];
		}
	}
	long ans = 0;
	long a1,a2,b1,b2 = 0;
	for (long i=1;i+a-1<n+1;i++)
	{
		for (long j=1;j+b-1<m+1;j++)
		{
			for (long k=i+1;k+c-1<i+a-1;k++)
			{
				for (long l=j+1;l+d-1<j+b-1;l++)
				{
					long tmp = sum[i+a-1][j+b-1]
						-sum[i+a-1][j-1]
						-sum[i-1][j+b-1]
						+sum[i-1][j-1];
					tmp -= sum[k+c-1][l+d-1]
						-sum[k+c-1][l-1]
						-sum[k-1][l+d-1]
						+sum[k-1][l-1];
					if (ans < tmp)
					{
						ans = tmp;
						a1 = i;
						a2 = j;
						b1 = k;
						b2 = l;
					}
				}
			}
		}
	}
	printf("%ld %ld\n%ld %ld\n",a2,a1,b2,b1);
	return 0;
}