question （单调栈）

10.9

首先我们要预处理以(I,j)为底的一个全1列最多有多高记作a[i][j]，那么如样例
0 0 0 0
0 1 1 0
0 1 1 0
0 0 0 0
就被转成
0 0 0 0
0 1 1 0
0 2 2 0
0 0 0 0
然后我们枚举以某一行为底时的最大矩阵。
对于一个位置（i，j），如果我们以所在行为底，高度为a[i][j]这个矩阵最大为多少？
这时我们要用到一种叫做单调栈的数据结构，原理如下：我们想要知道这个矩阵向左最多能延伸多少，需要维护一个单调增加的栈，例如1 5 3 4 2 9
1时，栈内为空，1进栈，
5时，栈内有1,5进栈不影响性质，
3时，栈内为1 5，先讲5出栈，再把3入栈。
这样我们就可以O(n)得到n个数据最靠近的第一个大于他的数，同理就可以得到向右能延伸到的最远的位置。

代码可能和题解在细节上有一些差别，但方法是一样的。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
#define LL long long 
#define N 1010
using namespace std;

int n, m, ans;
int mp[N][N], up[N][N], stk[N], ls[N];

void solve(int *cc){
    int tp = 0;
    for(int i=1; i<=m; ++i){
        if(cc[i] >= stk[tp]){//维护一个单调栈 
            stk[++tp] = cc[i];
            ls[tp] = 1;//以tp作为当前矩形的高向左做多可以覆盖的位置 
        }
        else{
            int len = 0;
            while(tp && stk[tp] > cc[i]){
                len += ls[tp];
                ans = max(ans, len * stk[tp--]);
            }
            stk[++tp] = cc[i];
            ls[tp] = len + 1;
        }
    }
    int len = 0;
    while( tp ){//处理剩下的 
        len += ls[tp];
        ans = max(ans, stk[tp--] * len);
    }
}

int main(){
    freopen( "question.in", "r", stdin );
    freopen( "question.out", "w", stdout );
    cin >> n >> m;
    for(int i=1; i<=n; i++)
        for(int j=1; j<=m; j++){
            scanf("%d", &mp[i][j]);
            if(mp[i][j]) up[i][j] = up[i-1][j] + 1;//从当前点出发向上连续的1的高度 
        }
    for(int i=1; i<=n; ++i) solve( up[i] );
    cout << ans << endl;
    return 0;
}