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

最新推荐文章于 2025-09-06 14:18:59 发布

原创最新推荐文章于 2025-09-06 14:18:59 发布 · 318 阅读

1 ·

本内容遵循CC 4.0 BY-SA版权协议

动态规划--状压DP 同时被 2 个专栏收录

15 篇文章

订阅专栏

位运算

9 篇文章

订阅专栏

本文探讨了一款游戏的关卡优化策略，通过分析游戏机制，采用状态压缩算法和动态规划来减少过关所需的小鸟数量。文章详细介绍了如何利用抛物线特性计算最优解，以及在特定条件下如何进一步简化问题。

题目描述

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

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

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

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

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

如果一只小鸟的飞行轨迹没有经过 $(xi,yi)\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$ ，则这个关卡将会满足：一定存在一种最优解，其中有一只小鸟消灭了至少 $⌊n/3⌋\lfloor n/3 \rfloor$ 只小猪。

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

上文中，符号 $⌈c⌉\lceil c \rceil$ 和 $⌊c⌋\lfloor c \rfloor$ 分别表示对 cc 向上取整和向下取整，例如： $⌈2.1⌉=⌈2.9⌉=⌈3.0⌉=⌊3.0⌋=⌊3.1⌋=⌊3.9⌋=3\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]+1)dp[\ s|state[i,j]\ ]=min(dp[\ s|state[i,j]\ ],dp[s]+1)$

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

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

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

最后我们要的答案就是 $d p [(1 < < n) - 1]$
而初始化就是对于所有 $i∈[1,n]i\in[1,n]$ ， $d p [1 < < i - 1] = 1$
以及每条抛物线单独消灭得到的状态res， $d p [r e s] = 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;
}