状压DP死活想不到啊,开始还以为是数据结构题,真是膜拜宋威大大。
包括点的矩形最小一定以这些点为顶点,所以枚举所有点对,构成所有状态。
由于不存在面积为0的矩形,所以当横纵坐标相同时得特判。
然后就是状压dp了。
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<vector>
#include<math.h>
#include<map>
#include<queue>
using namespace std;
struct node
{
int x,y;
}p[25];
int main()
{
int n;
while(cin>>n && n)
{
int state[120],dp[1<<15],area[120];
int cnt=0;
for(int i=0;i<n;i++)
cin>>p[i].x>>p[i].y;
for(int i=0;i<n;i++)
{
for(int j=i+1;j<n;j++)
{
state[cnt]=(1<<i|1<<j);
for(int k=0;k<n;k++)
{
if((p[i].x-p[k].x)*(p[k].x-p[j].x)>=0 && (p[i].y-p[k].y)*(p[k].y-p[j].y)>=0)//当一个点在这个矩形中时,把这个矩形的包含点的状态更新
state[cnt]|=1<<k;
}
if(p[i].x==p[j].x)//特判
area[cnt]=abs(p[i].y-p[j].y);
else if(p[i].y==p[j].y)
area[cnt]=abs(p[i].x-p[j].x);
else
{
area[cnt]=abs(p[i].y-p[j].y)*abs(p[i].x-p[j].x);
}
cnt++;
}
}
for(int i=0;i<(1<<n);i++)
dp[i]=0xfffffff;
dp[0]=0;
for(int i=0;i<(1<<n);i++)//好好想想为什么后枚举矩形?
{
for(int j=0;j<cnt;j++)
{
dp[i|state[j]]=min(dp[i|state[j]],dp[i]+area[j]);
}
}
cout<<dp[(1<<n)-1]<<endl;
}
return 0;
}