POJ3494 Largest Submatrix of All 1’s【单调栈】

对每一行进行考虑，均可以视为一个柱状图，
于是问题转化为为求柱状图的最大矩形面积，
这可以利用单调栈来处理，
利用单调栈在 $O(m)$ 时间内处理出所有元素左边 $/$ 右边第一个比其大 $/$ 小的元素位置。
在利用单调栈处理之前，
还需要求出 每个位置的 “柱高”，
其实只需要 对每一列进行遍历，
然后从上到下依次编号即可。
总的时间复杂度为 $O(n\cdot 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;
}