poj 3494 Largest Submatrix of All 1’s

  这道题目大意是求一个0-1矩阵的最大全1子矩阵的元素个数，开始想着动态规划什么的，结果一筹莫展，看来某大神题解，终于明白原来是用单调栈的，思路大概这样：可以先计算出每个元素向上的最大伸展值，在利用单调栈分别计算该元素向左和向右最大伸展到的位置，最后遍历每个元素，计算得最大面积，也就是最大元素个数。下面贴上代码：

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <stack> 
#define maxn 2000+5

using namespace std;
typedef pair<int,int> p;
int l[maxn][maxn],r[maxn][maxn],h[maxn][maxn],a[maxn][maxn];

int main()
{
//	freopen("input.txt","r",stdin);
	int n,m,i,j;
	while(scanf("%d%d",&n,&m)!=EOF){
		for(i=1;i<=n;i++)
		{
			for(j=1;j<=m;j++)
			{
				scanf("%d",&a[i][j]);
			}
		}
		memset(h,0,sizeof(h));
		for(j=1;j<=m;j++)
		{
			for(i=1;i<=n;i++)
			{
				if(a[i][j]) h[i][j]=h[i-1][j]+1;
			}
		}
		stack<p> s;
		memset(l,0,sizeof(l));
		memset(r,0,sizeof(r));
		for(i=1;i<=n;i++)
		{
			while(!s.empty()) s.pop();
			for(j=1;j<=m;j++)
			{
				while(!s.empty()&&s.top().first>=h[i][j]) s.pop();
				if(s.empty()) l[i][j]=1;
				else l[i][j]=s.top().second+1;
				s.push(p(h[i][j],j));
			}
		} 
		for(i=1;i<=n;i++)
		{
			while(!s.empty()) s.pop();
			for(j=m;j>=1;j--)
			{
				while(!s.empty()&&s.top().first>=h[i][j]) s.pop();
				if(s.empty()) r[i][j]=m;
				else r[i][j]=s.top().second-1;
				s.push(p(h[i][j],j));
			}
		}
		int ans=0;
		for(i=1;i<=n;i++)
		{
			for(j=1;j<=m;j++)
			{
				ans=max(ans,h[i][j]*(r[i][j]-l[i][j]+1));
			}
		}
		printf("%d\n",ans);
	}
	return 0;
}