算法一
O(S^2),S是障碍点的个数。
算法思想:以每个障碍点作为极大子矩形的左边界和右边界,枚举所有的极大子矩形,再找出最大的子矩形,最后还要算一下遗漏的以范围内的左边界和右边界为极大子矩形的左右边框的矩形大小。
P1578 奶牛浴场
题目描述
由于John建造了牛场围栏,激起了奶牛的愤怒,奶牛的产奶量急剧减少。为了讨好奶牛,John决定在牛场中建造一个大型浴场。但是John的奶牛有一个奇怪的习惯,每头奶牛都必须在牛场中的一个固定的位置产奶,而奶牛显然不能在浴场中产奶,于是,John希望所建造的浴场不覆盖这些产奶点。这回,他又要求助于Clevow了。你还能帮助Clevow吗?
John的牛场和规划的浴场都是矩形。浴场要完全位于牛场之内,并且浴场的轮廓要与牛场的轮廓平行或者重合。浴场不能覆盖任何产奶点,但是产奶点可以位于浴场的轮廓上。
Clevow当然希望浴场的面积尽可能大了,所以你的任务就是帮她计算浴场的最大面积。
#include <bits/stdc++.h>
using namespace std;
const int maxn = 1e5 +7;
typedef long long ll;
struct node {
int x, y;
bool operator < (node& a) const {
if(x == a.x) return y < a.y;
return x < a.x;
}
}p[maxn];
bool cmp(node a, node b) {
if(a.y == b.y) return a.x < b.x;
return a.y < b.y;
}
int main()
{
int L, W, n;
scanf("%d%d%d", &L, &W, &n);
for (int i = 1; i <= n; i++) scanf("%d%d", &p[i].x, &p[i].y);
p[++n].x = 0, p[n].y = 0;
p[++n].x = L, p[n].y = 0;
p[++n].x = 0, p[n].y = W;
p[++n].x = L, p[n].y = W;
sort(p + 1, p + 1 + n);
int ans = 0;
for (int i = 1; i <= n; i++) {//从左向右扫
int h = W, l = 0;
for (int j = i + 1; j <= n; j++) {
ans = max(ans, (p[j].x - p[i].x) * (h - l));
if(p[j].y < p[i].y) l = max(l, p[j].y);
else h = min(h, p[j].y);
}
}
for (int i = n; i >= 1; i--) {//从右向左扫
int h = W, l = 0;
for (int j = i - 1; j >= 1; j--) {
ans = max(ans, (p[i].x - p[j].x) * (h - l));
if(p[j].y < p[i].y) l = max(l, p[j].y);
else h = min(h, p[j].y);
}
}
sort(p + 1, p + 1 + n, cmp);//枚举遗漏的情况
for (int i = 1; i <= n - 1; i++) ans = max(ans, (p[i + 1].y - p[i].y) * L);
printf("%d\n", ans);
return 0;
}
算法二
O(nm)
算法思想:悬线法。预处理出每个空闲点能往上连续的空闲点的最上位置,再处理出以该点为底边的极大子矩形能扩展到左右的最远的空闲宽度,然后直接长乘宽计算即可,主要是递推的公式重要。
hdu1505
题意
给一块nm的地,F表示空闲,R表示占用,求最大的空闲的矩形面积。
#include <bits/stdc++.h>
using namespace std;
const int maxn = 1e3 + 7;
typedef long long ll;//g是原图 g[i][j] = 1表示所求的最大子矩形的组成元素
int g[maxn][maxn], u[maxn][maxn], l[maxn][maxn], r[maxn][maxn];
char s[100];//u是u[i][j]这个点的悬线的连续的1能延伸到最上面的地方
int main()//l是以g[i][j]这个点为底边的最大子矩形的左边界,r就是右边界
{
int t, n, m;
scanf("%d", &t);
while(t--) {
int ans = 0;
scanf("%d%d", &n, &m);
for (int i = 1; i <= n; i++) {
for (int j = 1; j <= m; j++) {
scanf("%s", s);
g[i][j] = s[0] == 'F' ? 1 : 0;
u[i][j] = l[i][j] = r[i][j] = 0;
if(g[i][j] == 1) u[i][j] = i, l[i][j] = r[i][j] = j;
}
}
for (int i = 1; i <= n; i++) {
for (int j = 1; j <= m; j++) {
if(g[i][j] == 0) continue;
if(i != 1) u[i][j] = g[i - 1][j] == 1 ? u[i - 1][j] : i;
if(j != 1) l[i][j] = g[i][j - 1] == 1 ? l[i][j - 1] : j;
}
for (int j = m; j >= 1; j--) {
if(g[i][j] == 0) continue;
if(j != m) r[i][j] = g[i][j + 1] == 1 ? r[i][j + 1] : j;
}
}
for (int i = 1; i <= n; i++) {
for (int j = 1; j <= m; j++) {
if(g[i][j] == 0) continue;
if(g[i][j] == 1 && g[i - 1][j] == 1) {
l[i][j] = max(l[i][j], l[i - 1][j]);
r[i][j] = min(r[i][j], r[i - 1][j]);
}
ans = max(ans, (r[i][j] - l[i][j] + 1) * (i - u[i][j] + 1));
}
}
printf("%d\n", ans*3);//ans就是最大子矩形的面积,乘3是题目要求
}
return 0;
}