POJ 2836状压dp

我们可以通过再草稿纸上模拟矩形的枚举方式发现，最优解矩形一定是两个点为对角线的矩形（特殊情况除外），那我们就通过枚举 每两个点来进行状态压缩，把这个两个点所包围的面积中的点压缩再里面，然后通过|表达式，如果最后点集|出来达到了（1<<n），那么就比较最小值。
主要是注意一下这里不能为零矩形的面积。
注意每一次cnt的意义。

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
#include<climits>
#include<stack>
#include<vector>
#include<queue>
#include<set>
#include<map>
#define up(i,a,b)  for(int i=a;i<b;i++)
#define dw(i,a,b)  for(int i=a;i>b;i--)
#define upd(i,a,b) for(int i=a;i<=b;i++)
#define dwd(i,a,b) for(int i=a;i>=b;i--)
//#define local
typedef long long ll;
const double esp = 1e-6;
const double pi = acos(-1.0);
const long long INF = 0x3f3f3f3f;
using namespace std;
typedef pair<int, int> pir;
int n;
struct { int x, y; }coo[101];
int dp[(1<<20)];
int cover[500];
int ss[500];
bool check(int i,int j,int k)//判断k点是否在i，j举行里面
{
	if (coo[k].x >= min(coo[i].x, coo[j].x) && coo[k].x <= max(coo[i].x, coo[j].x) && coo[k].y >= min(coo[i].y, coo[j].y) && coo[k].y <= max(coo[i].y, coo[j].y))
	{
		return true;
	}
	return false;
}
int fab1(int k)//poj不认同fabs
{
	return k < 0 ? -k : k;
}
int main()
{
	while (cin >> n && n)
	{
		up(i, 0, n)
		{
			scanf("%d %d", &coo[i].x, &coo[i].y);
		}
		int cnt = 0;
		up(i, 0, n)
		{
			up(j, 0, n)
			{
				if (j == i)continue;
				ss[cnt] = max(1, fab1(coo[i].x - coo[j].x))*max(1, fab1(coo[i].y - coo[j].y));//计算面积
			//	cout << ss[cnt] << endl;
				cover[cnt] = 0;
				up(k, 0, n)
				{
					if (check(i, j, k))
					{
						cover[cnt] |= (1 << k);//压缩点集的状态
					}
				}
				cnt++;//做完一遍再cnt++
			}
		}
		fill(dp, dp + (1 << 20), 1e9);//全部赋值成最大值，因为求的是最小值
		dp[0] = 0;//没有点的时候面积为零
		up(s, 0, (1 << n))
		{
			if (dp[s] != 1e9);//这里相当于背包问题，如果上一个状态没有被转移到，
			//一定是错误的状态
			up(i, 0, cnt)
			{
				int news = s | cover[i];
				if (news != s)//判断一下是否是无效的状态，相当于已经更新过了
				{
					dp[news] = min(dp[news],ss[i] + dp[s]);
				}
			}
		}
		cout << dp[(1 << n) - 1]<<endl;
	}
	return 0;
}