【PA2014Final】【BZOJ3726】Wykladzina
2016-03-28 17:11
274 查看
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题解.
细节上并不是很好理解…需要认真地体会一下..
有一块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; }
相关文章推荐
- 2015BJOI day1第二题 树的同构family
- #345 Div2 B. Beautiful Paintings
- CF621C - Wet Shark and Flowers
- CF624D - Array GCD
- 【BZOJ3251】树上三角形
- 【HNOI2014】【BZOJ3573】米特运输
- 【BZOJ4414】数量积
- 便利店选址
- Fuel Stops
- hdu 4472 Count (2012 ACM-ICPC 成都现场赛)
- hdu 4714 Tree2cycle
- HDU 4277 USACO ORZ
- HDU 4061 A Card Game
- HDU 4070 Phage War
- hdu 4696 Answers
- hdu 4782 Beautiful Soupz
- hdu 4790 Just Random
- URAL 1920 Titan Ruins: the Infinite Power of Magic
- ZOJ Candies 长沙
- ZOJ Goldbach 长沙