洛谷2831(NOIP2016)[愤怒的小鸟]--状压DP

【链接】
洛谷2831

【题目大意】

求经过所有的坐标的二次函数的个数最小值。

【解题报告】

状压DP。

定义$g_{i,j}$表示经过第$i$个点和第$j$个点的二次函数所能经过的节点状态。

定义$f_s$表示目前已经经过的节点状态。

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=20,maxm=(1<<18)+5;
int T,n,m,f[maxm],g[maxn][maxn];
struct wjd
{
    double x,y;
    bool operator < (const wjd &a) const{
        return x<a.x;
    }
}a[maxn];
int Fcmp(double x,double y) {if (fabs(x-y)<1e-10) return 0; if (x<y) return -1; return 1;}
void Work()
{
    scanf("%d%d",&n,&m); m=(1<<n)-1;
    for (int i=0; i<n; i++) scanf("%lf%lf",&a[i].x,&a[i].y);
    sort(a,a+n);
    memset(g,0,sizeof(g));
    for (int i=0; i<n-1; i++)
     for (int j=i+1; j<n; j++)
      if (Fcmp(a[i].x,a[j].x))
       {
        double A,B; A=(a[i].y-a[i].x*a[j].y/a[j].x)/(a[i].x*a[i].x-a[i].x*a[j].x);
        if (Fcmp(A,0)>=0) continue; B=(a[i].y-A*a[i].x*a[i].x)/a[i].x;
        for (int k=0; k<n; k++) if (!Fcmp(A*a[k].x*a[k].x+B*a[k].x,a[k].y)) g[i][j]|=1<<k;
       }
    memset(f,63,sizeof(f)); f[0]=0;
    for (int s=0; s<m; s++)
     for (int i=0; i<n; i++)
      if (!(s&(1<<i)))
       {
        f[s|(1<<i)]=min(f[s|(1<<i)],f[s]+1);
        for (int j=i+1; j<n; j++)
         f[s|g[i][j]]=min(f[s|g[i][j]],f[s]+1);
       }
    printf("%d\n",f[m]);
}
int main()
{
    freopen("2831.in","r",stdin);
    freopen("2831.out","w",stdout);
    scanf("%d",&T);
    while (T--) Work();
    return 0;
}