uva11008 附

感觉这个人写得很不错

#include <stdio.h>#include <string.h>#define INF 1000000000#define MAXN 20#define MAXS 1<<18int 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;}