2019牛客多校第二场 H Second Large Rectangle 【DP+单调栈】

题目链接：https://ac.nowcoder.com/acm/contest/882/H

题目大意：给出一个01矩阵，求由1构成的第二大矩阵的面积；

/*
dp[i][j]表示当前位置向上扩展连续是1的长度
开一个结构体存每个点构成最大矩阵的宽和高，单调队列里存高递增的点，
这样栈里只有比当前点高要低的元素（因为比当前高的也用不着），然后循
环一变栈里的元素更新第一大和第二大的矩阵
*/
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<iostream>
using namespace std;
typedef long long ll;
const int N=1e3+10;

struct node
{
    int l,h;
    node(int _l=0,int _h=0):l(_l),h(_h){}
}q[N];
int dp[N][N];
int s[N][N];
char mp[N][N];

int main()
{
    int n,m;
    scanf("%d %d",&n,&m);
    for(int i=1;i<=n;i++)
    {
        scanf("%s",mp[i]+1);
        for(int j=1;j<=m;j++)
        {
            s[i][j]=mp[i][j]-'0';
            dp[i][j]=s[i][j];
            dp[i][j]+=dp[i][j]*dp[i-1][j];
        }
    }
    int max1=0,max2=0;
    for(int i=1;i<=n;i++)
    {
        int top=0;
        for(int j=1;j<=m;j++)
        {
            if(s[i][j]==0) top=0;
            else
            {
                int tmp=j;
                while(q[top].h>dp[i][j]&&top) tmp=q[top--].l;
                //退出每个大于当前高的元素
                if(q[top].h!=dp[i][j]) q[++top]=node(tmp,dp[i][j]);
                //这里如果不判断，把相等的也添加进栈，会超时
                for(int k=1;k<=top;k++)
                {
                    //更新答案
                    int are=q[k].h*(j-q[k].l+1);
                    if(are>=max1)
                    {
                        max2=max1;
                        max1=are;
                    }
                    else if(are>max2)
                        max2=are;
                }
            }
        }
    }
    printf("%d\n",max2);
    return 0;
}