Type:悬线法,单调栈(未学,用悬线法做的)

提示

蓝书P51，最大子矩阵 O(mn)

Code

#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn=510;
int m,n;
int mat[maxn][maxn],up[maxn][maxn],left[maxn][maxn],right[maxn][maxn];

void print(){
    printf("Up:\n");
    for(int i=1;i<=m;++i){
        for(int j=1;j<=n;++j){
            printf("%d ",up[i][j]);
        }
        puts("");
    }
    printf("Left:\n");
    for(int i=1;i<=m;++i){
        for(int j=1;j<=n;++j){
            printf("%d ",left[i][j]);
        }
        puts("");
    }
    printf("Right:\n");
    for(int i=1;i<=m;++i){
        for(int j=1;j<=n;++j){
            printf("%d ",right[i][j]);
        }
        puts("");
    }
}

int main(){
    while(~scanf("%d%d",&m,&n)){
        for(int i=1;i<=m;++i){
            for(int j=1;j<=n;++j){
                scanf("%d",&mat[i][j]);
            }
        }
        int ans=0;
        for(int i=1;i<=m;++i){
            int lo=0,ro=n;
            for(int j=1;j<=n;++j){///从右往左扫描,维护up和left
                if(!mat[i][j]){
                    up[i][j]=left[i][j]=0;lo=j;
                }else{
                    up[i][j]=up[i-1][j]+1;
                    left[i][j]=max(left[i-1][j],lo+1);
                }
            }
            for(int j=n;j>=1;--j){///维护right
                if(!mat[i][j]){
                    right[i][j]=n+1;ro=j-1;
                }else{
                    right[i][j]=i==1?ro:min(right[i-1][j],ro);
                    ans=max(ans,up[i][j]*(right[i][j]-left[i][j]+1));
                }
            }
        }
        //print();
        printf("%d\n",ans);
    }
    return 0;
}