题目#
[NOIP2016 提高组] 愤怒的小鸟#
Kiana 最近沉迷于一款神奇的游戏无法自拔。
简单来说，这款游戏是在一个平面上进行的。
有一架弹弓位于 $(0,0)$ 处，每次 Kiana 可以用它向第一象限发射一只红色的小鸟，小鸟们的飞行轨迹均为形如 $y=ax^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$，表示第 $i$ 只小猪坐标为 $(x_i,y_i)$。数据保证同一个关卡中不存在两只坐标完全相同的小猪。
如果 $m=0$，表示Kiana输入了一个没有任何作用的指令。
如果 $m=1$，则这个关卡将会满足：至多用 $\lceil n/3 + 1 \rceil$ 只小鸟即可消灭所有小猪。
如果 $m=2$，则这个关卡将会满足：一定存在一种最优解，其中有一只小鸟消灭了至少 $\lfloor n/3 \rfloor$ 只小猪。
保证 $1\leq n \leq 18$，$0\leq m \leq 2$，$0 < x_i,y_i < 10$，输入中的实数均保留到小数点后两位。
上文中，符号 $\lceil c \rceil$ 和 $\lfloor c \rfloor$ 分别表示对 $c$ 向上取整和向下取整，例如：$\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$。
$(n \leq 18)$
题解#
看到 n 小于等于 18，要么搜索，要么状压
使用状压 dp，设f[S]表示二进制状态下要撞到状态为S集合的猪至少要几次抛物线
预处理出任意两点组成的抛物线能打到哪些猪
状态枚举的时候，~ 根据题解所说～，只需要对每一个状态集合，枚举它第一个没打到的猪即可
因为打猪的顺序并不影响结果
问题#
1.double 判断是否为 0——fabs(x) < eps
2.double 判断正负正常和 0 比较
3. 循环特别多的时候注意初始化别放错位置
代码#
#include<bits/stdc++.h>
using namespace std;
int T, n, m, a, b;
const double eps = 1e-6;
const int N = 19;
double x[N], y[N];
int lg[1 << N];
int lowbit[1 << N];
int f[1 << N];
int curve[N][N];
//int f[1 << N];
void getab(double &a, double &b, double x1, double y1, double x2, double y2)
{
a = (y2 * x1 - y1 * x2)/ (x1 * x2 * (x2 - x1));
b = (y2 * x1 * x1 - y1 * x2 * x2)/ (x1 * x2 * (x1 - x2));
}
int main()
{
for(int i = 2;i <= (1 << 18);i++)
{
lg[i] = lg[i >> 1] + 1;
}
for(int i = 0;i < (1 << 18);i++)
{
lowbit[i] = (~i & (-(~i)));
}
scanf("%d", &T);
while(T --)
{
scanf("%d%d", &n, &m);
for(int i = 1;i <= n;i++)
{
scanf("%lf%lf", &x[i], &y[i]);
}
double a, b;
memset(f, 0x3f, sizeof f);
memset(curve, 0, sizeof curve);
for(int i = 1;i <= n;i++)
{
for(int j = 1;j <= n;j++)
{
if(i == j || fabs(x[i] - x[j]) < eps) continue;
getab(a, b, x[i], y[i], x[j], y[j]);
if(a > 0) continue;
for(int k = 1;k <= n;k++)
{
if(fabs(a * x[k] * x[k] + b * x[k] - y[k]) < eps)
{
curve[i][j] |= (1 << (k - 1));
// cout << "qwq " << i <<"and"<< j << "shares" << k<< " "<< curve[i][j]<<endl;
}
}
}
}
for(int i = 1;i <= n;i++)
curve[i][i] = (1 << (i - 1));
f[0] = 0;
for(int s = 0;s < (1 << n) - 1;s++)
{
for(int j = 1;j <= n;j++)
{
f[s | curve[lg[lowbit[s]] + 1][j]] = min(f[s | curve[lg[lowbit[s]] + 1][j]], f[s] + 1);
}
}
printf("%d\n", f[(1 << n) - 1]);
}
return 0;
}