NOIP2016愤怒的小鸟

愤怒的小鸟

【题目分析】
比赛的时候，觉得这是近年来noip的常规套路，又是暴搜？可是自己估一下复杂度，就感觉药丸。但是没有办法啊，只会敲一个50分的暴力。
后来大神跟我说，这个是状态压缩…没有学过…
后来在网上自学了一下状压，基本是什么意思懂了：dp就是从一个原来的01状态，向现在的状态转移。只不过因为如果要存每只鸟是否被打过，数组要开17维。

【代码】

#define M 20
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
double a,b,x[M],y[M];
int dp[1<<M],g[M][M];
int n,m,cas;
bool chk(double x,double y){
    if(x>y)return x-y<(1e-6);
    return y-x<(1e-6);
}
void MIN(int &x,int y){
    if(x>y)x=y;
}
int main(){
    int i,j,k;
    scanf("%d",&cas);
    while(cas--){
        scanf("%d %d",&n,&m);
        for(i=1;i<=n;i++)scanf("%lf %lf",&x[i],&y[i]);
        memset(g,0,sizeof(g));
        for(i=1;i<=n;i++){
            for(j=i+1;j<=n;j++){
                if(chk(x[i],x[j]))continue;
                a=(y[j]*x[i]-y[i]*x[j])/((x[j]-x[i])*x[i]*x[j]);
                if(a>=0)continue;
                b=(y[i]-a*x[i]*x[i])/x[i];
                int t=0;
                for(k=1;k<=n;k++){
                    if(chk(a*x[k]+b,y[k]/x[k]))t+=(1<<(k-1));
                }g[i][j]=t;
            }
        }
        memset(dp,63,sizeof(dp));
        dp[0]=0;
        for(i=0;i<(1<<n);i++){
            for(j=1;j<=n;j++){
                if(i&(1<<(j-1)))continue;
                for(k=j+1;k<=n;k++)MIN(dp[i|g[j][k]],dp[i]+1);
                MIN(dp[i|(1<<(j-1))],dp[i]+1);
            }
        }printf("%d\n",dp[(1<<n)-1]);
    }return 0;
}