UVA - 11008 Antimatter Ray Clearcutting

状压dp，值得注意的地方还是很多的。

开始时始终思路不畅，不太清楚究竟该怎么定义状态，所以搁置了很久。后来才明白，由于最后剩下的树的状态有很多种，而开始所有树还在的状态只有一种，所以可以把状态定义为，当前还剩这么些树的时候，要使剩下的树小于n-m棵至少需要的射线数，这样，开始所有树都在的情况就是最大的父状态。这里要注意总结，一般最特殊的状态总是作为最大的父状态。

然后是状态压缩上一些编程技巧和能力，这里不再说了。

注意的是，搜索的时候，可能i和树j只剩下一颗，要处理这种情况，而不能要求i和j都在，不然会wa。

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#define MAX 16
#define INF 0x7f7f7f7f
using namespace std;

int n,m,s[MAX][2],line[MAX][MAX],dp[1<<MAX];

void init(){
	cin>>n>>m;
	for(int i=0;i<n;i++)cin>>s[i][0]>>s[i][1];
	for(int i=0;i<n-1;i++){
		for(int j=i+1;j<n;j++){
			line[i][j]=0;
			line[i][j]+=1<<i,line[i][j]+=1<<j;
			for(int t=j+1;t<n;t++){
				if((s[j][1]-s[i][1])*(s[t][0]-s[i][0])==(s[t][1]-s[i][1])*(s[j][0]-s[i][0]))
					line[i][j]+=1<<t;
			}
		}
	}
	memset(dp,0x7f,sizeof(dp));
}

void dfs(int state){
	//cout<<'a';
	if(dp[state]!=INF)
		return ;
	int left=0;
	for(int i=0;i<n;i++){
		if(state&(1<<i))
			left++;
	}
	if(left<=n-m){
		dp[state]=0;
		return;
	}
	for(int i=0;i<n-1;i++){
		//if(!(state&(1<<i)))
		//	continue;
		for(int j=i+1;j<n;j++){
			if(!(state&(1<<j))&&!(state&(1<<i)))
				continue;
			dfs(state&(~line[i][j]));
			dp[state]=min(dp[state],dp[state&(~line[i][j])]+1);
		}
	}
}

int main(){
	int T,cases=1;
	cin>>T;
	while(T--){
		init();
		if(cases!=1)
			cout<<endl;
		cout<<"Case #"<<cases++<<":"<<endl;
		if(m==0){
			cout<<0<<endl;
			continue;
		}
		if(n<=2){
			cout<<1<<endl;
			continue;
		}
		dfs((1<<n)-1);
		cout<<dp[(1<<n)-1]<<endl;
	}
	return 0;
}