原题:https://www.luogu.org/problemnew/show/P1578
题解:给一些障碍物,求最大矩形的面积,不能覆盖障碍物。悬线概念:
有效竖线:除了两个端点外,不覆盖任何障碍点的竖直线段。 悬线:上端点覆盖了一个障碍点或达到整个矩形上端的有效竖线。
每个悬线都与它底部的点一 一对应。每个悬线可以向左右扩展就会得到一个矩形。最大矩形一定在这些矩形中。设:
h[i][j] 为点(i,j)悬线的长度
l[i][j] 为点(i,j)悬线向左最远能达到的位置
r[i][j] 为点(i,j)悬线向右最远能达到的位置
考虑状态转移
如果(i-1,j)是障碍点。(i,j)对应的悬线长度1,左右能移动到的位置是整个矩形的左右边界。
如果不是。悬线长应该是h[i-1][j]+1。左端点应该是 max{l[i-1][j] ,(i-1,j)左边最近的障碍点 }右端点应是min{r[i-1][j],右边最近的障碍物}
细节:
1.当长宽很大时要离散化。
2.没有把上下边界形成的矩形算上,要进行特判。
3.本题递推时要用滚动数组,并且表示位置的数组要用short,不然mle
4.tag数组最外一圈是边界,离散化是坐标要+1
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
using namespace std;
const int N=5005;
int x[N],y[N],sx[N],sy[N];
short pre[N][N],nxt[N][N];
short h[2][N],l[2][N],r[2][N];
bool tag[N][N];
int xcnt,ycnt,n,ans;
inline int rd(){
int x=0;int f=1;char s=getchar();
while(!isdigit(s)){if(s=='-')f=-1; s=getchar();}
while(isdigit(s)) x=(x<<1)+(x<<3)+s-'0',s=getchar();
return x*f;
}
int main(){
// freopen("testdata.in","r",stdin);
sx[++xcnt]=0;sx[++xcnt]=rd();
sy[++ycnt]=0;sy[++ycnt]=rd();
n=rd();
for(int i=1;i<=n;i++){
x[i]=rd();y[i]=rd();
sx[++xcnt]=x[i];
sy[++ycnt]=y[i];
}
sort(sx+1,sx+xcnt+1);xcnt=unique(sx+1,sx+xcnt+1)-(sx+1);
sort(sy+1,sy+ycnt+1);ycnt=unique(sy+1,sy+ycnt+1)-(sy+1);
ans=0;
//特判 和上下边界的矩型
for(int j=2;j<=ycnt;j++) ans=max(ans,(sy[j]-sy[j-1])*sx[xcnt]);
for(int i=1;i<=n;i++){
x[i]=lower_bound(sx+1,sx+xcnt+1,x[i])-sx; // tag[1][]为边界
y[i]=lower_bound(sy+1,sy+ycnt+1,y[i])-sy;
tag[x[i]][y[i]]=1;
}
for(int i=1;i<=xcnt;i++) tag[i][1]=tag[i][ycnt]=1;
for(int j=1;j<=ycnt;j++) tag[1][j]=tag[xcnt][j]=1;
for(int i=1;i<=xcnt;i++){
int last=1;
for(int j=2;j<=ycnt;j++) {
if(!tag[i][j]) continue;
for(int k=last+1;k<j;k++) nxt[i][k]=j;
for(int k=last+1;k<j;k++) pre[i][k]=last;
last=j;
}
}
int t=1;
for(int j=1;j<=ycnt;j++) l[t][j]=1,r[t][j]=ycnt;
t^=1;
for(int i=2;i<=xcnt;i++,t^=1){
for(int j=1;j<=ycnt;j++){
if(tag[i-1][j]){
h[t][j]=1;
l[t][j]=1;
r[t][j]=ycnt;
}else{
h[t][j]=h[t^1][j]+1;
l[t][j]=max(l[t^1][j],pre[i-1][j]);
r[t][j]=min(r[t^1][j],nxt[i-1][j]);
}
ans=max(ans,(sx[i]-sx[i-h[t][j]])*(sy[r[t][j]]-sy[l[t][j]]));
}
}
printf("%d\n",ans);
return 0;
}