洛谷P2831 [Noip2016]愤怒的小鸟

题目描述

• $Kiana$最近沉迷于一款神奇的游戏无法自拔。
• 简单来说，这款游戏是在一个平面上进行的。
• 有一架弹弓位于$(0,0)$处，每次$Kiana$可以用它向第一象限发射一只红色的小鸟，小鸟们的飞行轨迹均为形如$y=ax^2+bx$的曲线，其中$a,b$是$Kiana$指定的参数，且必须满足$a<0$。
• 当小鸟落回地面(即$x$轴）时，它就会瞬间消失。
• 在游戏的某个关卡里，平面的第一象限中有$n$只绿色的小猪，其中第$i$只小猪所在的坐标为$(x_i,y_i)$。
• 如果某只小鸟的飞行轨迹经过了$(x_i,y_i)$,那么第$i$只小猪就会被消灭掉，同时小鸟将会沿着原先的轨迹继续飞行；
• 如果一只小鸟的飞行轨迹没有经过$(x_i,y_i)$，那么这只小鸟飞行的全过程就不会对第$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)$，表示第$i$只小猪坐标为$(x_i,y_i)$。数据保证同一个关卡中不存在两只坐标完全相同的小猪。
• 如果$m=0$，表示$Kiana$输入了一个没有任何作用的指令。
• 如果$m=1$，则这个关卡将会满足：至多用$\lceil\frac{n}{3} + 1\rceil$只小鸟即可消灭所有小猪。
• 如果$m=2$，则这个关卡将会满足：一定存在一种最优解，其中有一只小鸟消灭了至少$\lfloor \frac{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$分别表示对$c$向上取整和向下取整

输出格式

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

输入输出样例

输入样例#1

• 2
  2 0
  1.00 3.00
  3.00 3.00
  5 2
  1.00 5.00
  2.00 8.00
  3.00 9.00
  4.00 8.00
  5.00 5.00

输出样例#1

• 1
  1

输入样例#2

• 3
  2 0
  1.41 2.00
  1.73 3.00
  3 0
  1.11 1.41
  2.34 1.79
  2.98 1.49
  5 0
  2.72 2.72
  2.72 3.14
  3.14 2.72
  3.14 3.14
  5.00 5.00

输出样例#2

• 2
  2
  3

输入样例#3

• 1
  10 0
  7.16 6.28
  2.02 0.38
  8.33 7.78
  7.68 2.09
  7.46 7.86
  5.77 7.44
  8.24 6.72
  4.42 5.11
  5.42 7.79
  8.15 4.99

输出样例#3

• 6

说明

样例解释1

• 这组数据中一共有两个关卡。
• 第一个关卡与问题描述中的情形相同，$2$只小猪分别位于$(1.00,3.00)$和 $(3.00,3.00)$，只需发射一只飞行轨迹为$y = -x^2 + 4x$的小鸟即可消灭它们。
• 第二个关卡中有$5$只小猪，但经过观察我们可以发现它们的坐标都在抛物线 $y = -x^2 + 6x$上，故$Kiana$只需要发射一只小鸟即可消灭所有小猪。

数据范围

原题地址

洛谷P2831 UOJ#265

分析 状态压缩 + 记忆化搜索 + 剪枝

• 很容易想到暴搜枚举两个未被消灭的点来确定一条抛物线，并判断这条抛物线能消灭多少个点，消灭完所有点即为一种方案（当然，也可以每次直接选择消灭某一个点）
• 这样裸搜肯定会超时，我们考虑如何优化：
  
  1. 当我们枚举确定抛物线的第一个点的时候，记录当前枚举到的点$i$，则下一次只要从点$i + 1$开始枚举即可
  2. 最多只有$18$个点，可以将其压缩成一个$18$位的二进制数，若第$i$位为$1$表示第$i$个点已被消灭，那么枚举点的时候直接用当前二进制状态判断是否可选即可
  3. 我们记$f[s]$表示二进制状态为$s$的当前最优解，如果发现目前已发射的小鸟数$G \ge f[s]$，则继续搜下去肯定时不优的，否则就令$f[s] = G$
  4. 当$m = 1$时，我们可以在初始时直接将$f$数组全部赋为$\lceil \frac{n}3 + 1\rceil$（并不懂得m = 2的命令怎么用……）
  5. 记数组$mth[i][j]$表示选取两点$i, j$确定抛物线能消灭多少个点，同样用二进制数存储，考虑到$mth[i][j]$中可能包含一些之前已经被消灭的点，合并状态时应为$(s$&$mth[i][j]) + (s | mth[i][j])$
  6. 改变搜索顺序：先搜索两点确定抛物线的做法，然后再搜索单独消灭一点。一般情况下这样会先计算出较优的解，以避免不优解的重复搜索
• 最后的答案即为$f[2^n - 1]$

注意

• $UOJ$有特殊点专门卡精度，以下为在$UOJ$提交的$AC$代码，注意$long\_double$的使用和$eps$的取值

代码

#include <cstdio>
#include <cmath>
#include <cstring>

using namespace std;
typedef long double ld;
const ld eps = 1e-10;
const int Maxn = 0x3f3f3f3f;
const int N = 20;
ld x[N], y[N]; int n, T, m;
int f[(1 << 20) + 5], mth[N][N];

inline void CkMin(int &a, const int &b) {if (a > b) a = b;}

inline void Dfs(const int &t, const int &G, const int &s)
{
    if (G >= f[s]) return ; f[s] = G;
    if (s == ((1 << n) - 1)) return ;
    for (int i = t; i < n; ++i)
    if (!(s & (1 << i - 1)))
    { 
        for (int j = i + 1; j <= n; ++j) 
        if (!(s & (1 << j - 1)))
        {
            if (!mth[i][j]) continue;
            Dfs(i + 1, G + 1, (s & mth[i][j]) + (s | mth[i][j]));
        }       
    } 
    for (int i = t; i <= n; ++i)
     if (!(s & (1 << i - 1)))
      Dfs(i + 1, G + 1, s | (1 << i - 1));
}

inline int get()
{
    char ch; int res = 0; 
    while ((ch = getchar()) < '0' || ch > '9');
    res = ch - '0';
    while ((ch = getchar()) >= '0' && ch <= '9')
     res = (res << 3) + (res << 1) + ch - '0';
    return res; 
}

inline void put(int x)
{
    if (x > 9) put(x / 10);
    putchar(x % 10 + 48);
}

int main()
{
    T = get();
    while (T--)
    {
        n = get(); m = get();
        if (m & 1) for (int i = 0; i < (1 << n); ++i) f[i] = ceil((ld)n / 3.0 + 1);
         else memset(f, Maxn, sizeof(f));
        for (int i = 1; i <= n; ++i) scanf("%Lf %Lf", &x[i], &y[i]);
        for (int i = 1; i < n; ++i)
         for (int j = i + 1; j <= n; ++j)
         {  
            mth[i][j] = 0;
            ld a = (y[i] / x[i] - y[j] / x[j]) / (x[i] - x[j]),
               b = y[i] / x[i] - a * x[i];
            if (a >= -eps) continue;
            for (int k = 1; k <= n; ++k)
             if (fabs(a * x[k] * x[k] + b * x[k] - y[k]) < eps) 
              mth[i][j] |= (1 << k - 1);
         }
        Dfs(1, 0, 0); put(f[(1 << n) - 1]), putchar('\n');
    }
}