愤怒的小鸟 状压dp

今天第一次接触了状压dp，思路就是运用二进制运算来压缩运算空间

**

#include <iostream>
#include "cstring"
#include "vector"
#include "math.h"
using namespace std;
//状态压缩DP问题（状压DP）

#define mm(i) memset(i,0,sizeof(i))
#define exp 1e-6 //精度
int mark[20][20];//mark[i][j]记录的是猪的二进制状态（1表示已经打中，0表示未打中），代表了一条二次曲线，表示该曲线在打中i猪和j猪的同时还会打中哪些猪，
// 例如：mark[1][2]=11，换算成二进制是00001011，则代表打中1猪和2猪的同时也打中了4猪
void PreProcess(vector<double> x,vector<double> y,int n){
    for(int i=1;i<=n;i++){
        mark[i][i]|=1<<(i-1);
        for(int j=i+1;j<=n;j++){
            double a=(y[i]/x[i]-y[j]/x[j])/(x[i]-x[j]);
            if(a>=0){
                continue;
            }
            double b=y[i]/x[i]-a*x[i];
            mark[i][j]|=1<<(i-1);//打中了i猪，将i位标1
            mark[i][j]|=1<<(j-1);//打中了j猪，将j位标1
            for(int k=j+1;k<=n;k++){
                if(fabs(y[k]-a*x[k]*x[k]-b*x[k])<exp){//注意浮点类型求相等时的精度问题，如果不设置绝对值与exp比大小，
                    //直接用if(y[k]==a*x[k]*x[k]+b*x[k])的话，oj会报错
                    mark[i][j]|=1<<(k-1);//打中猪和j猪的同时，打中了k猪，将k位标1
                }
            }
            mark[j][i]=mark[i][j];
        }
    }
}

void dfs(vector<int> status, int n){
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            for(int k=0;k<1<<n;k++){
                status[k | mark[i][j]]=min(status[k | mark[i][j]], status[k] + 1);
                //状态k通过曲线mark[i][j]达到了状态k|mark[i][j]，要么就是通过mark[i][j]这条曲线达到的，这时加1；
                //要么就是之前状态k就已经和k|mark[i][j]相同了，这时就没必要mark[i][j]这条曲线了
                //例如k是0111（前三头猪都被打中了），通过曲线mark[i][j]达到的状态是0011（打中前两头猪），
                // 这时k|mark[i][j]就是0111，所以没必要使用mark[i][j]这条曲线
            }
        }
    }
    printf("%d\n",status[(1 << n) - 1]);
}

int main() {
    int T;
    while(scanf("%d",&T)!=EOF){
        for(int i=0;i<T;i++){
            int n,m;
            scanf("%d%d",&n,&m);
            vector<double> x(n+1),y(n+1);
            mm(mark);
            for(int j=1;j<=n;j++){
                scanf("%lf%lf",&x[j],&y[j]);
            }
            vector<int> status(1<<n);
            for(int k=0;k<1<<n;k++){
                status[k]=0x3f;
            }
            status[0]=0;
            PreProcess(x,y,n);
            dfs(status, n);
        }
    }
    return 0;
}