[NOIp2016提高组]愤怒的小鸟

题目大意：
　　平面直角坐标系的第一象限有n(n<=18)个点，你可以每次给出一个形如y=ax^2+bx的函数把经过这条函数的点消掉，问消掉所有的点至少要多少函数？

思路：
　　枚举每两个点对，可以唯一确定一条函数，再枚举第三个点，判断一下是否会经过这条函数。
　　状态压缩一下记录每条函数能消掉那些点。
　　然后就是一个简单的状压DP。
　　一开始由于把fabs打成了abs，样例都过不去，调了半个下午。
　　洛谷上随便A，UOJ上被Extra数据hack掉了，把eps从1e-7改成1e-10就A了。

#include<cmath>
#include<cstdio>
#include<cctype>
#include<algorithm>
inline int getint() {
    register char ch;
    while(!isdigit(ch=getchar()));
    register int x=ch^'0';
    while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
    return x;
}
const int N=18;
const double eps=1e-10;
struct Point {
    double x,y;
};
Point p[N];
int sit[N*N],s;
int f[1<<N];
inline double pow(const double &x) {
    return x*x;
}
int main() {
    for(register int T=getint();T;T--) {
        const int n=getint();getint();
        for(register int i=0;i<n;i++) {
            scanf("%lf%lf",&p[i].x,&p[i].y);
        }
        for(register int i=0;i<n;i++) {
            sit[i]=1<<i;
        }
        s=n;
        for(register int i=0;i<n;i++) {
            for(register int j=i+1;j<n;j++) {
                const double x1=p[i].x,x2=p[j].x,y1=p[i].y,y2=p[j].y;
                if(pow(x1)*x2==pow(x2)*x1) continue;
                const double a=(x2*y1/x1-y2)/(x1-x2)/x2,b=((pow(x2)*y1/pow(x1)-y2)/(pow(x2)/x1-x2));
                if(a>-eps) continue;
                sit[s]=(1<<i)|(1<<j);
                for(register int k=j+1;k<n;k++) {
                    if(fabs(a*pow(p[k].x)+b*p[k].x-p[k].y)<eps) sit[s]|=1<<k;
                }
                s++;
            }
        }
        for(register int i=1;i<(1<<n);i++) {
            f[i]=1e9;
        }
        for(register int i=0;i<(1<<n);i++) {
            for(register int j=0;j<s;j++) {
                f[i|sit[j]]=std::min(f[i|sit[j]],f[i]+1);
            }
        }
        printf("%d
",f[(1<<n)-1]);
    }
    return 0;
}