POJ3493 Largest Submatrix of All 1’s（单调栈）

题目给一个01矩阵，求最大的1子矩阵。

先用dp预处理出每一行的每一列的1能向上按连续的1延伸多少，然后枚举每一行作为子矩阵的底，那样对于每一行的答案就是POJ2559这个经典问题了。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define MAXN 2222
int stack[MAXN],top,l[MAXN],r[MAXN];
int calc(int *a,int n){
    a[++n]=-1; top=0;
    for(int i=1; i<=n; ++i){
        l[i]=r[i]=i;
        while(top && a[stack[top]]>a[i]){
            r[stack[top]]=i-1;
            l[i]=l[stack[top]];
            --top;
        }
        if(top && a[stack[top]]==a[i]) l[i]=l[stack[top]];
        stack[++top]=i;
    }
    int res=0;
    for(int i=1; i<n; ++i){
        res=max(res,a[i]*(r[i]-l[i]+1));
    }
    return res;
}
int d[MAXN][MAXN];
int main(){
    int n,m;
    while(~scanf("%d%d",&n,&m)){
        for(int i=1; i<=n; ++i){
            for(int j=1; j<=m; ++j) scanf("%d",&d[i][j]);
        }
        for(int i=1; i<=n; ++i){
            for(int j=1; j<=m; ++j){
                if(d[i][j]) d[i][j]=d[i-1][j]+1;
            }
        }
        int res=0;
        for(int i=1; i<=n; ++i){
            res=max(res,calc(&d[i][0],m));
        }
        printf("%d
",res);
    }
    return 0;
}