P4147玉蟾宫——最大子矩阵

悬线法裸题。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int const MAXN=,MAXM=;
int n,m,lf[MAXN][MAXM],rt[MAXN][MAXM],ans,a[MAXN][MAXM],sta[MAXM],top;
bool d[MAXN][MAXM];
int main()
{
    scanf("%d%d",&n,&m);
    char dc;
    for(int i=;i<=n;i++)
        for(int j=;j<=m;j++)
        {
            cin>>dc;
            if(dc=='R')d[i][j]=;
            else
            {
                for(int k=i-;k;k--)
                    if(d[k][j])
                    {
                        a[i][j]=i-k;
                        break;
                    }
                if(!a[i][j])a[i][j]=i;
//                printf("a[%d][%d]=%d\n",i,j,a[i][j]);
            }
        }
    for(int i=;i<=n;i++)
    {
        memset(sta,,sizeof sta);
        top=;
        for(int j=;j<=m;j++)
        {
            while(a[i][j]<=a[i][sta[top]]&&top)top--;
            lf[i][j]=sta[top];
            sta[++top]=j;
        }
        memset(sta,,sizeof sta);
        top=;sta[top]=m+;
        for(int j=m;j;j--)
        {
            while(a[i][j]<=a[i][sta[top]]&&top)top--;
            rt[i][j]=sta[top];
            sta[++top]=j;
        }
//        for(int j=1;j<=m;j++)
//            printf("lf[%d][%d]=%d rt[%d][%d]=%d\n",i,j,lf[i][j],i,j,rt[i][j]);
    }
    for(int i=;i<=n;i++)
        for(int j=;j<=m;j++)
        {
            int l,r;
            if(d[i][j])continue;
            if(!d[i][j]&&i>)
            {
                l=max(lf[i][j],lf[i-][j]);
                r=min(rt[i][j],rt[i-][j]);
                ans=max(ans,a[i][j]*(r-l-));
            }
            else ans=max(ans,rt[i][j]-lf[i][j]-);
        }
    printf("%d",*ans);
    return ;
}