P1578 奶牛浴场
由于n,m很大但是障碍点的数量很小,所以可以枚举障碍点从而找到最大子矩阵。
复杂度为障碍点的数量的平方,一般用于处理n,m很大时的最大子矩阵问题,枚举思想很简单,排序,然后以第一个点为左边界,依次枚举右边界,然后再以第二个点,第三个点为右边界进行枚举。
由于边界问题,细节比较多,实现起来很麻烦,但是也可以当作模板。
#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
struct skt{
int x,y;
}a[5100];
bool cmp1(skt a,skt b){
return a.x<b.x;
}
bool cmp2(skt a,skt b){
return a.y<b.y;
}
int main(void) {
int n,L,W,i,j;
scanf("%d%d%d",&L,&W,&n);
for(i=0;i<n;i++)
scanf("%d%d",&a[i].x,&a[i].y);
a[n].x=0;a[n].y=0;
a[n+1].x=0;a[n+1].y=W;
a[n+2].x=L;a[n+2].y=0;
a[n+3].x=L;a[n+3].y=W;
n+=4;
sort(a,a+n,cmp1);
int ans=0,up,down,len;
for(i=0;i<n;i++){
up=0;down=W;len=L-a[i].x;
for(j=i+1;j<n;j++){
if(len*(down-up)<=ans)
break;
ans=max(ans,(down-up)*(a[j].x-a[i].x));
if(a[j].y>=a[i].y)
down=min(down,a[j].y);
else
up=max(down,a[j].y);
}
}
sort(a,a+n,cmp2);
int l,r;
for(i=0;i<n;i++){
l=0;r=L;len=W-a[i].y;
for(j=i+1;j<n;j++){
if(len*(r-l)<=ans)
break;
ans=max(ans,(r-l)*(a[j].y-a[i].y));
if(a[j].x>=a[i].x)
r=min(r,a[j].x);
else
l=max(l,a[j].x);
}
}
cout<<ans<<"\n";
return 0;
}