1158 全是1的最大子矩阵

1158 全是1的最大子矩阵

基准时间限制：1 秒 空间限制：131072 KB

给出1个M*N的矩阵M1，里面的元素只有0或1，找出M1的一个子矩阵M2，M2中的元素只有1，并且M2的面积是最大的。输出M2的面积。

Input

第1行:2个数m,n中间用空格分隔(2 <= m,n <= 500)
第2 - N + 1行：每行m个数，中间用空格分隔，均为0或1。

Output

输出最大全是1的子矩阵的面积。

Input示例

3 3
1 1 0
1 1 1
0 1 1

Output示例

4
思路：单调栈；
先预处理出每行中连续1的个数，然后看每列，像这样用单调栈维护当前点能扩展到的最大范围，然后面积是（R[i]-L[i]）*ans[i];
然后更新最大值；复杂度O（n*m）；

#include<stdio.h>
#include<algorithm>
#include<iostream>
#include<string.h>
#include<stdlib.h>
#include<queue>
#include<set>
#include<vector>
#include<map>
#include<stack>
using namespace std;
typedef long long LL;
typedef struct node
{
        int x;
        int id;
} ss;
stack<ss>sta;
int a[600][600];
int b[600][600];
int L[600];
int R[600];
int main(void)
{
        int n,m;
        while(scanf("%d %d",&n,&m)!=EOF)
        {
                int i,j;
                for(i = 1; i <= n; i++)
                {
                        for(j = 1; j <= m; j++)
                        {
                                scanf("%d",&a[i][j]);
                        }
                }
                memset(b,0,sizeof(b));
                for(i = 1; i <= n; i++)
                {
                        for(j = 1; j <= m; j++)
                        {
                                if(a[i][j])
                                {
                                        b[i][j] = b[i][j-1]+1;
                                }
                                //printf("%d ",b[i][j]);
                        }
                }
                int maxx = 0;
                for(j = 1; j <= m; j++)
                {
                        while(!sta.empty())
                                sta.pop();
                        for(i = 1; i <= n; i++)
                        {
                                while(true)
                                {
                                        ss ak ;
                                        ak.x = b[i][j];
                                        ak.id = i;
                                        if(sta.empty())
                                        {
                                                L[i] = 1;
                                                sta.push(ak);
                                                break;
                                        }
                                        else
                                        {
                                                ss ac = sta.top();
                                                if(ac.x > b[i][j])
                                                {
                                                        sta.pop();
                                                        R[ac.id] = i-1;
                                                }
                                                else
                                                {
                                                        L[i] = ac.id;
                                                        sta.push(ac);
                                                        break;
                                                }
                                        }
                                }
                        }
                        while(!sta.empty())
                        {
                            ss al = sta.top();
                            sta.pop();
                            R[al.id] = n;
                        }
                        for(i = 1; i <= m; i++)
                        {
                                int xx = R[i]-L[i]+1;
                                maxx = max(maxx,xx*b[i][j]);
                        }
                }
                printf("%d
",maxx);
        }
        return 0;
}