「NOIP2016提高组」 愤怒的小鸟

题目描述

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想知道，对于每一个关卡，至少需要发射多少只小鸟才能消灭所有的小猪。由于她不会算，所以希望由你告诉她。

$1\le n\le 18,0\le m\le 2,T\le 30$

题解

状压DP

首先解释一下题目给的$m$其实是没有任何用处的，可以直接忽略

看到$n\le 18$，可以猜到正解不是暴力就是状压

而每个猪要么被消灭要么还没被消灭，只有两种状态，就更加可以肯定这题采用状压

我们先预处理任意两个点确定的抛物线所能经过哪些点，也用二进制表示，记为$fly[i][j]$，注意横坐标相等的情况

设$f[s]$表示目前被消灭的猪的情况是$s$的最少使用的鸟的数量

那么转移方程就可以简单推出：

$f[s|fly[i][j]]=min(f[s|fly[i][j]],f[s]+1)$

特判一下

时间复杂度$O(2^n{n}^2)$，有点接近极限

那么看一下能不能优化

想到如果找到一个还没有被消灭的猪，我们就可以直接退出了

为什么呢？

因为如果当前$i$这个猪没有消灭，而直接去消灭后面的猪，那么$i$最后还是要消灭，因此之前的转移其实都是无用转移，所以我们只用转移后面第一个没有消灭的猪就好了

时间复杂度降到了$O(2^{n}n)$！！！

Code

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#define db double
#define mxdb 1e-6
using namespace std;
int t,n,m,fly[20][20],f[300000];
db a,b,x[20],y[20];
int main()
{
	freopen("angrybirds.in","r",stdin);
	freopen("angrybirds.out","w",stdout);
	scanf("%d",&t);
	while (t--)
	{
		scanf("%d%d",&n,&m);
		for (int i=1;i<=n;++i)
			scanf("%lf%lf",&x[i],&y[i]);
		memset(fly,0,sizeof(fly));
		for (int i=1;i<n;++i)
			for (int j=i+1;j<=n;++j)
			{
				if (fabs(x[i]-x[j])<mxdb) continue;
				a=(y[i]/x[i]-y[j]/x[j])/(x[i]-x[j]);
				if (a>=0) continue;
				b=y[i]/x[i]-x[i]*a;
				for (int k=1;k<=n;++k)
					if (fabs(y[k]/x[k]-(a*x[k]+b))<mxdb) fly[i][j]|=1<<(k-1);//预处理fly[i][j]
			}
		memset(f,0X3f3f3f3f,sizeof(f));
		f[0]=0;//初始化
		for (int s=0;s<(1<<n);++s)
			for (int i=1;i<=n;++i)
				if (!(s&(1<<(i-1))))
				{
					for (int j=1;j<=n;++j)
					{
						if (i==j) f[s|(1<<(i-1))]=min(f[s|(1<<(i-1))],f[s]+1);
						if (fabs(x[i]-x[j])<mxdb) continue;
						f[s|fly[i][j]]=min(f[s|fly[i][j]],f[s]+1);
					}
					break;//找到第一只没有消灭的猪就可以直接返回
				}
		printf("%d\n",f[(1<<n)-1]);
	}
	return 0;
}