Second Large Rectangle
最大矩形问题
题目描述
给出一个N×M二元01矩阵。求由1组成的第二大矩形的大小。
输入描述
输入的第一行包含两个空格分隔的整数N和M,后面的N行每行包含M个字符c i j _{ij} ij
1 ≤ \leq ≤ N ≤ \leq ≤ 14
M × \times × N ≥ \geq ≥ 2
c i j _{ij} ij ∈ \in ∈ “01”
输出描述
输出一行,用一个整数表示答案。如果由1组成的矩形少于2个,则输出0
输入样例
1 2
01
1 3
101
输出样例
0
1
解题思路
解决这题首先我们要知道如何求最大矩形
step1求高度矩阵
在N × \times ×M的01矩阵中,先求出给出矩阵的高度矩阵(对于矩阵中的一格c i j _{ij} ij,其高度为从该格往上(i方向)数连续的1的个数)。
时间O(n 2 ^2 2)。
step2求每一格所能扩展出的最大矩阵
将高度矩阵中的每一行取出,用其j做x轴,高度h做y轴作出对应的高度直方图,默认其宽度为1。那么对于每一格,按照其高度向左右两边扩展,直到左右两边的高度都小于它(再往前扩展就无法构成矩形了)或者到达边界,最终得到该格对应的最大矩形的左右边界l,r。这样该格对应的最大矩形的面积可以表示为该格的高度h乘以该格的宽度[l,r]。
这里可以用一个单调栈以O(n) 的时间遍历每一行(每个元素只入栈出栈一次)总共有n行,所以这一步的时间是 O(n 2 ^2 2)
step3比较每格的矩形面积求出最大矩形
时间是 O(n 2 ^2 2)
step4次大矩形
这里在最大矩阵的基础上稍加处理就行了。
对于每一格,其对应的次大矩形可能是其对应的最大矩形减少一个高度或者减少一个宽度,比较的时候加上就可以了。
另外还有一个细节问题,如果最大和次大的面积相同,一定要保证它们是不同的矩形。
如果两个矩形的面积和左右边界都相同,那么我们可以认为它们是同一个矩形。
举个例子:
101
这里的两个1 × \times × 1矩形是不同的,其左右边界不同
111
这里有三个相同的矩形,c 00 _{00} 00,c 01 _{01} 01,c 02 _{02} 02,它们对应的最大矩形都是1 × \times × 3 的111。如果其中一个做了最大矩形,那么剩下的都不能作为次大矩形,所以在实际计算时需要去重。
AC代码
#include<iostream>
#include<stack>
#include<vector>
#include<algorithm>
#include<cstring>
using namespace std;
const int N = 1005;
int n,m;
int map[N][N];
int h[N][N];
int l[N],r[N];
string str;
int s1,s2;//第一大,第二大
struct Node{//用来存矩形的面积,和它在直方图中的左右边界
int area;
int left;
int right;
Node(int a,int l,int r){
area=a;
left=l;
right=r;
}
};
bool cmp(Node a,Node b){
return a.area > b.area;
}
bool isEqual(Node a,Node b){//判断两个矩形是否相等
if(a.area==b.area&&a.left == b.left && a.right == b.right )
return true;
return false;
}
void solve(){
vector<Node> ss;//存面积最大的两个矩阵
ss.push_back(Node(0,0,0));
ss.push_back(Node(0,0,0));
//初始化直方图高度矩阵
int height=0;
for(int j=0;j<m;j++){
height = 0;
for(int i=0;i<n;i++){
if(map[i][j]==1)
h[i][j] = ++height;
else height = 0;
}
}
int s[3];
for(int i=0;i<n;i++){
//对于每个格子,找出其拓展出的最大矩阵的l,r
stack<int> L;
stack<int> R;
for(int j=0;j<m;j++){
if(j==0){
L.push(j);
l[j] = j;
}
else{//注意empty()要在前面判断,否则当栈为空时访问top会出错
while(!L.empty()&&h[i][j]<=h[i][L.top()]){//当前高度比左边小,可以往左边扩展
L.pop();
}
if(L.empty()) l[j] = 0;
else l[j] = L.top()+1;
//此时高度比左边大,压栈,维持升序栈
L.push(j);
}
}
for(int j=m-1;j>=0;j--){
if(j==m-1){
R.push(j);
r[j]=j;
}else{
while(!R.empty() && h[i][j]<=h[i][R.top()] ){//当前高度比右边小,可以往右边扩展
R.pop();
}
if(R.empty()) r[j] = m-1;
else r[j] = max(0,R.top()-1);
//此时高度比左边大,压栈,维持升序栈
R.push(j);
}
}
for(int j=0;j<m;j++){//当前格扩展出的最大矩形的面积
s[0] = h[i][j]*(r[j]-l[j]+1);
s[1] = (h[i][j]-1)*(r[j]-l[j]+1);//高度减一
s[2] = h[i][j]*(r[j]-l[j]);//宽度减一
ss.push_back(Node(s[0],l[j],r[j]));
ss.push_back(Node(s[1],l[j],r[j]));
ss.push_back(Node(s[2],l[j],r[j]-1));
//排序
sort(ss.begin(),ss.end(),cmp);
//去除重复
int kk=0;
for(int k=0;k<ss.size();k++)
if(!isEqual(ss[k],ss[k-1]))
ss[kk++]=ss[k];
//只保留前两个
ss.erase(ss.begin()+2,ss.end());
}
}
cout<<ss[1].area<<endl;
}
int main(){
ios::sync_with_stdio(false);
while(cin>>n>>m){
for(int i=0;i<n;i++){
cin>>str;
for(int j=0;j<m;j++)
map[i][j]=str[j]-'0';
}
solve();
}
return 0;
}
总结
这题所用到的知识:
最大矩形问题
单调栈
排序去重