Second Large Rectangle（牛客多校第二场）

Second Large Rectangle

本题大意为给一个N*M由01构成的矩形，求其中的第二大的全一矩阵的面积。

可以想到，第二大矩形的面积肯定比第一大矩形的小，所以我们可以找到每一个极大全1矩形，然后用h和l代表每个全1矩阵的长和高，把hl，h(l-1),(h-1)*l存进一个规格为2的优先队列，然后把第一个弹出输出就可以了（本弱鸡垃圾做法，不过容易理解），重点就是怎样找到每一个极大全1矩阵，这里需要用到单调栈的知识。

对序列中每个元素 i ，以它为最小值向两边扩展，当前的贡献值为扩展长度 len 乘 元素值 val, 求出最大的贡献值。
如果要一直扩展下去，必须当前元素值不小于扩展元素值。
如果我们从第 1 个元素计算到第 n 个元素，计算到元素 i 时，
如果我们能快速的找出左边第一个比 i 小的元素，就可以完成向左的扩展，
而右扩展，这个过程本来就是一直在向右走，所以直接扩展即可。
用单调递减栈来维护这个过程，栈内储存的为坐标i。
在每个元素入栈的时候，更新入栈元素向左扩展的最小位置
在每个元素出栈的时候，更新栈内元素向右扩展的最大位置
（引自https://www.jianshu.com/p/fde143beb0ee）

#include <iostream>
#include <cstdio>
#include <cstring>
#include <vector>
#include <algorithm>
#include <map>
#include <string>
#include <queue>
#include <stack>
#define inf 0x3f3f3f3f
using namespace std;
typedef long long ll;
const int maxn=2005;
int n,m,h[maxn][maxn],l[maxn],r[maxn];
map<int,map<int,int>>flag;
char a[maxn][maxn];
stack<int>s;
priority_queue<int,vector<int>,greater<int> >q;//存面积的优先队列
int main()
{
    scanf("%d%d",&n,&m);
    memset(h,0,sizeof(h));
    for(int i=1; i<=n; i++)
        for(int j=1; j<=m; j++)
        {
            scanf(" %c",&a[i][j]);
            if(a[i][j]=='0')h[i][j]=0;
            else h[i][j]=h[i-1][j]+1;//h[i][j]表示从（i，j）这个点往上有多少连续的1
        }
    for(int p=1; p<=n; p++)//每一行单独算，复杂度有点高的,n,m<=1000可以用
    {
        for(int i=1; i<=m; i++)
        {
            while(!s.empty()&&h[p][i]<h[p][s.top()])//找到最近的一个比当前点的高度低的点
            {
                r[s.top()]=i-1;//当前点比栈中的点小的话就把战中的点弹出并且把栈中点的r标为当前点下标减一
                s.pop();
            }
            if(s.empty())l[i]=1;
            else if(h[p][s.top()]==h[p][i])l[i]=l[s.top()];//与当前点相同则继承一下
            else l[i]=s.top()+1;//比当前点低则加一记为矩形的l
            s.push(i);//每一次把当前点存到栈里
        }
        while(!s.empty())
        {
            r[s.top()]=m;//最后栈中剩下的全部的r标为最右边
            s.pop();
        }
        for(int i=1; i<=m; i++)
            if(!flag[l[i]][r[i]])//确保不会有重复的，比如2和3扩散的矩形为一个
            {
                flag[l[i]][r[i]]=1;
                int p1=h[p][i],p2=(r[i]-l[i]+1);
                q.push(p1*p2);
                q.push(p1*(p2-1));
                q.push((p1-1)*p2);
            }
        while(q.size()>2)q.pop();
        flag.clear();
    }
    printf("%d\n",q.top());
    while(!q.empty())q.pop();
    return 0;
}