UVa 11008 - Antimatter Ray Clearcutting（状态压缩dp）

题目链接
思路：dp[s]表示砍树的状态为s的最小花费，那么每次枚举其中两个点，以这两个点为直线上的点有多少个，那么出去这些点可以得到状态s1，那么dp[s]可以由dp[s1]+1得到。
处于一条直线上的点可以首先与处理出来，然后对于只有一个点的状态特殊处理一下

#include<bits/stdc++.h>
using namespace std;
const int maxn=110;
const int INF=0x3f3f3f3f;
int N,M;
int dp[1<<16];
struct Node
{
    int x,y;
}p[maxn];
int a[maxn];
int g[maxn][maxn];
void getK(int i,int j,pair<int,int> &k)
{
    if(p[i].x==p[j].x)k=make_pair(INF,INF);
    else if(p[i].y==p[j].y)k=make_pair(0,0);
    else
    {
        int g=__gcd(abs(p[i].y-p[j].y),abs(p[i].x-p[j].x));
        k=make_pair((p[i].y-p[j].y)/g,(p[i].x-p[j].x)/g);
    }
}
void solve()
{
    memset(dp,INF,sizeof(dp));
    dp[0]=0;
    for(int S=1;S<(1<<N);S++)
    {
        int cnt=0;
        for(int i=0;i<N;i++)
            if(S&(1<<i))a[cnt++]=i;
        if(cnt==1)
        {
            dp[(1<<(a[0]))]=1;
            continue;
        }
        for(int i=0;i<cnt;i++)
            for(int j=i+1;j<cnt;j++)
            {
                int tmp=(g[a[i]][a[j]]&S);
                tmp=(S^tmp);
                dp[S]=min(dp[S],dp[tmp]+1);
            }
    }
    int ans=INF;
    for(int S=0;S<(1<<N);S++)
    {
        int cnt=0;
        for(int j=0;j<N;j++)
            if(S&(1<<j))cnt++;
        if(cnt>=M)ans=min(dp[S],ans);
    }
    printf("%d\n",ans);
}
int main()
{
    int T,cas=1;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&N,&M);
        for(int i=0;i<N;i++)scanf("%d%d",&p[i].x,&p[i].y);
        memset(g,0,sizeof(g));
        for(int i=0;i<N;i++)
        {
            for(int j=i+1;j<N;j++)
            {
               g[i][j]|=(1<<i);g[i][j]|=(1<<j);
               pair<int,int> K;
               getK(i,j,K);
               for(int k=0;k<N;k++)
               {
                   if(k==i||k==j)continue;
                   pair<int,int> tmp;
                   getK(i,k,tmp);
                   if(tmp==K)g[i][j]|=(1<<k);
               }
            }
        }
        if(cas!=1)printf("\n");
        printf("Case #%d:\n",cas++);
        solve();
    }
    return 0;
}