Luogu 2831 [NOIP2016] 愤怒的小鸟

第一眼看成爆搜的状压dp，膜Chester大神犇。

考虑到三个不在同一直线上的点可以确定一条抛物线，而固定点$(0, 0)$和不在同一直线上这两个条件是题目中给定的，所以我们只要枚举两个点然后暴力算抛物线，然后chk一遍观察一下多少点在这一条抛物线上就行了。

想到状压之后状态和方程就显然了。

注意判解出来的抛物线$a \leqslant 0$的情况， eps开小一点。

因为决策集合的大小是$O(n ^ {2})$级别的，所以时间复杂度为$O(Tn^{2}2^{n})$。

Code：

#include <cstdio>
#include <cstring>
using namespace std;
typedef double db;

const int N = ;
const int S = ( << ) + ;
const int inf = 0x3f3f3f3f;
const db eps = 1e-;

int testCase, n, m, q[N * (N + )], f[S];
bool vis[S];

struct Node {
    db x, y;
} a[N];

inline db fabs(db x) {
    return x >  ? x : -x;
}

inline bool chk0(db x) {
    return fabs(x) < eps;
}

inline void chkMin(int &x, int y) {
    if(y < x) x = y;
}

inline int calc(int u, int v) {
    db a1 = a[u].x * a[u].x, b1 = a[u].x, c1 = a[u].y;
    db a2 = a[v].x * a[v].x, b2 = a[v].x, c2 = a[v].y;
    if(chk0((b1 * a2 - b2 * a1)) || chk0((a1 * b2 - b1 * a2))) return ;
    db y = (c1 * a2 - c2 * a1) / (b1 * a2 - b2 * a1);
    db x = (c1 * b2 - c2 * b1) / (a1 * b2 - b1 * a2);
    if(x >  ||chk0(x)) return ;
    int res = ( << (u - )) | ( << (v - ));
    for(int i = ; i <= n; i++) {
        if(i == u || i == v) continue;
        if(chk0(x * a[i].x * a[i].x + y * a[i].x - a[i].y)) res |= ( << (i - ));
    }
    return res;
}

int main() {
    for(scanf("%d", &testCase); testCase--; ) {
        scanf("%d%d", &n, &m);
        for(int i = ; i <= n; i++) scanf("%lf%lf", &a[i].x, &a[i].y);

        m = ;
        memset(vis, , sizeof(vis));
        for(int i = ; i <= n; i++)
            for(int j = ; j < i; j++) {
                int s = calc(i, j);
                if(vis[s] || s == ) continue;
                vis[s] = ;
                q[++m] = s;
            }
        for(int i = ; i <= n; i++) q[++m] = ( << (i - ));

/*        for(int i = 1; i <= m; i++)
            printf("%d ", q[i]);
        printf("\n");   */

        memset(f, 0x3f, sizeof(f));
        f[] = ;
        for(int s = ; s < ( << n); s++)
            for(int i = ; i <= m; i++)
                chkMin(f[s | q[i]], f[s] + );

        printf("%d\n", f[( << n) - ]);
    }
    return ;
}