【PA2014Final】【BZOJ3726】Wykladzina

最新推荐文章于 2021-03-11 10:08:21 发布

原创最新推荐文章于 2021-03-11 10:08:21 发布 · 1.1k 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#乱搞

随便搞搞专栏收录该内容

403 篇文章

订阅专栏

本文介绍了一种在地毯中寻找面积最大的完美子矩形的方法，该子矩形内至多包含一个缺陷格子。通过算法分析，实现复杂度优化，确保效率和准确性。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

Description

有一块w*s的地毯，其中有些格子是坏的。从中剪出一块面积尽可能大的子矩形，使得其中至多包含一个坏格。

Input

第一行两个整数w,s(1<=w,s<=2000)。
接下来w行，每行有s个字符，’.’表示这格是好的，’#’表示这格是坏的。

Output

输出最大面积。

Sample Input
Sample Output

12
HINT

Source

鸣谢Jcvb

感觉每次写这种题多半会变成通过自己的blog免费给Claris的题解打广告.
一开始感觉此题可以悬线法?但是悬线法从来没学过,不会用.
所以我们还是走传送门吧.传送至Claris题解
主要的问题在于此题的限制是:最多一个坏点,所以考虑对有一个坏点的边界,没有坏点的边界都找出来,然后我们可以使用有一个坏点的边界和没有坏点的边界互相替换推理一下,最后得到答案.详细内容看Claris题解.
细节上并不是很好理解…需要认真地体会一下..

//连代码都变成Claris风的辣鸡代码.
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MAXN 2010
using namespace std;
int n,m,top,ans;
int a[MAXN],b[MAXN];
int l0[MAXN],r0[MAXN],l1[MAXN],r1[MAXN];
int L[MAXN],R[MAXN];
char ch[MAXN];
int pre[MAXN],nxt[MAXN];
struct edge {   int to; edge *next; }e[MAXN],*prev[MAXN];
inline void insert(int u,int v) {   e[++top].to=v;e[top].next=prev[u];prev[u]=&e[top];  }
inline void del(int x)  {   pre[nxt[x]]=pre[x];nxt[pre[x]]=nxt[x];  }
int main()
{
    scanf("%d%d",&n,&m);
    for (int i=1;i<=m;i++)  l0[i]=l1[i]=1,r0[i]=r1[i]=m;
    for (int T=1;T<=n;T++)
    {
        scanf("%s",ch+1);int i,j;
        for (i=j=1;i<=m;i++)
            if (ch[i]=='.') ++a[i],++b[i],l0[i]=max(l0[i],j),l1[i]=max(l1[i],j);
            else    b[i]=a[i]+1,a[i]=0,l1[i]=max(l0[i],j),l0[i]=1,j=i+1;
        for (i=j=m;i;i--)
            if (ch[i]=='.') r0[i]=min(r0[i],j),r1[i]=min(r1[i],j);
            else    r1[i]=min(r0[i],j),r0[i]=m,j=i-1;
        memset(prev,0,sizeof(edge*)*(T+2));top=0;
        for (i=m;i;i--) insert(a[i],i),pre[i]=i-1,nxt[i]=i+1;
        pre[1]=nxt[m]=0;
        for (int x=T;x>=0;x--)  for (edge *I=prev[x];I;I=I->next)   L[I->to]=pre[pre[I->to]],del(I->to);
        memset(prev,0,sizeof(edge*)*(T+2));top=0;
        for (i=m;i;i--) insert(a[i],i),pre[i]=i-1,nxt[i]=i+1;
        pre[1]=nxt[m]=0;
        for (int x=T;x>=0;x--)  for (edge *I=prev[x];I;I=I->next)   R[I->to]=nxt[nxt[I->to]],del(I->to);
        for (i=1;i<=m;i++)
        {
            ans=max(ans,(r1[i]-l1[i]+1)*b[i]);
            if(!a[i])   continue;
            j=l0[i];
            if (j>1&&b[j-1]>=a[i])  j=L[i]?L[i]+1:j-1;
            ans=max(ans,(r0[i]-j+1)*a[i]);
            j=r0[i];
            if (j<m&&b[j+1]>=a[i])  j=R[i]?R[i]-1:j+1;
            ans=max(ans,(j-l0[i]+1)*a[i]);
        }
    }
    cout<<ans<<endl;
}