POJ 3494 Largest Submatrix of All 1’s

意思和做法和HDU的1505可以说完全一样了。

这里有个问题。你可能会发现我这个代码，若某段连续区间是等高的，则除了这段区间最左边一个，其他的求出来的左区间的界限是错的，那为什么这段代码还能AC呢，

其实是这样的，对于这段连续区间来说若是区间界限都求正确的话，那么对于这个区间的任意一点来说，作为中间点的值都是一样的，对于以下的求法来说，只要最左边

的点是求正确的，其他点都是错的，就并不影响最终答案的正确性。

#include <cstdio>
#include <cstring>
#include <string>
#include <iostream>
#include <map>
#include <vector>
#include <cmath>
#include <stack>
#include <queue>
#include <cstdlib>
#include <algorithm>
using namespace std;
typedef __int64 int64;
typedef long long ll;
#define M 50005
#define N 1000005
#define max_inf 0x7f7f7f7f
#define min_inf 0x80808080
#define mod 1000000007

int m , n , top;
int q[2005] , lef[2005] , rig[2005];

int main()
{
	int h[2005];
	int ans;
	while (~scanf("%d%d",&m,&n))
	{
		int i , j;
		ans = 0;
		memset(h , 0 , sizeof h);
		for (i = 0 ; i < m ; i++)
		{
			top = 0;
			for (j = 0 ; j < n ; j++)
			{
				int x;
				scanf("%d",&x);
				h[j]++;
				h[j] = x ? h[j] : 0;
				while (top > 0 && h[q[top-1]] > h[j])
				{
					rig[q[top-1]] = j-1;
					top--;
				}
				if (top <= 0)lef[j] = 0;
				else lef[j] = q[top-1]+1;
				q[top++] = j;
			}
			while (top)rig[q[top-1]] = n-1,top--;

			for (j = 0 ; j < n ; j++)
			{
				int temp = (rig[j]-lef[j]+1)*h[j];
				ans = max(ans,temp);
			}
		}
		printf("%d\n",ans);
	}
	return 0;
}