uva11008(dp + 状态压缩)

题目的的意思就是有ｎ棵树，现在你要砍掉其中的ｍ棵．

你有一把枪，可以砍掉一条直线上的所有树，问砍掉ｍ棵最少几枪．

首先感谢Titanium的博客提供的帮助

http://www.cnblogs.com/scau20110726/archive/2012/09/28/2707866.html

参考了一下，虽然之前也做过状态压缩的题，但是就是没有往这方面想，看了下别人的题解才茅塞顿开．

状态压缩主要是应用于数量不多，而且状态只有两种的情况．

像这题，树最多１６棵，只有已经被砍，和没有被砍两种两种状态，我们可以用一个二进制数来代表现在这些数的状态．

大家可以参考一下上面那篇博文．

ＡＣ代码：

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;

const int N = 20;
const int MAX = 1 << 19;
const int INF = 0x3f3f3f3f;
int n,m,left;
struct point {
	int x;
	int y;
}p[N];
int f[MAX];
int s[N][N];
int dp(int st) {
	int& ans = f[st];
	if(ans != -1)
		return ans;
	int num = 0 ;
	for(int i = 0 ; i < n ; i++) {
		if((1 << i)&st)
			num++;
	}
	if(num <= left)
		return ans = 0;
	if(num == 1) {
		return ans = 1;
	}
	ans = INF;
	for (int i = 0 ; i < n ;i++) {
		for (int j = 0 ; j < n ;j++) {
			if(((1 << i) & st) &&((1 << j) & st) ) {
				ans = min(ans , dp(st & (~s[i][j])) + 1);
			}
		}
	}
	return ans;
}
int main () {
	int t;
	int cas = 1;
	scanf("%d",&t);
	while(t--) {
		memset(f, -1 , sizeof(f));
		scanf("%d%d",&n,&m);
		for(int i = 0 ; i < n ; i++) {
			scanf("%d%d",&p[i].x , &p[i].y);
		}
		for(int i = 0 ;  i < n ; i++) {
			for(int j = 0 ; j < n ; j++) {
				if(i == j)
					continue;
				for(int k = n -1 ; k >= 0 ;k--) {
					s[i][j] <<= 1;
					if((p[j].y - p[i].y)*(p[k].x - p[i].x) == (p[k].y - p[i].y)*(p[j].x - p[i].x)) {
						s[i][j]++;
					}
				}
			}
		}
		left = n - m;
		int res = dp((1 << n) - 1);
		printf("Case #%d:\n%d\n",cas++ , res);
		if(t)
			printf("\n");
	}
}