POJ2836

状压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;
}