【BZOJ3672】【UOJ7】【NOI2014】购票-CSDN博客

本文链接：https://blog.csdn.net/qq_39972971/article/details/79410975

【题目链接】

BZOJ
UOJ

【思路要点】

令$f_i$表示从$i$出发到根节点的最小花费，$s_i$表示根节点到节点$i$的路径长度之和。
那么，显然有方程：$$f_i=min_{s_i-l_i≤s_j}\{f_j+(s_i-s_j)p_i+q_i\}$$$$f_i=min_{s_i-l_i≤s_j}\{f_j-s_jp_i\}+s_ip_i+q_i$$
其中$j$是$i$的一个祖先。
令$k>j$即$s_k>s_j$，则$f_{i,k}≤f_{i,j}$的充要条件为：$$f_k-s_kp_i≤f_j-s_jp_i$$$$p_i≥\frac{f_k-f_j}{s_k-s_j}$$
如果除去$s_i-l_i≤s_j$的限制，那么我们可以维护一个凸壳，在上面二分来寻找最优决策点。
但加上距离限制之后，我们会对凸壳进行切割，这导致最优解有可能不在凸壳上，而在凸壳上相邻两个决策点之间。
笔者在自己做这道题时没有想到有效的解决方法，只是在暴力扫描了可能需要的相邻两个决策点之间的决策点，但在BZOJ和UOJ上均测试通过，并没有被Extra Test卡掉，这说明针对这个做法的数据极难构造，对于几乎所有数据，它具有$O(NLogN)$的时间复杂度，是一种优秀的考场做法。
正确的做法应该是树分治：每次分治时优先处理分治重心父亲所在的子树，再用该子树的信息更新其它子树。具体来讲，对其它子树的节点的$s_i-l_i$排序，然后我们只考虑其父亲所在子树中的决策点对其它子树的影响，这样，我们就只需要维护一个支持加点的凸包，并在上面二分决策点就可以了。
时间复杂度$O(NLog^2N)$，但笔者太懒没有写，下面的是A了题的暴力。

【代码】

/*BruteForce which got AC*/
#include<bits/stdc++.h>
using namespace std;
const int MAXN = 200005;
const int MAXLOG = 20;
template <typename T> void read(T &x) {
	x = 0; int f = 1;
	char c = getchar();
	for (; !isdigit(c); c = getchar()) if (c == '-') f = -f;
	for (; isdigit(c); c = getchar()) x = x * 10 + c - '0';
	x *= f;
}
template <typename T> void write(T x) {
	if (x < 0) x = -x, putchar('-');
	if (x > 9) write(x / 10);
	putchar(x % 10 + '0');
}
template <typename T> void writeln(T x) {
	write(x);
	puts("");
}
vector <int> t[MAXN];
int n, type, root;
long double rate[MAXN];
int father[MAXN], up[MAXN][MAXLOG];
long long s[MAXN], l[MAXN], f[MAXN], p[MAXN], q[MAXN];
long double calc(int x, int y) {
	return 1.0 * (f[x] - f[y]) / (s[x] - s[y]);
}
void work(int pos, int fa) {
	s[pos] += s[fa];
	int tmp = fa;
	for (int i = MAXLOG - 1; i >= 0; i--)
		if (s[pos] - s[up[tmp][i]] <= l[pos] && p[pos] <= rate[up[tmp][i]]) tmp = up[tmp][i];
	if (s[pos] - s[up[tmp][0]] <= l[pos] && p[pos] <= rate[tmp]) tmp = up[tmp][0];
	f[pos] = f[tmp] + (s[pos] - s[tmp]) * p[pos] + q[pos];
	if (s[pos] - s[up[tmp][0]] > l[pos]) {
		tmp = father[tmp];
		while (tmp && s[pos] - s[tmp] <= l[pos]) {
			f[pos] = min(f[pos], f[tmp] + (s[pos] - s[tmp]) * p[pos] + q[pos]);
			tmp = father[tmp];
		}
	}
	tmp = fa;
	for (int i = MAXLOG - 1; i >= 0; i--)
		if (calc(pos, up[tmp][i]) <= rate[up[tmp][i]]) tmp = up[tmp][i];
	if (calc(pos, tmp) <= rate[tmp]) tmp = up[tmp][0];
	up[pos][0] = tmp;
	for (int i = 1; i < MAXLOG; i++)
		up[pos][i] = up[up[pos][i - 1]][i - 1];;
	if (pos != 1) rate[pos] = calc(pos, tmp);
	for (unsigned i = 0; i < t[pos].size(); i++)
		work(t[pos][i], pos);
}
int main() {
	read(n), read(type);
	for (int i = 2; i <= n; i++) {
		read(father[i]), t[father[i]].push_back(i);
		read(s[i]), read(p[i]), read(q[i]), read(l[i]);
	}
	work(1, 0);
	for (int i = 2; i <= n; i++)
		writeln(f[i]);
	return 0;
}