愤怒的小鸟（angrybirds）

题目描述

【问题描述】

kiana最近沉迷于一款神奇的游戏无法自拔。

简单来说,这款游戏是在一个平面上进行的。

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

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

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

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

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

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

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

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

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

 

【输入格式】

从文件angrybirds.in 中读入数据。

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

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

如果 m = 0 ,表示kiana输入了一个没有任何作用的指令。

如果 m = 1 ,则这个关卡将会满足:至多用 Ln/3 + 1」只小鸟即可消灭所有小猪。

如果 m = 2 ,则这个关卡将会满足:一定存在一种最优解,其中有一只小鸟消灭 了至少「n/31 只小猪。

保证 1 < n < 18 , 0 < m < 2 , 0 < xi , yi  < 10 ,输入中的实数均保留到小数点后两 位。

上文中 ’符号 Lc」和「c1 分别表示对 c 向上取整和向下取整 ’例如: L2.1」= L2.9」= L3.0」=「3.01 =「3.11 =「3.91 = 3 。

 

【输出格式】

输出到文件angrybirds.out 中。

对每个关卡依次输出一行答案。

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

 

【样例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

 

 【样例1说明】

这组数据中一共有两个关卡。

第一个关卡与【问题描述】 中的情形相同’ 2 只小猪分别位于 (1.00，3.00) 和 (3.00，3.00) ’只需发射一只飞行轨迹为 y = -x2  + 4x 的小鸟即可消灭它们。

第二个关卡中有 5 只小猪’但经过观察我们可以发现它们的坐标都在抛物线 y = -x2  + 6x 上’故kiana只需要发射一只小鸟即可消灭所有小猪。

【样例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

 

 

【子任务】

数据的一些特殊规定如下表：

 

测试点编号	n	m	T
1	< 2	= 0	< 10
2			< 30
3	< 3		< 10
4			< 30
5	< 4		< 10
6			< 30
7	< 5		< 10
8	< 6
9	< 7
10	< 8
11	< 9		< 30
12	< 10
13	< 12	= 1
14		= 2
15	< 15	= 0	< 15
16		= 1
17		= 2
18	< 18	= 0	< 5

代码

#include <bits/stdc++.h>
using namespace std;
#define FOR(i,a,b) for(int i=(a),i_##END_=(b);i<=i_##END_;++i)
#define REP(i,a,b) for(int i=(a),i_##BEGIN_=(b);i>=i_##BEGIN_;--i)
#define M 20
#define e 1e-6
typedef double db;
int T,n,m;
struct node{
    db x,y;
}P[M];
int tol,Sit[M*M];
int dp[1<<M];
void Init(){
    FOR(i,1,n)FOR(j,i+1,n){
        if(P[i].x==P[j].x)continue;
        db tmp1=P[i].x*P[i].x,res1=P[i].x;
        db tmp2=P[j].x*P[j].x,res2=P[j].x;
        db aa=res2/res1;
        db a=(P[i].y*aa-P[j].y)/(tmp1*aa-tmp2);
        db b=(P[i].y-a*tmp1)/res1;
        tol++;
        if(a>=0)continue;
        Sit[tol]|=1<<(i-1);
        Sit[tol]|=1<<(j-1);
        FOR(k,1,n){
            if(k==i||k==j)continue;
            if(abs(P[k].x*P[k].x*a+P[k].x*b-P[k].y)<=e)Sit[tol]|=1<<(k-1);
        }
    }
}
int cmp(node a,node b){
    return a.x<b.y;
}
int main(){
    cin>>T;
    while(T--){
        tol=0;
        scanf("%d%d",&n,&m);
        FOR(i,1,n)scanf("%lf%lf",&P[i].x,&P[i].y);
        memset(Sit,0,sizeof(Sit));
        FOR(i,1,n)Sit[++tol]=1<<(i-1);
        Init();
        FOR(i,0,(1<<n)-1)dp[i]=1e9;
        dp[(1<<n)-1]=0;
        REP(i,(1<<n)-1,0){
            FOR(j,1,tol){
                dp[i]=min(dp[i],dp[i|Sit[j]]+1);
            }
        }
        printf("%d\n",dp[0]);
    }
    return 0;
}