题目传送门

https://lydsy.com/JudgeOnline/problem.php?id=4009

题解

考虑怎样的情况就会有一个链覆盖另一个链。

设被覆盖的链为 \(a - b\),覆盖的链为 \(x - y\)。假设有 \(dfn[a] < dfn[b], dfn[x] < dfn[y]\)

那么如果 \(a\) 是 \(b\) 的祖先,那么令 \(g\) 为 \(a\) 的子树中包含 \(b\) 的点,那么 \(x, y\) 中有一个点在 \(g\) 的子树外面,一个在 \(b\) 子树里面,即 \(x \in [1, dfn[g] - 1], y \in [dfn[b], dfn[b] + siz[b] - 1]\),\(x \in [dfn[b], dfn[b] + siz[b] - 1], y \in [dfn[g] + siz[g], n]\)。

否则,那么 \(x, y\) 必然一个在 \(a\) 子树中,一个在 \(b\) 子树中。即 \(x \in [dfn[a], dfn[a] + siz[a] - 1]\),\(y \in [dfn[b], dfn[b] + siz[b] - 1]\)。

那么我们不妨把上面的问题看成求出包含一个点的矩阵的权值的第 \(k\) 小。

这个问题似乎可以用扫描线+动态主席树来做。(就是带修改区间 \(k\) 小值啊)

但是应该整体二分更好写一写,整体二分以后转化为求出一个点被几个矩形覆盖,也是扫描线+树状数组维护。


但是我似乎已经忘了整体二分怎么写了,然后学了很久,写的时候写出了一堆的 bug,调了一整天才调完。

不过这个题目把询问转化为矩形维护的相关操作是一个很好的思想。


时间复杂度 \(O(q\log^2n)\)。

#include<bits/stdc++.h>

#define fec(i, x, y) (int i = head[x], y = g[i].to; i; i = g[i].ne, y = g[i].to)
#define dbg(...) fprintf(stderr, __VA_ARGS__)
#define File(x) freopen(#x".in", "r", stdin), freopen(#x".out", "w", stdout)
#define fi first
#define se second
#define pb push_back template<typename A, typename B> inline char smax(A &a, const B &b) {return a < b ? a = b, 1 : 0;}
template<typename A, typename B> inline char smin(A &a, const B &b) {return b < a ? a = b, 1 : 0;} typedef long long ll; typedef unsigned long long ull; typedef std::pair<int, int> pii; template<typename I> inline void read(I &x) {
int f = 0, c;
while (!isdigit(c = getchar())) c == '-' ? f = 1 : 0;
x = c & 15;
while (isdigit(c = getchar())) x = (x << 1) + (x << 3) + (c & 15);
f ? x = -x : 0;
} const int N = 40000 + 7; #define lowbit(x) ((x) & -(x)) int n, m, Q, dfc, cnt;
int dep[N], f[N], siz[N], son[N], top[N], dfn[N], pre[N], ans[N];
std::priority_queue<std::pair<int, pii> > q;
struct Matrix {
int x1, x2, y1, y2, v;
inline bool operator < (const Matrix &b) const { return x1 < b.x1; }
} a[N << 1], a1[N << 1], a2[N << 1];
struct Point {
int x, y, k, *ans;
inline bool operator < (const Point &b) const { return x < b.x; }
} b[N], b1[N], b2[N]; struct Edge { int to, ne; } g[N << 1]; int head[N], tot;
inline void addedge(int x, int y) { g[++tot].to = y, g[tot].ne = head[x], head[x] = tot; }
inline void adde(int x, int y) { addedge(x, y), addedge(y, x); } inline void dfs1(int x, int fa = 0) {
dep[x] = dep[fa] + 1, f[x] = fa, siz[x] = 1;
for fec(i, x, y) if (y != fa) dfs1(y, x), siz[x] += siz[y], siz[y] > siz[son[x]] && (son[x] = y);
}
inline void dfs2(int x, int pa) {
top[x] = pa, dfn[x] = ++dfc, pre[dfc] = x;
if (!son[x]) return; dfs2(son[x], pa);
for fec(i, x, y) if (y != f[x] && y != son[x]) dfs2(y, y);
}
inline int gson(int x, int p) {
int g = 0;
while (top[x] != top[p]) g = top[x], x = f[g];
return x == p ? g : son[p];
}
inline bool intr(int x, int p) { return dfn[x] >= dfn[p] && dfn[x] <= dfn[p] + siz[p] - 1; } namespace BIT {
int s[N];
inline void qadd(int x, int k) { for (; x <= n; x += lowbit(x)) s[x] += k; }
inline int qsum(int x) {
int ans = 0;
for (; x; x -= lowbit(x)) ans += s[x];
return ans;
}
inline void qadd(int l, int r, int k) { qadd(l, k), qadd(r + 1, -k); }
}
using BIT::qadd;
using BIT::qsum; inline void solve(int L, int R, int st, int ed, int l, int r) {
if (l > r) return;
if (L == R) {
for (int i = l; i <= r; ++i) *b[i].ans = L;
return;
}
int M = (L + R) >> 1, n1 = 0, n2 = 0, j = l - 1;
for (int i = st; i <= ed; ++i) {
if (a[i].v <= M) {
qadd(a[i].y1, a[i].y2, 1);
q.push(std::make_pair(-a[i].x2, pii(a[i].y1, a[i].y2)));
}
if (i != ed && a[i].x1 == a[i + 1].x1) continue; while (j < r && (i == ed || b[j + 1].x < a[i + 1].x1)) {
++j;
while (!q.empty() && -q.top().fi < b[j].x) qadd(q.top().se.fi, q.top().se.se, -1), q.pop();
int cnt = qsum(b[j].y);
if (cnt >= b[j].k) b1[++n1] = b[j];
else b2[++n2] = b[j], b2[n2].k -= cnt;
}
}
while (!q.empty()) qadd(q.top().se.fi, q.top().se.se, -1), q.pop();
int m1 = 0, m2 = 0;
for (int i = st; i <= ed; ++i) if (a[i].v <= M) a1[++m1] = a[i]; else a2[++m2] = a[i];
assert(m1 + m2 == ed - st + 1), assert(n1 + n2 == r - l + 1);
std::copy(a1 + 1, a1 + m1 + 1, a + st), std::copy(a2 + 1, a2 + m2 + 1, a + st + m1);
std::copy(b1 + 1, b1 + n1 + 1, b + l), std::copy(b2 + 1, b2 + n2 + 1, b + l + n1);
solve(L, M, st, st + m1 - 1, l, l + n1 - 1), solve(M + 1, R, st + m1, ed, l + n1, r);
} inline void work() {
std::sort(a + 1, a + cnt + 1);
std::sort(b + 1, b + Q + 1);
solve(0, 1e9, 1, cnt, 1, Q);
for (int i = 1; i <= Q; ++i) printf("%d\n", ans[i]);
} inline void init() {
read(n), read(m), read(Q);
int x, y;
for (int i = 1; i < n; ++i) read(x), read(y), adde(x, y);
dfs1(1), dfs2(1, 1);
for (int i = 1; i <= m; ++i) {
int x, y, k;
read(x), read(y), read(k);
assert(x != y);
if (dfn[x] > dfn[y]) std::swap(x, y);
if (intr(y, x)) {
x = gson(y, x);
if (dfn[x] > 1) a[++cnt] = (Matrix){ 1, dfn[x] - 1, dfn[y], dfn[y] + siz[y] - 1, k };
if (dfn[x] + siz[x] - 1 < n) a[++cnt] = (Matrix){ dfn[y], dfn[y] + siz[y] - 1, dfn[x] + siz[x], n, k };
} else a[++cnt] = (Matrix){ dfn[x], dfn[x] + siz[x] - 1, dfn[y], dfn[y] + siz[y] - 1, k };
}
for (int i = 1; i <= Q; ++i) {
int x, y, k;
read(x), read(y), read(k);
if (dfn[x] > dfn[y]) std::swap(x, y);
b[i] = (Point){ dfn[x], dfn[y], k, ans + i };
}
} int main() {
#ifdef hzhkk
freopen("hkk.in", "r", stdin);
#endif
init();
work();
fclose(stdin), fclose(stdout);
return 0;
}

bzoj4009 [HNOI2015]接水果 整体二分+扫描线+树状数组+dfs序的更多相关文章

  1. BZOJ 4009: [HNOI2015]接水果 (整体二分+扫描线 树状数组)

    整体二分+扫描线 树状数组 具体做法看这里a CODE #include <cctype> #include <cstdio> #include <cstring> ...

  2. 【BZOJ4009】[HNOI2015]接水果 DFS序+整体二分+扫描线+树状数组

    [BZOJ4009][HNOI2015]接水果 Description 风见幽香非常喜欢玩一个叫做 osu!的游戏,其中她最喜欢玩的模式就是接水果.由于她已经DT FC 了The big black, ...

  3. [bzoj4009] [HNOI2015]接水果 整体二分+扫描线+dfs序+树状数组

    Description 风见幽香非常喜欢玩一个叫做 osu!的游戏,其中她最喜欢玩的模式就是接水果. 由于她已经DT FC 了The big black, 她觉得这个游戏太简单了,于是发明了一个更 加 ...

  4. 【BZOJ】2434: [Noi2011]阿狸的打字机 AC自动机+树状数组+DFS序

    [题意]阿狸喜欢收藏各种稀奇古怪的东西,最近他淘到一台老式的打字机.打字机上只有28个按键,分别印有26个小写英文字母和'B'.'P'两个字母. 经阿狸研究发现,这个打字机是这样工作的: l 输入小写 ...

  5. [BZOJ4009][HNOI2015]接水果(整体二分)

    [HNOI2015]接水果 时间限制:60s      空间限制:512MB 题目描述 风见幽香非常喜欢玩一个叫做 osu!的游戏,其中她最喜欢玩的模式就是接水果. 由于她已经DT FC 了The b ...

  6. BZOJ.4009.[HNOI2015]接水果(整体二分 扫描线)

    LOJ BZOJ 洛谷 又是一个三OJ rank1!=w= \(Description\) (还是感觉,为啥非要出那种题目背景啊=-=直接说不好么) 给定一棵树和一个路径集合(每条路径有一个权值).\ ...

  7. Luogu3527 POI2011 Meteors 整体二分、树状数组、差分

    传送门 比较板子的整体二分题目,时限有点紧注意常数 整体二分的过程中将时间在\([l,mid]\)之间的流星使用树状数组+差分进行维护,然后对所有国家查看一遍并分好类,递归下去,记得消除答案在\([m ...

  8. [bzoj2527][Poi2011]Meteors_整体二分_树状数组

    Meteors bzoj-2527 Poi-2011 题目大意:题目链接. 注释:略. 想法: 首先答案可以离线,且具有单调性. 这里的单调性就是随着时间的推移,每个国家收集的陨石数增加. 不难想到整 ...

  9. [bzoj2738]矩阵乘法_整体二分_树状数组

    矩阵乘法 bzoj-2738 题目大意:给定一个$n*n$的矩阵.每次给定一个矩阵求矩阵$k$小值. 注释:$1\le n\le 500$,$1\le q\le 6\cdot 10^4$. 想法: 新 ...

随机推荐

  1. Splinter 的认识和基础应用

    Splinter 是一个使用Python开发的开源web应用测试程序,它可以帮助我们实现自动浏览站点和与其进行交互.它是依赖于其它python插件或拓展进行的,所以我们使用它之前需要安装一系列的依赖包 ...

  2. Linux内核调试方法总结之backtrace

    backtrace [用途]用户态或者内核态程序异常退出时回溯堆栈信息 [原理]通过对当前堆栈的分析,回溯上层函数在当前栈中的帧地址,直至顶层函数.帧地址是指在栈中存在局部变量.上一级函数返回地址.寄 ...

  3. 基于GTID模式MySQL主从复制

    基于GTID模式MySQL主从复制 GTID复制原理:基于GTID的复制是MySQL 5.6后新增的复制方式GTID (global transaction identifier) 即全局事务ID, ...

  4. 爬虫相关概念和https加密

    一.爬虫的相关概念 1.什么是爬虫 互联网:由网络设备(网线,路由器,交换机,防火墙)和一台台计算机连接而成,像一张网一样. 互联网建立目的:互联网的核心价值在与数据的共享/传递:数据是存放在一台台机 ...

  5. 接口自动化之ddt

    接口自动化会用到数据驱动模式,也就是一个ddt模块 目录 1.环境准备 2.调用时标准格式 3.应用(结合excle来传值) 1.环境准备 首先,需要安装ddt模块 pip install ddt 2 ...

  6. 接口自动化之get/post请求

    本篇旨在对get和post请求做一个总结 目录 1.get请求 2.post请求 1.get请求 get请求一般形式相同,仅带url即可发送请求 对于https协议的请求,加一个特殊处理即可(veri ...

  7. Java 基础-IO、stream 流、文件操作

    输入输出流的分类 在 java.io 包中,包含了输入输出操作所需的类. I/O 流可以安装不同的标准分类: 按照流的方向分类: 输入流:将信息从代码外部输入代码 输出流:将代码得到的数据输出到文件. ...

  8. Vue的生命周期(在其他地方看到一份非常好又详细的详解)

    链接地址:https://segmentfault.com/a/1190000011381906 首先,每个Vue实例在被创建之前都要经过一系列的初始化过程,这个过程就是vue的生命周期.首先看一张图 ...

  9. Nacos 配置中心原理分析

    我们从原生SDK代码中入手,可以发现最核心的两行代码: ConfigService configService=); 首先我们先来看 NacosFactory.createConfigService ...

  10. [BZOJ 2199] [USACO11JAN] 大陆议会The Continental Cowngress(2-SAT)

    [BZOJ 2199] [USACO11JAN] 大陆议会The Continental Cowngress(2-SAT) 题面 题面较长,略 分析 考虑把问题转化成一个依赖性问题 我们把每只奶牛投出 ...