对每一行进行考虑,均可以视为一个柱状图,
于是问题转化为为求柱状图的最大矩形面积,
这可以利用单调栈来处理,
利用单调栈在
O
(
m
)
O(m)
O(m) 时间内处理出所有元素左边
/
/
/ 右边第一个比其大
/
/
/ 小的元素位置。
在利用单调栈处理之前,
还需要求出 每个位置的 “柱高”,
其实只需要 对每一列进行遍历,
然后从上到下依次编号即可。
总的时间复杂度为
O
(
n
⋅
m
)
O(n⋅m)
O(n⋅m)。
代码:
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<stack>
using namespace std;
int n,m,a[2003][2003],l[100010],r[100010];
void ddz(int i)
{
stack<int>z;
for(int j=1; j<=m; j++)
{
while(!z.empty()&&a[i][j]<=a[i][z.top()])
z.pop();
if(!z.empty())
l[j]=z.top();
else
l[j]=0;
z.push(j);
}
while(!z.empty())
z.pop();
for(int j=m; j>=1; j--)
{
while(!z.empty()&&a[i][j]<=a[i][z.top()])
z.pop();
if(!z.empty())
r[j]=z.top();
else
r[j]=m+1;
z.push(j);
}
}
int main()
{
while(cin>>n>>m)
{
for(int i=1; i<=n; i++)
for(int j=1; j<=m; j++)
scanf("%d",&a[i][j]);
for(int j=1; j<=m; j++)
for(int i=1,cnt=0; i<=n; i++)
if(a[i][j]==1)
cnt++,a[i][j]=cnt;
else
cnt=0;
int ans=0;
for(int i=1; i<=n; i++)
{
ddz(i);
for(int j=1; j<=m; j++)
ans=max(ans,a[i][j]*(r[j]-l[j]-1));
}
printf("%d\n",ans);
}
return 0;
}