【状态压缩DP】愤怒的小鸟

题目来源

点我进入ACwing官网提交此题

题目描述

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 想知道，对于每一个关卡，至少需要发射多少只小鸟才能消灭所有的小猪。

由于她不会算，所以希望由你告诉她。

注意:本题除 NOIP 原数据外，还包含加强数据。

输入格式
第一行包含一个正整数 T，表示游戏的关卡总数。

下面依次输入这 T 个关卡的信息。

每个关卡第一行包含两个非负整数 n,m，分别表示该关卡中的小猪数量和 Kiana 输入的神秘指令类型。

接下来的 n 行中，第 i 行包含两个正实数 (xi,yi)，表示第 i 只小猪坐标为 (xi,yi)，数据保证同一个关卡中不存在两只坐标完全相同的小猪。

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

如果 m=1，则这个关卡将会满足：至多用 ⌈n/3+1⌉ 只小鸟即可消灭所有小猪。

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

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

上文中，符号 ⌈c⌉ 和 ⌊c⌋ 分别表示对 c 向上取整和向下取整，例如 ：⌈2.1⌉=⌈2.9⌉=⌈3.0⌉=⌊3.0⌋=⌊3.1⌋=⌊3.9⌋=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。

void dfs(已经走过了多少点， 总条数){
	所有点都走过了， 更新答案。

	枚举每一种之前没有包含的曲线
		把这个曲线包含的所有点做个标记，代表已经被包含在曲线中。进行DFS递归求解。
}

• 再根据暴力代码寻找思路，或者优化

在暴力代码中可能同一条直线被用多次，造成时间上的浪费，所以可以根据这个点进行突破，代表有没有走过可以用0/1来表示，那么我们就会想到用二进制数来记录路径。那么就会想到使用状态压缩DP， 那么尝试使用状态压缩DP思想来进行分析。

• DP分析

因为是直线被重复利用了，那么我们就先将所有直线都预处理出来。
一般考虑动态规划是从前一个状态进行划分，如果我要选择这个状态，那么前一个状态是什么呢？
很显然有很多种情况，数据范围比较小，我们可以都枚举出来。然后利用前一个状态来更新我现在的这个状态。

AC代码

#include<bits/stdc++.h>
using namespace std;

#define _for(i, a, b) for (int i = (a); i < (b); ++i)
#define _rep(i, a, b) for (int i = (a); i <= (b); ++i)
#define For(i, a, b) for (int i = (a); i >= (b); --i)
#define debug(a) cout << #a << " = " << a << endl;
#define mod(x) (x) % MOD
#define ENDL "\n"
#define x first
#define y second
typedef long long ll;
typedef pair<double, double> pdd;
typedef vector<int> vi;

const int N = 18 + 5, M = 1 << 18, INF = 0x3f3f3f3f;
const double eps = 1e-8;
int f[M], route[N][N];
pdd alls[N];

int cmp(double x, double y) {
    if (fabs(x - y) <= eps) return 0;
    else if (x < y) return -1;
    return 1;
}

int main()
{
#ifdef LOCAL
    freopen("data.in", "r", stdin);
#endif
    ios::sync_with_stdio(false); // 取消cin与stdin 的同步
    cout.tie(0), cin.tie(0);

    int T;
    cin >> T;
    while (T--) {
        int n, m;
        cin >> n >> m;

        _for(i, 0, n) cin >> alls[i].x >> alls[i].y;

        memset(route, 0, sizeof route);
        _for(i, 0, n) {
            route[i][i] = 1 << i;
            _for(j, 0, n) {
                double x1 = alls[i].x, y1 = alls[i].y;
                double x2 = alls[j].x, y2 = alls[j].y;
                if (!cmp(x1, x2)) continue;
                double a = (y1 / x1 - y2 / x2) / (x1 - x2);
                double b = y1 / x1 - a * x1;
                if (cmp(a, 0) >= 0) continue;

                int path = 0;
                _for(k, 0, n) {
                    double x = alls[k].x, y = alls[k].y;
                    if (!cmp(a * x * x + b * x, y)) path += 1 << k;
                }
                route[i][j] = path;
            }
        }

        memset(f, 0x3f, sizeof f);
        f[0] = 0;
        _for(i, 0, (1 << n) - 1) {
            int u = 0;
            _for(j, 0, n) if (!(i >> j & 1)) {
                u = j;
                break;
            }
            _for(j, 0, n) f[i | route[u][j]] = min(f[i | route[u][j]], f[i] + 1);
        }
        cout << f[(1 << n) - 1] << ENDL;
    }

    return 0;
}