【YBT2023寒假Day7 B】打怪兽（cdq分治）（斜率优化）

最新推荐文章于 2023-06-01 20:46:22 发布

原创最新推荐文章于 2023-06-01 20:46:22 发布 · 387 阅读

0 ·

CC 4.0 BY-SA版权

要转发的话记得附上链接哦~（不过也不会有人转发的啦）

文章标签：

#cdq分治 #斜率优化

分治同时被 2 个专栏收录

41 篇文章

订阅专栏

斜率优化

9 篇文章

订阅专栏

文章讨论了一个与游戏策略相关的算法问题，其中涉及对怪物的攻击力和血量进行处理，通过计算最小血量来战胜所有怪物。关键策略包括将血量转换为击败次数，解决经典排序问题，以及利用CDQ分治和斜率优化来降低复杂度。代码部分展示了实现这一算法的C++程序。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

打怪兽

题目链接：YBT2023寒假Day7 B

题目大意

有 n 个怪，每个怪有攻击力和血量。
你每次可以选一个怪打 b 的伤害，如果一个怪的血量小于等于 0 就死了。
然后每次你打完之后所有活着的怪会各打你一次。
一开始你可以选两个怪秒杀，问你打死所有怪需要的最小血量。

思路

首先把血量 $c$ 改成要打死的次数 $(c + b - 1) / b$ （就是 $\left\lceil\dfrac{c}{b}\right\rceil$ ）

首先不考虑秒杀，其实是一个经典的排序问题。
（你直接看谁先就直接假设只有这两个怪，谁先打掉的血少谁就放前面）
（至于原因把式子列一下就看得出来）

接下来考虑看秒杀，考虑秒掉之后减少的费用。
我们假设 $i, j (i < j)$ 为秒掉的两个，那减少的费用是：
$a_i(\sum\limits_{k=1}^ic_k-1)+c_i\sum\limits_{k=i+1}^na_k+a_j(\sum\limits_{k=1}^jc_k-c_i-1)+c_j\sum\limits_{k=j+1}^na_k$
（第一个和第三个是它打你的次数减少了，第二个和第四个是你打它的时间省了使得后面的怪打你的次数少了）

会发现可以写成 $A_i+B_j-a_jc_i$ 的形式（其实 $A_i,B_i$ 是同一个数组，都是 $a_i(\sum\limits_{k=1}^ic_i-1)+c_i\sum\limits_{k=i+1}^na_i$ ，前缀后缀预处理一下就行）

那如果我们枚举一个 $j$ ，它就会变成要求 $A_i-a_jc_i$ 的最大值，那这个东西是线段，可以维护上凸壳 $j$ 变大，在凸壳上二分（也可以直接指针找）最大值加上点。
但是一个问题是 $a, c$ 不能同时递增。

于是这种两个变量的考虑通过 cdq 分治减少一个。
大概是给前面一半按 $c_i$ 排序，后面一半按 $a_i$ 排序。
然后推 $A_{i1}-a_jc_{i1}>A_{i2}-a_jc_{i2}$ 这个式子得到比较方式。（就是斜率优化）

代码

#include<cstdio>
#include<iostream>
#include<algorithm> 
#define ll long long

using namespace std;

const int N = 3e5 + 100;
struct node {
	ll a, d, A;
}a[N];
int n, b, tot;
ll ans, sta[N];

bool cmp(node x, node y) {
	return x.a * (x.d - 1) + y.a * (x.d + y.d - 1) < y.a * (y.d - 1) + x.a * (y.d + x.d - 1);
}

bool cmpd(node x, node y) {return x.d < y.d;}
bool cmpa(node x, node y) {return x.a < y.a;}

void cdq(int l, int r) {
	if (l == r) return ;
	int mid = (l + r) >> 1;
	cdq(l, mid); cdq(mid + 1, r);
	sort(a + l, a + mid + 1, cmpd); sort(a + mid + 1, a + r + 1, cmpa);
	tot = 1; int L = 1; sta[0] = sta[1] = 0;
	for (int i = l; i <= mid; i++) {
		while (L < tot && (a[i].A - a[sta[tot]].A) * (a[sta[tot]].d - a[sta[tot - 1]].d) >= (a[sta[tot]].A - a[sta[tot - 1]].A) * (a[i].d - a[sta[tot]].d)) tot--;
		sta[++tot] = i;
	}
	for (int i = mid + 1; i <= r; i++) {
		while (L < tot && (a[sta[L + 1]].A - a[sta[L]].A) >= a[i].a * (a[sta[L + 1]].d - a[sta[L]].d)) L++;
		ans = max(ans, a[i].A + a[sta[L]].A - a[i].a * a[sta[L]].d);
	}
}

int main() {
	freopen("fittest.in", "r", stdin);
	freopen("fittest.out", "w", stdout);
	
	scanf("%d %d", &n, &b);
	for (int i = 1; i <= n; i++) {
		scanf("%lld %lld", &a[i].a, &a[i].d); a[i].d = (a[i].d + b - 1) / b;
	}
	sort(a + 1, a + n + 1, cmp);
	
	ll all = 0, tmp = 0;
	for (int i = 1; i <= n; i++) {
		tmp += a[i].d;
		all += a[i].a * (tmp - 1);
		a[i].A = a[i].a * (tmp - 1);
	}
	tmp = 0;
	for (int i = n; i >= 1; i--) {
		a[i].A += a[i].d * tmp;
		tmp += a[i].a;
	}
	
	cdq(1, n);
	printf("%lld", all - ans);
	
	return 0;
}