感觉这个人写得很不错
#include <stdio.h>
#include <string.h>
#define INF 1000000000
#define MAXN 20
#define MAXS 1<<18
int n,m,c;
int x[MAXN],y[MAXN];
int s[MAXN][MAXN];
int f[MAXS];
int dp(int ss)
{
int count,i,j,ans;
if(f[ss]!=-1)
return f[ss];
for(count=0,i=0; i<n; i++)
if( (1<<i) & ss ) count++;
//位运算,其实是找出在n位二进制数中还有多少个1,即还有多少棵树没有砍
if(count<=c) return f[ss]=0;
//如果剩下的已经小于等于c,那么就是不用再砍的,当前这个状态需要的抢数为0
if(count==1) return f[ss]=1;
//剩下一棵树,那么就用一把枪去消灭
//上述情况都不符合那说明还需要砍树,那么就进行DP构建找出最优方案
f[ss]=INF;
for(i=0; i<n; i++) if( (1<<i)&ss ) //第i棵树还存在
for(j=i+1; j<n; j++) if( (1<<j)&ss ) //第j棵树还存在
{
ans=dp( ss&(~s[i][j]) )+1; //~按位取反运算
//ss&(~s[i][j]) 这个运算就能表示砍掉这条直线上的数并得到下一个状态
if(ans<f[ss]) f[ss]=ans;
}
return f[ss];
}
void solve(int CASE)
{
int i,j,k,ans;
//进行状态压缩,记录在同一条直线上的点
memset(s,0,sizeof(s));
for(i=0; i<n; i++)
for(j=0; j<n; j++)
if(i!=j)
{
for(k=n-1; k>=0; k--)
{
s[i][j]<<=1; //左移一位,为当前的点k腾出位置
if( (y[j] - y[i]) * (x[k] - x[i]) == (y[k] - y[i])*(x[j] - x[i]) )
s[i][j]++; //其实相当于在长度为n的二进制的数对应的位置上标记1
//其实是通过判断斜率来确定,原型为 (y[j]-y[i])/(x[j]-x[i]) == (y[k]-y[i])/(x[k]-x[i])
//若相等可以容易证明i,j,k三点一线
//写成乘法可以消掉特殊情况例如斜率为0或者斜率不存在等等
}
}
c=n-m; //当剩下的树等于小于c的时候说明已经消灭完成
memset(f,-1,sizeof(f));
ans=dp( (1<<n)-1 );
//一开始有2^n-1棵树,剩下c棵或以下的树时构建结束
printf("Case #%d:\n",CASE);
printf("%d\n",ans);
}
int main()
{
int i,T,t;
scanf("%d",&T);
for(t=1; t<=T; t++)
{
scanf("%d%d",&n,&m);
for(i=0; i<n; i++)
scanf("%d%d",&x[i],&y[i]);
solve(t);
if(t!=T) printf("\n");
}
return 0;
}