https://ac.nowcoder.com/acm/contest/882/H
题意:给定一个
n
×
m
n×m
n×m的01矩阵(
1
≤
n
,
m
≤
1000
1≤n,m≤1000
1≤n,m≤1000),求面积第二大的全1矩阵的面积。
求最大全1矩阵面积的方法参考这个博客http://www.luyixian.cn/news_show_102909.aspx
对于每一行的第
j
j
j个位置,用
h
[
j
]
h[j]
h[j]表示从该位置向上连续1的数量,用
l
[
j
]
l[j]
l[j]表示从该位置向左连续1的最左位置,用
r
[
j
]
r[j]
r[j]表示从该位置向右连续1的最左位置,这样最大全1矩阵面积就是
m
a
x
(
h
[
j
]
∗
(
r
[
j
]
−
l
[
j
]
+
1
)
)
max(h[j]*(r[j]-l[j]+1))
max(h[j]∗(r[j]−l[j]+1))。
求
l
[
j
]
l[j]
l[j]和
r
[
j
]
r[j]
r[j]的方法是,先令其初始值为
j
j
j,如果
h
[
j
]
≤
h
[
l
[
j
]
−
1
]
h[j]≤h[l[j] - 1]
h[j]≤h[l[j]−1] (用图形来说就是j左边的矩形向上扩展的大,那么左边的矩阵就可以并入当前位置),那么
l
[
j
]
=
l
[
l
[
j
]
−
1
]
l[j]=l[l[j]-1]
l[j]=l[l[j]−1],同理可得:如果
h
[
j
]
≤
h
[
r
[
j
]
+
1
]
h[j]≤h[r[j]+1]
h[j]≤h[r[j]+1],那么
r
[
j
]
=
r
[
r
[
j
]
+
1
]
r[j]=r[r[j]+1]
r[j]=r[r[j]+1]。
这个题要求第二大的面积,就要维护前二大的面积,维护过程中要注意去重,去重的方法就是根据左右端点(
l
[
j
]
l[j]
l[j]和
r
[
j
]
r[j]
r[j])判断是否是同一个矩形。另外因为通过该方法求出来的是最大面积的矩阵,而该题要求第二大面积,所以在维护时不仅要插入
h
[
j
]
∗
(
r
[
j
]
−
l
[
j
]
+
1
)
h[j]*(r[j]-l[j]+1)
h[j]∗(r[j]−l[j]+1)这个面积,而且要插入长-1和宽-1的面积。
#include<iostream>
#include<cstdio>
#include<set>
using namespace std;
int n,m,ans1,ans2;
int h[2010],r[2010],l[2010];
char a[2010][2010];
struct node{
int lx,rx;
node(int _lx,int _rx)
{
lx=_lx;rx=_rx;
}
bool operator<(const node &b) const
{
return lx<b.lx||(lx==b.lx&&rx<b.rx);
}
};
set<node> s;
void maxans(int x)
{
if(x>=ans1)
{
ans2=ans1;
ans1=x;
}
else if(x>ans2)
{
ans2=x;
}
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
scanf("%s",a[i]+1);
for(int i=1;i<=n;i++)
{
s.clear();
for(int j=1;j<=m;j++)
if(a[i][j]=='1') h[j]++;
else h[j]=0;
for(int j=1;j<=m;j++)
{
l[j]=j;
while(l[j]>1&&h[l[j]-1]>=h[j])
l[j]=l[l[j]-1];
}
for(int j=m;j>=1;j--)
{
r[j]=j;
while(r[j]<m&&h[r[j]+1]>=h[j])
r[j]=r[r[j]+1];
}
for(int j=1;j<=m;j++)
{
//printf("%d:%d %d %d ",j,l[j],r[j],h[j]);
if(s.find(node(l[j],r[j]))==s.end())
{
s.insert(node(l[j],r[j]));
maxans(h[j]*(r[j]-l[j]+1));
maxans(h[j]*(r[j]-l[j]));
maxans((h[j]-1)*(r[j]-l[j]+1));
}
}
}
printf("%d",ans2);
return 0;
}