[SDOI2013]森林 主席树+启发式合并
这题的想法真的很妙啊。
看到题的第一眼,我先想到树链剖分,并把\(DFS\)序当成一段区间上主席树。但是会发现在询问的时候,可能会非常复杂,因为你需要把路径拆成很多条轻链和重链,它们还不一定连续,很难做(这个做法貌似可以用于子树第\(k\)大问题)。
于是我们换一个思路,让某个点的从它的父亲继承信息,也就是让这个结点对应的主席树维护一条从它到根的链。
那查询该如何是好?假设要查询的结点为\(x,y\),我们利用树上差分的思想让\(x,y,lca(x,y),fa[lca(x,y)]\)对应的四棵主席树一起向下跳就行了。
有连边操作时,就要用到启发式合并了,也就是把小的那棵树合并到大的中,暴力重构一下小树的倍增数组和主席树就行了。
代码:
#include <bits/stdc++.h>
using namespace std;
#define N 100000
#define pb push_back
#define mp make_pair
#define mid ((l+r)>>1)
#define pii pair<int, int>
#define D 17
int n0, n, m, q, T, root[N+5], sz[N+5], d[N+5], me[N+5];
int w0[N+5], w[N+5], as[N+5], f[N+5][20];
int nid, sumv[100*N+5], lson[100*N+5], rson[100*N+5];
int lastans;
vector<int> G[N+5];
void build(int &o, int l, int r)
{
o = ++nid;
if(l == r) return ;
build(lson[o], l, mid), build(rson[o], mid+1, r);
}
int LCA(int x, int y)
{
if(d[x] < d[y]) swap(x, y);
for(int i = D; i >= 0; --i)
if(d[f[x][i]] >= d[y]) x = f[x][i];
if(x == y) return x;
for(int i = D; i >= 0 && f[x][0] != f[y][0]; --i)
if(f[x][i] != f[y][i]) x = f[x][i], y = f[y][i];
return f[x][0];
}
void insert(int o, int &u, int l, int r, int x, int k)
{
u = ++nid;
lson[u] = lson[o], rson[u] = rson[o], sumv[u] = sumv[o]+k;
if(l == r) return ;
if(x <= mid) insert(lson[o], lson[u], l, mid, x, k);
else insert(rson[o], rson[u], mid+1, r, x, k);
}
int query(int x, int y, int k) //x-->y 第k大
{
int lca = LCA(x, y), fa = f[lca][0], l = 1, r = n0;
x = root[x], y = root[y], lca = root[lca], fa = root[fa];
while(l != r)
{
int t = sumv[lson[x]]+sumv[lson[y]]-sumv[lson[fa]]-sumv[lson[lca]];
if(t >= k) r = mid, x = lson[x], y = lson[y], fa = lson[fa], lca = lson[lca];
else k -= t, l = mid+1, x = rson[x], y = rson[y], fa = rson[fa], lca = rson[lca];
}
return as[l];
}
void dfs(int u, int fa, int anc)
{
d[u] = d[fa]+1, sz[u] = 1, f[u][0] = fa, me[u] = anc;
for(int i = 1; i <= D; ++i) f[u][i] = f[f[u][i-1]][i-1];
insert(root[fa], root[u], 1, n0, w[u], 1);
for(int i = 0, v; i < G[u].size(); ++i)
{
v = G[u][i];
if(v == fa) continue;
dfs(v, u, anc);
sz[u] += sz[v];
}
}
void link(int x, int y)
{
G[x].pb(y), G[y].pb(x);
if(sz[me[x]] < sz[me[y]]) dfs(x, y, me[y]);
else dfs(y, x, me[x]);
}
int main()
{
#ifndef ONLINE_JUDGE
freopen("testdata.in", "r", stdin);
freopen("testdata.out", "w", stdout);
#endif
scanf("%d", &T);
scanf("%d%d%d", &n, &m, &q);
for(int i = 1; i <= n; ++i) scanf("%d", &w0[i]), w[i] = w0[i]; //离散化
sort(w0+1, w0+n+1);
n0 = unique(w0+1, w0+n+1)-w0-1;
for(int i = 1, t; i <= n; ++i)
{
t = lower_bound(w0+1, w0+n0+1, w[i])-w0;
as[t] = w[i], w[i] = t;
}
for(int i = 1, x, y; i <= m; ++i)
{
scanf("%d%d", &x, &y);
G[x].pb(y), G[y].pb(x);
}
build(root[0], 1, n0);
for(int i = 1; i <= n; ++i)
if(!d[i]) dfs(i, 0, i);
char c;
for(int i = 1, x, y, z; i <= q; ++i)
{
cin >> c;
scanf("%d%d", &x, &y);
x ^= lastans, y ^= lastans;
if(c == 'Q') scanf("%d", &z), z ^= lastans, printf("%d\n", lastans = query(x, y, z));
else link(x, y);
}
return 0;
}
[SDOI2013]森林 主席树+启发式合并的更多相关文章
- Bzoj 3123: [Sdoi2013]森林(主席树+启发式合并)
3123: [Sdoi2013]森林 Time Limit: 20 Sec Memory Limit: 512 MB Description Input 第一行包含一个正整数testcase,表示当前 ...
- [bzoj3123] [SDOI2013]森林 主席树+启发式合并+LCT
Description Input 第一行包含一个正整数testcase,表示当前测试数据的测试点编号.保证1≤testcase≤20. 第二行包含三个整数N,M,T,分别表示节点数.初始边数.操作数 ...
- BZOJ 3123: [Sdoi2013]森林 [主席树启发式合并]
3123: [Sdoi2013]森林 题意:一个森林,加边,询问路径上k小值.保证任意时刻是森林 LCT没法搞,树上kth肯定要用树上主席树 加边?启发式合并就好了,小的树dfs重建一下 注意 测试点 ...
- luoguP3302 [SDOI2013]森林 主席树 启发式合并
题目链接 luoguP3302 [SDOI2013]森林 题解 本来这题树上主席树暴力启发式合并就完了 结果把lca写错了... 以后再也不这么写了 复杂度\(O(nlog^2n)\) "f ...
- [BZOJ3123][Sdoi2013]森林 主席树+启发式合并
3123: [Sdoi2013]森林 Time Limit: 20 Sec Memory Limit: 512 MB Description Input 第一行包含一个正整数testcase,表示当 ...
- 【BZOJ 3123】 [Sdoi2013]森林 主席树启发式合并
我们直接按父子关系建主席树,然后记录倍增方便以后求LCA,同时用并查集维护根节点,而且还要记录根节点对应的size,用来对其启发式合并,然后每当我们合并的时候我们都要暴力拆小的一部分重复以上部分,总时 ...
- 【BZOJ-3123】森林 主席树 + 启发式合并
3123: [Sdoi2013]森林 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 2738 Solved: 806[Submit][Status] ...
- P3302 [SDOI2013]森林(主席树+启发式合并)
P3302 [SDOI2013]森林 主席树+启发式合并 (我以前的主席树板子是错的.......坑了我老久TAT) 第k小问题显然是主席树. 我们对每个点维护一棵包含其子树所有节点的主席树 询问(x ...
- BZOJ_3123_[Sdoi2013]森林_主席树+启发式合并
BZOJ_3123_[Sdoi2013]森林_主席树+启发式合并 Description Input 第一行包含一个正整数testcase,表示当前测试数据的测试点编号.保证1≤testcase≤20 ...
随机推荐
- Spring MVC(四)文件上传
文件上传步骤 1.写一个文件上传的页面 2.写一个文件上传的控制器 注意: 1.method="post" 2.enctype="multipart/form-data& ...
- animate-queue和step-animate
Step-animate: 分为3部分:{配置},{step:function(){...},duration:1000} <div id="warpper" style=& ...
- 生鲜配送管理系统_升鲜宝V2.0 价格组功能 操作说明_15382353715
价格组功能是B端供应链系统,必不可少的一个功能,其主要实现不同的客户不同的价格,B端系统有一个最大的不同就是,有些商品后台下单人员能看到的.有些商品在销售的那一瞬间,还不知道价格.所以这些商品只有后台 ...
- 【Android】用Cubism 2制作自己的Live2D——android sdk样本的下载与Android studio编译!
前言- 在浏览Live2d说明书的时候我无意中发现了一个有趣的东西,就是android sdk中居然自带动态壁纸!那就让我们来试试吧,说明书此页的网址连接——中文版||日文版 Android开发所必需 ...
- event 和delegate的分别
突然想起delegate委托是支持+= 和-=操作的,然后研究一下究竟这个是怎么做到的,好模仿一下.一开始以为是+=的运算符重载,但是在类库参考中并没有这个运算符重载,只有!= 和==运算符重载.有点 ...
- Ngnix负载均衡安装及配置
1.ngnix概念 Nginx是一款高性能的http 服务器/反向代理服务器及电子邮件(IMAP/POP3)代理服务器.由俄罗斯的程序设计师Igor Sysoev所开发,官方测试nginx能够支支撑5 ...
- 导入python库失败时的方法
出现以下错误如何解决: e.g. cmd: pip install Django -i http://mirrors.aliyun.com/pypi/simple/ --trusted-host ...
- 浏览器仿EXCEL表格插件 版本更新 - 智表ZCELL产品V1.3.1更新
智表(zcell)是一款浏览器仿excel表格jquery插件.智表可以为你提供excel般的智能体验,支持双击编辑.设置公式.设置显示小数精度.下拉框.自定义单元格.复制粘贴.不连续选定.合并单元格 ...
- java.util.Arrays.useLegacyMergeSort=true 作用
(原) 今天看了一下现场的环境,发现有个其它部门的项目用到了这样一个参数: -Djava.util.Arrays.useLegacyMergeSort=true 于是查看了一下什么作用. 在JDK1. ...
- CentOS 安装 ceph 单机版(luminous版本)
一.环境准备 CentOS Linux release 7.4.1708 (Core)一台,4块磁盘(sda.sdb,.sdc.sdd) 192.168.27.130 nceph 二.配置环境 1.修 ...