SDOI2013森林
主席树启发式合并,每次连边维护并查集,集合大小,求lca所需信息,合并两个树上的主席树, 重点看代码.
#include <iostream>
#include <algorithm>
#include <cassert>
#include <cstdio>
#include <cstring>
using namespace std;
const int N = 8e4 + 1, M = 2e7 + 1;
int n, m, SIZE, Ht[N], a[N], T;
int head[N], nxt[N<<1], ver[N<<1], tot;
void add(int u, int v)
{ ver[++tot] = v, nxt[tot] = head[u], head[u] = tot; }
int f[N][21], dep[N];
int LCA(int u, int v) {
if (dep[u] < dep[v]) swap(u, v);
for (int i = 20; i >= 0; -- i) {
if (dep[f[u][i]] >= dep[v]) {
u = f[u][i];
}
}
assert(dep[u] == dep[v]);
if (u == v) return u;
for (int i = 20; i >= 0; -- i) {
if (f[u][i] != f[v][i]) {
u = f[u][i];
v = f[v][i];
}
}
assert(f[u][0] == f[v][0]);
return f[u][0];
}
int sum[M], ls[M], rs[M], cnt, rt[N];
void insert(int &p, int l, int r, int val) {
p = ++cnt;
if (l == r) return;
int mid = l + r >> 1;
if (val <= mid) insert(ls[p], l, mid, val);
else insert(rs[p], mid + 1, r, val);
}
void update(int last, int p, int l, int r, int val) {
assert(p);
sum[p] = sum[last] + 1;
if (l ==r) return;
int mid = l + r >> 1;
if (val <= mid) update(ls[last], ls[p], l, mid, val), rs[p] = rs[last];
else update(rs[last], rs[p], mid + 1, r, val), ls[p] = ls[last];
}
int query(int last1, int last2, int p1, int p2, int l, int r, int k) {
if (l == r) return Ht[l];
int s = sum[ls[p1]] - sum[ls[last1]] + sum[ls[p2]] - sum[ls[last2]];
int mid = l + r >> 1;
if (k <= s) return query(ls[last1], ls[last2], ls[p1], ls[p2], l, mid, k);
else return query(rs[last1], rs[last2], rs[p1], rs[p2], mid + 1, r, k - s);
}
bool vis[N];
void DFS(int u, int fa) {
vis[u] = 1;
update(rt[fa], rt[u], 1, SIZE, a[u]);
dep[u] = dep[fa] + 1;
f[u][0] = fa;
for (int i = 1; i <= 20; ++ i)
f[u][i] = f[ f[u][i - 1] ][i - 1];
for (int i = head[u]; i; i = nxt[i]) {
if (ver[i] != fa) DFS(ver[i], u);
}
}
int fa[N], size[N];
int find(int x)
{ return x == fa[x] ? x : fa[x] = find(fa[x]); }
void link(int &u, int &v) {
add(u, v); add(v, u);
int Fu = find(u), Fv = find(v);
if (size[Fu] > size[Fv]) swap(u, v), swap(Fu, Fv);
fa[Fv] = Fu; size[Fu] += size[Fv];
}
int calc(int u, int v, int k) {
int lca = LCA(u, v);
return query(rt[lca], rt[f[lca][0]], rt[u], rt[v], 1, SIZE, k);
}
int main() {
#ifndef ONLINE_JUDGE
freopen("P3302.in", "r", stdin);
freopen("P3302.out", "w", stdout);
#endif
int _;
scanf("%d %d %d %d", &_, &n, &m, &T);
for (int i = 1; i <= n; ++ i) scanf("%d", a + i), Ht[i] = a[i];
sort(Ht + 1, Ht + 1 + n);
SIZE = unique(Ht + 1 ,Ht + 1+ n) - Ht - 1;
for (int i = 1 ;i <= n; ++ i) a[i] = lower_bound(Ht + 1, Ht + 1 + SIZE, a[i]) - Ht;
for (int i = 1; i <= n; ++ i) fa[i] = i, size[i] = 1, insert(rt[i], 1, SIZE, a[i]);
for (int i = 1; i <= m; ++ i) {
int u, v;
scanf("%d %d", &u, &v);
link(u, v);
}
for (int i = 1; i <= n; ++ i)
if (!vis[find(i)]) DFS(find(i), 0);
int ans = 0;
char op[4];
while (T --) {
scanf("%s", op);
if (op[0] == 'Q') {
int x, y, k; scanf("%d%d%d", &x, &y, &k); x ^= ans, y ^= ans, k ^= ans; printf("%d\n", ans = calc(x, y, k));
} else {
int x, y; scanf("%d%d", &x, &y); x ^= ans, y ^= ans; link(x, y); DFS(x, y);
}
}
}
SDOI2013森林的更多相关文章
- BZOJ 3123: [Sdoi2013]森林 [主席树启发式合并]
3123: [Sdoi2013]森林 题意:一个森林,加边,询问路径上k小值.保证任意时刻是森林 LCT没法搞,树上kth肯定要用树上主席树 加边?启发式合并就好了,小的树dfs重建一下 注意 测试点 ...
- luoguP3302 [SDOI2013]森林 主席树 启发式合并
题目链接 luoguP3302 [SDOI2013]森林 题解 本来这题树上主席树暴力启发式合并就完了 结果把lca写错了... 以后再也不这么写了 复杂度\(O(nlog^2n)\) "f ...
- P3302 [SDOI2013]森林(主席树+启发式合并)
P3302 [SDOI2013]森林 主席树+启发式合并 (我以前的主席树板子是错的.......坑了我老久TAT) 第k小问题显然是主席树. 我们对每个点维护一棵包含其子树所有节点的主席树 询问(x ...
- [BZOJ3123][Sdoi2013]森林 主席树+启发式合并
3123: [Sdoi2013]森林 Time Limit: 20 Sec Memory Limit: 512 MB Description Input 第一行包含一个正整数testcase,表示当 ...
- BZOJ3123: [Sdoi2013]森林(启发式合并&主席树)
3123: [Sdoi2013]森林 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 4813 Solved: 1420[Submit][Status ...
- 【BZOJ3123】[Sdoi2013]森林 主席树+倍增LCA+启发式合并
[BZOJ3123][Sdoi2013]森林 Description Input 第一行包含一个正整数testcase,表示当前测试数据的测试点编号.保证1≤testcase≤20. 第二行包含三个整 ...
- 洛谷 P3302 [SDOI2013]森林 解题报告
P3302 [SDOI2013]森林 题目描述 小\(Z\)有一片森林,含有\(N\)个节点,每个节点上都有一个非负整数作为权值.初始的时候,森林中有\(M\)条边. 小Z希望执行\(T\)个操作,操 ...
- 3123: [Sdoi2013]森林
3123: [Sdoi2013]森林 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 3336 Solved: 978[Submit][Status] ...
- bzoj 3123: [Sdoi2013]森林(45分暴力)
3123: [Sdoi2013]森林 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 4184 Solved: 1235[Submit][Status ...
- AC日记——[Sdoi2013]森林 bzoj 3123
3123: [Sdoi2013]森林 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 3216 Solved: 944[Submit][Status] ...
随机推荐
- Jquery Mobile 常用参数
data-role 属性值: data-role参数表: page 页面容器,其内部的mobile元素将会继承这个容器上所设置的属性 header 页面标题容器,这个容器内部 ...
- Python 学习笔记(八)Python列表(二)
列表函数 追加和扩展 list.append() 在列表末尾追加新的对象 >>> dir(list) #dir 查看列表的函数 ['__add__', '__class__', '_ ...
- jzoj5196. 【NOIP2017提高组模拟7.3】B (数论,巧妙的暴力)
Description
- shell的命令格式
参考高峻峰 著 循序渐进Linux(第二版) command [options] [arguments] command:表示命令的名称 options:表示命令的选项 arguments:表示命令的 ...
- PostgreSQL异步主从流复制搭建
1 总体规划 Master库 Slave库 操作系统 CentOS Linux release 7.5.1804 CentOS Linux release 7.5.1804 处理器 1 1 内存 ...
- JAVA日期类图解
- Asp.Net Core使用Nginx实现反向代理
---恢复内容开始--- 前两篇文章介绍了使用Docker作为运行环境,用Nginx做反向代理的部署方法,这篇介绍一下使用Nginx配合.Net Core运行时直接在Liunx上裸奔的方法. 一.安装 ...
- mysql计算排名
mysql计算排名,获取行号rowno 学生成绩表数据 SELECT * FROM table_score ORDER BY score DESC; 获取某个学生成绩排名并计算该学生和上一名学生成绩差 ...
- JavaSE 第二次学习随笔(关于内存的小题)
class HelloA { public HelloA() { System.out.println("HelloA"); } { System.out.println(&quo ...
- vi-vim常用命令
vi-vim常用命令 1 简介 在UNIX系统中,创建和修改配置文件.shell脚本.初始化文件.编写程序都离不开VI. 1 vi[1]属于两个主要的UNIX规范:POSIX和单一UNIX规 ...