P1578 奶牛浴场 悬线法+dp

原题: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;
}