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;
}