codevs 1159 最大全0子矩阵

 

1159 最大全0子矩阵

 时间限制: 1 s

 空间限制: 128000 KB

题目描述 Description

在一个0,1方阵中找出其中最大的全0子矩阵，所谓最大是指O的个数最多。

输入描述 Input Description

输入文件第一行为整数N，其中1<=N<=2000，为方阵的大小，紧接着N行每行均有N个0或1，相邻两数间严格用一个空格隔开。

输出描述 Output Description

输出文件仅一行包含一个整数表示要求的最大的全零子矩阵中零的个数。

样例输入 Sample Input

5
0 1 0 1 0
0 0 0 0 0
0 0 0 0 1
1 0 0 0 0
0 1 0 0 0

样例输出 Sample Output

9

这道题刚开始不会，后来看了一篇国家集训队的论文就比较简单了，是王知昆的《浅谈用极大化思想解决最大子矩形问题》，

用h[j]表示第j列从第1行到第i行的最大高度

l[j]表示h[j]的左边界，r[j]表示h[j]的右边界，最后就ans=max(ans,h[j]*(r[i] - l[i]  + 1)

代码：

#include<cstdio>
#include<iostream>
using namespace std;

int n;
int map[2000 + 5][2000 + 5];
int h[2000 + 5];
int l[2000 + 5];
int r[2000 + 5];
int ans = 0;

int main()
{
	scanf("%d",&n);
	for(int i = 1; i <= n; ++i){
		for(int j = 1; j <= n; ++j)
		scanf("%d",&map[i][j]);
	}
	for(int i = 1; i <= n; ++i)h[i] = 0,l[i] = 1,r[i] = n;
	for(int i = 1; i <= n; ++i){
		int lo = 0,ro = n + 1;
		for(int j = 1; j <= n; ++j){
			if(map[i][j]){
				h[j] = 0;
				lo = j;
				l[j] = 1;
			}
			else{
				++h[j];
				l[j] = max(l[j],lo + 1);
			}
		}
		for(int j = n; j >= 1; --j){
			if(map[i][j]){
				ro = j;
				h[j] = 0;
				r[j] = n;
			}
			else {
				r[j] = min(r[j], ro - 1);
				ans = max(ans,h[j] * (r[j] - l[j] + 1));
			}
		}
	}
	printf("%d",ans);
	return 0;
}