BZOJ3726: PA2014Final Wykladzina

先不考虑可以包含坏格
如果求不包含坏格的最大子矩形，有一种办法是枚举矩形的下界，用一个单调队列扫过去，队列里维护左端点和对应的高度，显然随左端点递增，高度递增，预处理u[i][j]代表(i,j)上方第一个坏格的位置，枚举右端点j扫过去时维护队列里的高度，如果当前队尾不合法，那他在j-1的时候合法且面积最大，更新一下答案，更新一下队列

那么包含一个坏格的最大子矩形就多维护一个队列就行了，每次从不包含坏格队列里出来的子矩形可以把他扔进第二个队列，用u[u[i][j]][j]更新一下它的高度

细节有点多qwq，建议拍一下再交

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

const int maxn = 2100;

int n,m;
char str[maxn];
bool v[maxn][maxn];

int L[maxn][maxn],U[maxn][maxn];

struct node{int p,u;}q[maxn],q2[maxn],t[maxn]; int tp;
int head,tail,head2,tail2;

int re;
int cal(int l,int r,int d,int u) { return (r-l+1)*(u-d+1); }

int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
    {
        scanf("%s",str+1);
        for(int j=1;j<=m;j++) if(str[j]=='#') v[i][j]=true;
    }
    for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)
    {
        L[i][j]=L[i][j-1];
        if(v[i][j]) L[i][j]=j;
    }
    for(int j=1;j<=m;j++) for(int i=1;i<=n;i++)
    {
        U[i][j]=U[i-1][j];
        if(v[i][j]) U[i][j]=i;
    }

    for(int i=1;i<=n;i++)
    {
        head=1,tail=0;
        head2=1,tail2=0;
        for(int j=1;j<=m;j++)
        {
            int u=U[i][j],uu=0;
            if(u) uu=U[u-1][j];
            int la=0;
            while(head2<=tail2&&q2[tail2].u<=u) 
            {
                re=max(re,cal(q2[tail2].p,j-1,q2[tail2].u,i));
                la=q2[tail2--].p;
            }
            if(!(head2<=tail2&&q[tail2].u==u+1)&&la) q2[++tail2]=(node){la,u+1};

            la=0; int la2=0;
            while(head<=tail&&q[tail].u<=u)
            {
                re=max(re,cal(q[tail].p,j-1,q[tail].u,i));
                if(q[tail].u>uu) t[++tp]=q[tail];
                else la2=q[tail].p;
                la=q[tail--].p;
            }
            if(!(head<=tail&&q[tail].u==u+1)&&la) q[++tail]=(node){la,u+1};
            else q[++tail]=(node){j,u+1};

            while(tp) 
            {
                if(!la2||t[tp].p<la2) 
                {
                    while(head2<=tail2&&q2[tail2].p>=t[tp].p) tail2--;
                    q2[++tail2]=t[tp];
                }
                tp--;
            }
            if(!(head2<=tail2&&q2[tail2].u==uu+1)&&la2) 
            {
                while(head2<=tail2&&q2[tail2].p>=la2) tail2--;
                q2[++tail2]=(node){la2,uu+1};
            }
            else q2[++tail2]=(node){j,uu+1};
        }
        while(head<=tail) re=max(re,cal(q[head].p,m,q[head].u,i)),head++;
        while(head2<=tail2) re=max(re,cal(q2[head2].p,m,q2[head2].u,i)),head2++;
    }
    printf("%d\n",re);

    return 0;
}