洛谷P2831 [NOIP2016] 愤怒的小鸟【状压DP】

题目描述

Kiana 最近沉迷于一款神奇的游戏无法自拔。
简单来说，这款游戏是在一个平面上进行的。

有一架弹弓位于 (0,0)处，每次 Kiana 可以用它向第一象限发射一只红色的小鸟，小鸟们的飞行轨迹均为形如$y=a{x}^2+bx$的曲线，其中a,b是Kiana 指定的参数，且必须满足a < 0，a,b都是实数。

当小鸟落回地面（即$x$ 轴）时，它就会瞬间消失。

在游戏的某个关卡里，平面的第一象限中有n只绿色的小猪，其中第 $i$只小猪所在的坐标为 $\left(x_i,y_i\right)$

如果某只小鸟的飞行轨迹经过了$\left(x_i,y_i\right)$，那么第$i$ 只小猪就会被消灭掉，同时小鸟将会沿着原先的轨迹继续飞行；

如果一只小鸟的飞行轨迹没有经过 $\left(x_i,y_i\right)$，那么这只小鸟飞行的全过程就不会对第 i 只小猪产生任何影响。

例如，若两只小猪分别位于(1,3) 和(3,3) ，Kiana 可以选择发射一只飞行轨迹为 $y=-x^2+4x$的小鸟，这样两只小猪就会被这只小鸟一起消灭。

而这个游戏的目的，就是通过发射小鸟消灭所有的小猪。

这款神奇游戏的每个关卡对 Kiana来说都很难，所以Kiana还输入了一些神秘的指令，使得自己能更轻松地完成这个游戏。这些指令将在【输入格式】中详述。

假设这款游戏一共有$ T$ 个关卡，现在 Kiana想知道，对于每一个关卡，至少需要发射多少只小鸟才能消灭所有的小猪。由于她不会算，所以希望由你告诉她。

输入格式：

第一行包含一个正整数 $T$，表示游戏的关卡总数。

下面依次输入这 $T$ 个关卡的信息。每个关卡第一行包含两个非负整数 $n,m$ ，分别表示该关卡中的小猪数量和 Kiana 输入的神秘指令类型。接下来的 $n $行中，第$i$ 行包含两个正实数 $x_i,y_i$，表示第 ii 只小猪坐标为$(x_i,y_i)$。数据保证同一个关卡中不存在两只坐标完全相同的小猪。

如果 $m=0$ ，表示Kiana输入了一个没有任何作用的指令。
如果$m=1$，则这个关卡将会满足：至多用$ \lceil n/3 + 1 \rceil$只小鸟即可消灭所有小猪。
如果 $m=2$ ，则这个关卡将会满足：一定存在一种最优解，其中有一只小鸟消灭了至少$\lfloor n/3\rfloor$只小猪。

保证$1\le n\le 18$，$0\le m\le 2$，$0<x_i,y_i<10$ ，输入中的实数均保留到小数点后两位。

上文中，符号$\lceil c\rceil$ 和 $\lfloor c\rfloor$ 分别表示对 cc 向上取整和向下取整，例如：$\lceil 2.1\rceil =\lceil 2.9\rceil =\lceil 3.0\rceil =\lfloor 3.0\rfloor =\lfloor 3.1\rfloor =\lfloor 3.9\rfloor =3$

输出格式：

对每个关卡依次输出一行答案。
输出的每一行包含一个正整数，表示相应的关卡中，消灭所有小猪最少需要的小鸟数量。

---

题目分析

8.31update

回头看这篇题解的时候发现了一个小错，虽然不影响正确性，但时间效率就低了很多

//关于判断状态i的第j位是否为1，前一个是之前写的，相当于没用
//把1左移再&运算怎么都不可能出1是吧，当时的我真的很naive
if( (i&(1<<(j-1)))== 1 ) continue;//错误的
if( ((i>>(j-1))&1)== 1 ) continue;//正确的

---

n<=18的范围马上就想到状压

题目中的任意一条抛物线只用两个坐标就可以确定
所以我们把小猪两两组合计算抛物线
若计算得出的a，b符合题目要求(a<0)
则将n个小猪遍历一次，记录这条抛物线还能射中哪些小猪
也就是用一个二进制数表示
第i位为1表示这条抛物线能消灭第i个小猪，为0则不能

用dp[s]表示当前小猪被消灭状态为s
第i位为1表示已消灭第i个小猪，为0则还没有

对于当前状态s
枚举任意两个还没消灭的小猪i，j
设他们的抛物线能消灭的小猪为state[i,j]
那么有状态转移$dp[s|state[i,j]]=min(dp[s|state[i,j]],dp[s]+1)$

但是考虑到有些小猪只能单独被消灭
即他不能和其他小猪算出合法抛物线
或有些小猪被单独消灭可以有更优方案消灭剩余小猪
上述转移并不完全

对于当前状态s
还要枚举还未消灭的小猪并单独消灭(设为i)
那么有状态转移$dp[s|(1<<i-1)]=min(dp[s|(1<<i-1)],dp[s]+1)$
这一步的枚举转移可以与上一步的合并

特别的，对于只能被单独消灭的小猪
可以特判不用继续枚举和他组合的小猪了

最后我们要的答案就是$dp[(1<<n)-1]$
而初始化就是对于所有$i\in [1,n]$，$dp[1<<i-1]=1$
以及每条抛物线单独消灭得到的状态res，$dp[res]=1$

另外这题精度卡的我好烦呐，
非要开long double，
卡常也是十分的刁钻=_=

---

#include<iostream>
#include<cmath>
#include<algorithm>
#include<queue>
#include<cstring>
#include<cstdio>
using namespace std;
typedef double dd;

int read()
{
    int f=1,x=0;
    char ss=getchar();
    while(ss<'0'||ss>'9'){if(ss=='-')f=-1;ss=getchar();}
    while(ss>='0'&&ss<='9'){x=x*10+ss-'0';ss=getchar();}
    return f*x;
}

const int maxn=300010;
int t,n,m;
struct node{dd x,y;}pi[20];
int st[20][20],judge[20];
int dp[maxn];

int calc(int d1,int d2)
{
    dd x1=pi[d1].x,y1=pi[d1].y;
    dd x2=pi[d2].x,y2=pi[d2].y;
    if(x1==x2) return 0;

    dd tx=x1*x1*x2-x2*x2*x1,ty=y1*x2-y2*x1;
    dd a=ty/tx,b=(y1-a*x1*x1)/x1;
    if(a>=0) return 0;//着两只小猪不能被一起消灭

    int res=0;//记录这条抛物线能消灭的小猪
    judge[d1]=judge[d2]=1;
    for(int i=1;i<=n;++i)
    {
        dd x=pi[i].x,y=pi[i].y;
        if(fabs(a*x*x+b*x-y)<1e-7) 
        res|=1<<(i-1),dp[1<<(i-1)]=dp[res]=1;
    }
    return res;
}

int main()
{
    t=read();
    while(t--)
    {
        memset(dp,67,sizeof(dp));
        memset(judge,0,sizeof(judge));
        n=read();m=read();
        for(int i=1;i<=n;++i)
        {
            scanf("%lf%lf",&pi[i].x,&pi[i].y);
            dp[1<<(i-1)]=1;//初始化只射一只小猪
        }

        for(int i=1;i<=n;++i)
        for(int j=1;j<i;++j)
        st[i][j]=calc(i,j);//每两只小猪一起消灭的抛物线

        for(int i=1;i<=(1<<n)-1;++i)//枚举每个状态
        {
            for(int j=1;j<=n;++j)//枚举下一个没被消灭的小猪
            {
                if( ((i>>(j-1))&1)== 1 ) continue;
                if(!judge[j]){dp[i|(1<<(j-1))]=min(dp[i|(1<<(j-1))],dp[i]+1);continue;}
                //只能被单独消灭
                for(int k=1;k<j;++k)//枚举另一个小猪一起组成抛物线
                {
                    if( (i&(1<<(k-1))) == 1 )continue;
                    dp[i|st[j][k]]=min(dp[i|st[j][k]],dp[i]+1);
                }
                dp[i|(1<<(j-1))]=min(dp[i|(1<<(j-1))],dp[i]+1);
                //可能存在单独消灭更优
            }
        }
        printf("%d\n",dp[(1<<n)-1]);
    }
    return 0;
}