P1034 矩形覆盖

第一反应就是大法师：枚举横着切和竖着切。

感觉这样就已经对了，但后来想想貌似存在这样一种状况。



其实想想这种情况也不会是最优的。但是这题k等于四，总觉得不讨论一下就很不舒服。

然后就A掉了。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
#define rep(i,j,k) for(i=j;i<=k;++i)
#define per(i,j,k) for(i=j;i>=k;--i)
#define ll long long
#define db double
#define mkp(x,y) make_pair(x,y)
#define pii pair<int,int>
#define X first
#define Y second
const int N=51,M=500;
int n,k,ans;pii v
;
bool cmp(pii x,pii y){
return x.Y==y.Y?x.X<y.X:x.Y<y.Y;
}
int cal(int l,int r,int d,int u){
return (r-l)*(u-d);
}
void f1(pii v1,pii v2){
int h,u1=v1.Y,l1=v1.X,d2=v2.Y,r2=v2.X,u3=0,d3=M,l3=M,r3=0,u4=0,d4=M,l4=M,r4=0;
rep(h,1,n)if(v[h]!=v1&&v[h]!=v2){
if(v[h].X<=v1.X&&v[h].Y>=v1.Y)
u1=max(u1,v[h].Y),l1=min(l1,v[h].X);
if(v[h].X>=v2.X&&v[h].Y<=v2.Y)
d2=min(d2,v[h].Y),r2=max(r2,v[h].X);
if(v[h].X>v1.X&&v[h].Y>v2.Y)
u3=max(u3,v[h].Y),d3=min(d3,v[h].Y),l3=min(l3,v[h].X),r3=max(r3,v[h].X);
if(v[h].X<v2.X&&v[h].Y<v1.Y)
u4=max(u4,v[h].Y),d4=min(d4,v[h].Y),l4=min(l4,v[h].X),r4=max(r4,v[h].X);
}
ans=min(ans,cal(l1,v1.X,v1.Y,u1)+cal(v2.X,r2,d2,v2.Y)+cal(l3,r3,d3,u3)+cal(l4,r4,d4,u4));
}
void f2(pii v1,pii v2){
int h,u1=v1.Y,r1=v1.X,d2=v2.Y,l2=v2.X,u3=0,d3=M,l3=M,r3=0,u4=0,d4=M,l4=M,r4=0;
rep(h,1,n)if(v[h]!=v1&&v[h]!=v2){
if(v[h].X>=v1.X&&v[h].Y>=v1.Y)
u1=max(u1,v[h].Y),r1=max(r1,v[h].X);
if(v[h].X<=v2.X&&v[h].Y<=v2.Y)
d2=min(d2,v[h].Y),l2=max(l2,v[h].X);
if(v[h].X<v1.X&&v[h].Y>v2.Y)
u3=max(u3,v[h].Y),d3=min(d3,v[h].Y),l3=min(l3,v[h].X),r3=max(r3,v[h].X);
if(v[h].X>v2.X&&v[h].Y<v1.Y)
u4=max(u4,v[h].Y),d4=min(d4,v[h].Y),l4=min(l4,v[h].X),r4=max(r4,v[h].X);
}
ans=min(ans,cal(v1.X,r1,v1.Y,u1)+cal(l2,v2.X,d2,v2.Y)+cal(l3,r3,d3,u3)+cal(l4,r4,d4,u4));
}
int F1(int le,int ri,int cur);
int F2(int le,int ri,int cur);
int main(){
int i,j,h,x,y;pii v1,v2;
scanf("%d%d",&n,&k);
rep(i,1,n){
scanf("%d%d",&x,&y);
v[i]=mkp(x,y);
}
sort(v+1,v+n+1);
j=0;rep(i,1,n)if(v[i]!=v[j])v[++j]=v[i];n=j;
if(k>=n){puts("0");return 0;}ans=M*M;
if(k==4){
per(i,n-1,1){
rep(j,i+1,n){
v1=v[i],v2=v[j];
if(v1.X==v2.X)f1(v1,v2),f2(v1,v2);
else if(v1.Y==v2.Y)f1(v1,v2),f2(v2,v1);
else{
x=min(v1.Y,v2.Y),y=max(v1.Y,v2.Y);
rep(h,1,n)
if(h!=i&&h!=j&&v[h].X>=v1.X&&v[h].X<=v2.X&&v[h].Y>=x&&v[h].Y<=y)
break;
if(h>n)
if(v1.Y<v2.Y)f1(v2,v1),f1(v1,v2);
else f2(v1,v2),f2(v2,v1);
}
}
}
}
ans=min(ans,F1(1,n,k));
printf("%d\n",ans);
return 0;
}
int F1(int le,int ri,int cur){
if(ri-le<cur)return 0;
if(cur==1){
int i,u=0,d=M,l=M,r=0;
rep(i,le,ri)
u=max(v[i].Y,u),d=min(v[i].Y,d),l=min(v[i].X,l),r=max(v[i].X,r);
return cal(l,r,d,u);
}
int i,j,rtn=M*M;
per(i,ri-1,le)
per(j,cur-1,1)
rtn=min(rtn,F1(le,i,j)+F1(i+1,ri,cur-j));
sort(v+le,v+ri+1,cmp);
per(i,ri-1,le)
per(j,cur-1,1)
rtn=min(rtn,F2(le,i,j)+F2(i+1,ri,cur-j));
sort(v+le,v+ri+1);
return rtn;
}
int F2(int le,int ri,int cur){
if(ri-le<cur)return 0;
if(cur==1){
int i,u=0,d=M,l=M,r=0;
rep(i,le,ri)
u=max(v[i].Y,u),d=min(v[i].Y,d),l=min(v[i].X,l),r=max(v[i].X,r);
return cal(l,r,d,u);
}
int i,j,rtn=M*M;
per(i,ri-1,le)
per(j,cur-1,1)
rtn=min(rtn,F2(le,i,j)+F2(i+1,ri,cur-j));
sort(v+le,v+ri+1);
per(i,ri-1,le)
per(j,cur-1,1)
rtn=min(rtn,F2(le,i,j)+F2(i+1,ri,cur-j));
sort(v+le,v+ri+1,cmp);
return rtn;
}