poj2836 状态压缩dp

自己的做法是枚举i,j作为顶点的矩形，然后再更新状态S，但是这种做法是错误的

正解是先把所有矩形对求出来，然后枚举状态S，每个处理每个状态时再枚举已经求出的矩形对，用旧状态更新新状态

#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[];
int main()
{
    int n;
    while(cin>>n && n)
    {
        int state[],dp[<<],area[];
        int cnt=;
        for(int i=;i<n;i++)
        cin>>p[i].x>>p[i].y;
        for(int i=;i<n;i++)
        {
            for(int j=i+;j<n;j++)
            {
                state[cnt]=(<<i|<<j);
                for(int k=;k<n;k++)
                {
                    if((p[i].x-p[k].x)*(p[k].x-p[j].x)>= && (p[i].y-p[k].y)*(p[k].y-p[j].y)>=)//当一个点在这个矩形中时，把这个矩形的包含点的状态更新
                    state[cnt]|=<<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=;i<(<<n);i++)
        dp[i]=0xfffffff;
        dp[]=;
        for(int i=;i<(<<n);i++)//好好想想为什么后枚举矩形？
        {
            for(int j=;j<cnt;j++)
            {
                dp[i|state[j]]=min(dp[i|state[j]],dp[i]+area[j]);
            }
        }
        cout<<dp[(<<n)-]<<endl;
    }
    return ;
}