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

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

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

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<algorithm>
 4 using namespace std;
 5 #define MAXN 2222
 6 int stack[MAXN],top,l[MAXN],r[MAXN];
 7 int calc(int *a,int n){
 8     a[++n]=-1; top=0;
 9     for(int i=1; i<=n; ++i){
10         l[i]=r[i]=i;
11         while(top && a[stack[top]]>a[i]){
12             r[stack[top]]=i-1;
13             l[i]=l[stack[top]];
14             --top;
15         }
16         if(top && a[stack[top]]==a[i]) l[i]=l[stack[top]];
17         stack[++top]=i;
18     }
19     int res=0;
20     for(int i=1; i<n; ++i){
21         res=max(res,a[i]*(r[i]-l[i]+1));
22     }
23     return res;
24 }
25 int d[MAXN][MAXN];
26 int main(){
27     int n,m;
28     while(~scanf("%d%d",&n,&m)){
29         for(int i=1; i<=n; ++i){
30             for(int j=1; j<=m; ++j) scanf("%d",&d[i][j]);
31         }
32         for(int i=1; i<=n; ++i){
33             for(int j=1; j<=m; ++j){
34                 if(d[i][j]) d[i][j]=d[i-1][j]+1;
35             }
36         }
37         int res=0;
38         for(int i=1; i<=n; ++i){
39             res=max(res,calc(&d[i][0],m));
40         }
41         printf("%d\n",res);
42     }
43     return 0;
44 }