九度1497：面积最大的全1子矩阵 (单调队列，单调栈)

题目描述：

在一个M * N的矩阵中，所有的元素只有0和1，从这个矩阵中找出一个面积最大的全1子矩阵，所谓最大是指元素1的个数最多。

输入：

输入可能包含多个测试样例。

对于每个测试案例，输入的第一行是两个整数m、n(1<=m、n<=1000)：代表将要输入的矩阵的大小。

矩阵共有m行，每行有n个整数，分别是0或1，相邻两数之间严格用一个空格隔开。

输出：

对应每个测试案例，输出矩阵中面积最大的全1子矩阵的元素个数。

样例输入：

2 2
0 0
0 0
4 4
0 0 0 0
0 1 1 0
0 1 1 0
0 0 0 0

样例输出：

0

4

思路：和之前zjnu1735差不多，只不过这里求的是最大面积，用q[ ][0]表示高度，q[ ][1]表示下标就行了。

#include<iostream>
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<vector>
#include<map>
#include<set>
#include<queue>
#include<stack>
#include<string>
#include<bitset>
#include<algorithm>
using namespace std;
typedef long long ll;
typedef long double ldb;
#define inf 99999999
#define pi acos(-1.0)
#define maxn 1005
int gra[maxn][maxn];
int h[maxn][maxn],l[maxn],r[maxn];
int q[111111][2];

int main()
{
int n,m,i,j;
while(scanf("%d%d",&n,&m)!=EOF)
{
memset(h,0,sizeof(h));
for(i=1;i<=n;i++){
for(j=1;j<=m;j++){
scanf("%d",&gra[i][j]);
if(gra[i][j]==0)h[i][j]=0;
else{
h[i][j]=1;
if(i>1){
h[i][j]+=h[i-1][j];
}
}
}
}

int front,rear;
int maxx=0;
for(i=1;i<=n;i++){
front=1;rear=0;
for(j=1;j<=m;j++){
while(front<=rear && q[rear][0]>=h[i][j]){
rear--;
}
if(rear==0){
l[j]=1;
}
else{
l[j]=q[rear][1]+1;
}
rear++;
q[rear][0]=h[i][j];
q[rear][1]=j;
}
front=1;rear=0;
for(j=m;j>=1;j--){
while(front<=rear && q[rear][0]>=h[i][j]){
rear--;
}
if(rear==0){
r[j]=m;
}
else{
r[j]=q[rear][1]-1;
}
maxx=max(maxx,h[i][j]*(r[j]-l[j]+1) );
rear++;
q[rear][0]=h[i][j];
q[rear][1]=j;
}
}
printf("%d\n",maxx);
}
return 0;
}