poj 3494 Largest Submatrix of All 1’s 单调栈/DP迭代法

题意：给出一个m*n的01矩阵，问在系数全为1的子矩阵中 系数和的最大值。

解法：DP中的迭代法，或者单调栈。

不论哪种方法核心都是对于每个高度为h的矩形，找到最左高度大于等于它的位置，和最右的位置。

#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<vector>
using namespace std;

#define all(x) (x).begin(), (x).end()
#define for0(a, n) for (int (a) = 0; (a) < (n); (a)++)
#define for1(a, n) for (int (a) = 1; (a) <= (n); (a)++)
typedef long long ll;
typedef pair<int, int> pii;
const int INF =0x3f3f3f3f;
const int maxn= 2000   ;


struct Node
{
    int h;
    int le;
};
int n,m,ans;
int h[maxn+10];

Node sta[maxn+10];
int main()
{
   while(~scanf("%d%d",&m,&n))
   {
       int x,ans=0;
       memset(h,0,(n+1)*sizeof h[0]);
       for1(i,m)
       {
           int top=0;
           for1(j,n)
          {
             scanf("%d",&x);
             if(!x) h[j]=0;
             else h[j]++;

             if(top&&sta[top].h>h[j])
             {
                while(top&&sta[top].h>h[j])
                {
                    ans=max(ans,  sta[top].h*(j-sta[top].le) );
                    top--;
                }
                sta[++top].h=h[j];
             }
             else
             {
                 sta[++top].h=h[j];
                 sta[top].le=j;
                if(top>=2&&sta[top-1].h==sta[top].h ) sta[top].le=sta[top-1].le;
             }
          }
          while(top)
          {
             ans=max(ans,sta[top].h*( n+1- sta[top].le) );
              top--;
          }
       }
        printf("%d\n",ans);

   }

   return 0;
}

/*
4 4
1 0 0 0
1 1 1 1
0 1 1 1
1 0 0 1

*/