洛谷 P2831 愤怒的小鸟（状压dp）

题目：https://www.luogu.org/problem/show?pid=2831
思路：
n^3预处理
二进制数s用来表示小猪的状态，1表示被消灭，
f[i][j]表示i点和j点确定的曲线可以干掉的小猪数
dp[s]表示s状态的小猪所需要的最少小鸟数
dp[s | x] = min(dp[s | x],dp[s] + 1)
代码(参考sdfzyhx)：

/*
n^3预处理
二进制数s用来表示小猪的状态，1表示被消灭，
f[i][j]表示i点和j点确定的曲线可以干掉的小猪数
dp[s]表示s状态的小猪所需要的最少小鸟数
dp[s | x] = min(dp[s | x],dp[s] + 1)
*/
#include <bits/stdc++.h>
using namespace std;

const double eps = 1e-10;
const int N = 25;

double x[N],y[N];
int f[N][N],dp[1<<18];
int n,m;

int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        //n^3预处理
        memset(f,0,sizeof(f));
        scanf("%d%d",&n,&m);
        for(int i = 0;i < n;i++)
            scanf("%lf%lf",&x[i],&y[i]);
        for(int i = 0;i < n;i++)
        {
            for(int j = i+1;j < n;j++)
            {
                if(x[i] == x[j])
                    continue;
                double a = (y[j]*x[i] - y[i]*x[j])/((x[j]*x[i]*(x[j]-x[i]));
                double b = (y[j]*x[i]*x[i]-y[i]*x[j]*x[j])/((x[i]-x[j])*x[i]*x[j]);
                if(a >= -eps)//a >= 0
                    continue;
                for(int k = i;k < n;k++)
                    if(fabs(a*x[k]*x[k]+b*x[k]-y[k]) <= eps)//在曲线上
                        f[i][j] |= 1<<k;
            }
        }

        memset(dp,0x3f,sizeof(dp));
        dp[0] = 0;
        for(int s = 0;s < (1<<n)-1;s++)
        {
            int i = 0;
            while(s & (1<<i))
                i++;
            dp[s | (1<<i)] = min(dp[s | (1<<i)],dp[s] + 1);
            for(int j = i;j < n;j++)
                dp[s | f[i][j]] = min(dp[s | f[i][j]],dp[s] + 1);
        }
        printf("%d\n",dp[(1<<n)-1]);
    }
    return 0;
}