NOIP2016D2T3-愤怒的小鸟

问题描述

Kiana 最近沉迷于一款神奇的游戏无法自拔。
简单来说，这款游戏是在一个平面上进行的。
有一架弹弓位于(0,0) 处，每次 Kiana 可以用它向第一象限发射一只红色的小鸟，小鸟们的飞行轨迹均为形如 $y=ax^2+bx$的曲线，其中a,b 是Kiana 指定的参数，且必须满足 a<0 ，a,b 都是实数。
当小鸟落回地面（即 x 轴）时，它就会瞬间消失。
在游戏的某个关卡里，平面的第一象限中有 n 只绿色的小猪，其中第 ii 只小猪所在的坐标为$\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还输入了一些神秘的指令，使得自己能更轻松地完成这个游戏。这些指令将在【输入格式】中详述。
假设这款游戏一共有 TT 个关卡，现在 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 ，则这个关卡将会满足：一定存在一种最优解，其中有一只小鸟消灭了至少$\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$。

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

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

思路分析

我用的是DFS，具体见注释

AC代码

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
const int MAX_N=20;
const double rc=1e-8;
int n,m,ans;
double zx[MAX_N],zy[MAX_N],dzx[MAX_N],dzy[MAX_N],aa[MAX_N],bb[MAX_N];

bool same_double(double a,double b){
    if(fabs(a-b)<rc)return true;
    return false;
}

//dfs(抛物线数量，第几只猪，单独猪)
void dfs(int p,int cur,int dd){
    if(p+dd>=ans)return;
    if(cur>n){
        ans=min(ans,p+dd);
        return;
    }   
    bool passed=false;
    for(int i=1;i<=p;i++){
        if(same_double(aa[i]*zx[cur]*zx[cur]+bb[i]*zx[cur],zy[cur])){passed=true;break;}
    }
    if(passed==true){
        dfs(p,cur+1,dd);
    }
    else{
        for(int i=1;i<=dd;i++){
            if(same_double(zx[cur],dzx[i])) continue;
            double a=(zy[cur]*dzx[i]-dzy[i]*zx[cur])/(zx[cur]*zx[cur]*dzx[i]-dzx[i]*dzx[i]*zx[cur]);//求参数a 
            double b=(zy[cur]-zx[cur]*zx[cur]*a)/zx[cur];//将a代入求参数b 
            if(a<0){
                aa[p+1]=a;//记录抛物线的参数a 
                bb[p+1]=b;//记录抛物线的参数b
                //将这一只单独的猪从数组中删去，因为它已经和当前猪组成抛物线，不再单独了 
                double q=dzx[i],w=dzy[i];
                for(int j=i;j<dd;j++)
                {
                    dzx[j]=dzx[j+1];
                    dzy[j]=dzy[j+1];
                }
                dfs(p+1,cur+1,dd-1);
                for(int j=dd;j>i;j--) 
                {
                    dzx[j]=dzx[j-1];
                    dzy[j]=dzy[j-1];
                }
                dzx[i]=q;
                dzy[i]=w;
            }
        }
        dzx[dd+1]=zx[cur];
        dzy[dd+1]=zy[cur];
        dfs(p,cur+1,dd+1);
    }
}

int main(){
    int T;
    cin>>T;
    while (T--){
        ans=0x3f3f3f3f;
        cin>>n>>m;
        for(int i=1;i<=n;i++){
        cin>>zx[i]>>zy[i];
        }
        dfs(0,1,0);
        cout<<ans<<endl;
    }
}