[bzoj3123][洛谷P3302] [SDOI2013]森林(树上主席树+启发式合并)
突然发现好像没有那么难……https://blog.csdn.net/stone41123/article/details/78167288
首先有两个操作,一个查询,一个连接
查询的话,直接在树上建主席树
然后难点在于连接
用启发式合并就可以了(想了半天都没想出来)
每次合并时,我们把小的树接到大的上,然后dfs一遍小的树,更新信息
然后注意数组……别太小也别太大……(被数组大小坑了好几次提交)
//minamoto
#include<bits/stdc++.h>
using namespace std;
#define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
char buf[<<],*p1=buf,*p2=buf;
inline int read(){
#define num ch-'0'
char ch;bool flag=;int res;
while(!isdigit(ch=getc()))
(ch=='-')&&(flag=true);
for(res=num;isdigit(ch=getc());res=res*+num);
(flag)&&(res=-res);
#undef num
return res;
}
char obuf[<<],*o=obuf;
void print(int x){
if(x>) print(x/);
*o++=x%+;
}
const int N=,M=N*;
int ver[N<<],Next[N<<],head[N];
int a[N],fa[N],sz[N],b[N];
int n,m,tot,q,size,ans;
void add(int u,int v){
ver[++tot]=v,Next[tot]=head[u],head[u]=tot;
ver[++tot]=u,Next[tot]=head[v],head[v]=tot;
}
int L[M],R[M],sum[M],rt[N],cnt;
void update(int last,int &now,int l,int r,int x){
sum[now=++cnt]=sum[last]+;
if(l==r) return;
int mid=(l+r)>>;
if(x<=mid) R[now]=R[last],update(L[last],L[now],l,mid,x);
else L[now]=L[last],update(R[last],R[now],mid+,r,x);
}
int query(int u,int v,int lca,int lca_fa,int l,int r,int k){
if(l>=r) return l;
int x=sum[L[v]]+sum[L[u]]-sum[L[lca]]-sum[L[lca_fa]];
int mid=(l+r)>>;
if(x>=k) return query(L[u],L[v],L[lca],L[lca_fa],l,mid,k);
else return query(R[u],R[v],R[lca],R[lca_fa],mid+,r,k-x);
}
inline int hash(int x){
return lower_bound(b+,b++size,x)-b;
}
int ff(int x){
return fa[x]==x?x:fa[x]=ff(fa[x]);
}
int st[N][],d[N],vis[N];
void dfs(int u,int father,int root){
st[u][]=father;
for(int i=;i<=;++i)
st[u][i]=st[st[u][i-]][i-];
++sz[root];
d[u]=d[father]+;
fa[u]=root;
vis[u]=;
update(rt[father],rt[u],,size,hash(a[u]));
for(int i=head[u];i;i=Next[i]){
int v=ver[i];
if(v==father) continue;
dfs(v,u,root);
}
}
int LCA(int x,int y){
if(x==y) return x;
if(d[x]<d[y]) swap(x,y);
for(int i=;i>=;--i){
if(d[st[x][i]]>=d[y]) x=st[x][i];
}
if(x==y) return x;
for(int i=;i>=;--i){
if(st[x][i]!=st[y][i])
x=st[x][i],y=st[y][i];
}
return st[x][];
}
int main(){
//freopen("testdata.in","r",stdin);
int t=read();
n=read(),m=read(),q=read();
for(int i=;i<=n;++i)
a[i]=b[i]=read(),fa[i]=i;
sort(b+,b++n);
size=unique(b+,b++n)-b-;
for(int i=;i<=m;++i){
int u=read(),v=read();
add(u,v);
}
for(int i=;i<=n;++i)
if(!vis[i]) dfs(i,,i);
while(q--){
char ch;int x,y;
while(!isupper(ch=getc()));
x=read()^ans,y=read()^ans;
if(ch=='Q'){
int k=read()^ans;
int lca=LCA(x,y);
ans=b[query(rt[x],rt[y],rt[lca],rt[st[lca][]],,size,k)];
print(ans),*o++='\n';
}
else{
add(x,y);
int u=ff(x),v=ff(y);
if(sz[u]<sz[v]) swap(x,y),swap(u,v);
dfs(y,x,u);
}
}
fwrite(obuf,o-obuf,,stdout);
return ;
}
[bzoj3123][洛谷P3302] [SDOI2013]森林(树上主席树+启发式合并)的更多相关文章
- P3302 [SDOI2013]森林(主席树+启发式合并)
P3302 [SDOI2013]森林 主席树+启发式合并 (我以前的主席树板子是错的.......坑了我老久TAT) 第k小问题显然是主席树. 我们对每个点维护一棵包含其子树所有节点的主席树 询问(x ...
- bzoj3123 [Sdoi2013]森林 树上主席树+启发式合并
题目传送门 https://lydsy.com/JudgeOnline/problem.php?id=3123 题解 如果是静态的查询操作,那么就是直接树上主席树的板子. 但是我们现在有了一个连接两棵 ...
- BZOJ_3123_[Sdoi2013]森林_主席树+启发式合并
BZOJ_3123_[Sdoi2013]森林_主席树+启发式合并 Description Input 第一行包含一个正整数testcase,表示当前测试数据的测试点编号.保证1≤testcase≤20 ...
- BZOJ2123 [Sdoi2013]森林 【主席树 + 启发式合并】
题目 输入格式 第一行包含一个正整数testcase,表示当前测试数据的测试点编号.保证1≤testcase≤20. 第二行包含三个整数N,M,T,分别表示节点数.初始边数.操作数.第三行包含N个非负 ...
- bzoj 3123 [Sdoi2013]森林(主席树+启发式合并+LCA)
Description Input 第一行包含一个正整数testcase,表示当前测试数据的测试点编号.保证1≤testcase≤20. 第二行包含三个整数N,M,T,分别表示节点数.初始边数.操作数 ...
- [luoguP3302] [SDOI2013]森林(主席树 + 启发式合并 + lca)
传送门 显然树上第k大直接主席树 如果连边的话,我们重构小的那一棵,连到另一棵上. 说起来简单,调了我一晚上. 总的来说3个错误: 1.离散化写错位置写到了后面 2."="写成了& ...
- 洛谷 P3302 [SDOI2013]森林 解题报告
P3302 [SDOI2013]森林 题目描述 小\(Z\)有一片森林,含有\(N\)个节点,每个节点上都有一个非负整数作为权值.初始的时候,森林中有\(M\)条边. 小Z希望执行\(T\)个操作,操 ...
- p3302 [SDOI2013]森林(树上主席树+启发式合并)
对着题目yy了一天加上看了一中午题解,终于搞明白了我太弱了 连边就是合并线段树,把小的集合合并到大的上,可以保证规模至少增加一半,复杂度可以是\(O(logn)\) 合并的时候暴力dfs修改倍增数组和 ...
- 洛谷 P3302 [SDOI2013]森林 Lebal:主席树 + 启发式合并 + LCA
题目描述 小Z有一片森林,含有N个节点,每个节点上都有一个非负整数作为权值.初始的时候,森林中有M条边. 小Z希望执行T个操作,操作有两类: Q x y k查询点x到点y路径上所有的权值中,第k小的权 ...
随机推荐
- node.js安装后出现环境变量错误找不到node
安装node.js和bower之后,运行bower出现/usr/bin/env: 'node': No such file or directory错误 这个错误是由于安装完node.js环境变量并没 ...
- 洛谷P1044 栈
之前看这题还是一头雾水,现在看:啊啊啊lydnb! 思考了一段时间,发现可以用DP. 令f[i]表示有i辆车时的方案数. 我一开始考虑的是在后面加车,可是这样搞不出状态转移方程来. 然后我考虑从前面加 ...
- A1102. Invert a Binary Tree
The following is from Max Howell @twitter: Google: 90% of our engineers use the software you wrote ( ...
- 跟我一起学习vue2(使用命令行搭建单页应用)[二]
第一步:运行git命令,全局安装 vue-cli $ cnpm install --global vue-cli 第二步: 创建一个基于 webpack 模板的新项目 $ vue init webpa ...
- H5新特性总结
Web前端最新的官方标准无疑就是es6和H5了,es6之前已经说过了(多是一些语法糖),现在来总结一下H5给我们来了哪些新“朋友”吧~ 1.video/radio 视频/音频 2.canvas 绘画 ...
- Nginx的alias的用法及与root的区别
以前只知道Nginx的location块中的root用法,用起来总是感觉满足不了自己的一些想法.然后终于发现了alias这个东西. 先看toot的用法 location /request_path/i ...
- springboot websocket 一篇足够了
什么是WebSocket WebSocket是一种在单个TCP连接上进行全双工通信的协议 … 为什么要实现握手监控管理 如果说,连接随意创建,不管的话,会存在错误,broken pipe 表面看单纯报 ...
- .Net进阶系列(15)-异步多线程(线程的特殊处理和深究委托赋值)(被替换)
1. 线程的异常处理 我们经常会遇到一个场景,开启了多个线程,其中一个线程报错,导致整个程序崩溃.这并不是我们想要的,我需要的结果是,其中一个线程报错,默默的记录下,其它线程正常进行,保证程序整体可以 ...
- CSS魔法(二)
# 文档类型<!DOCTYPE> <!DOCTYPE html> # 字符集 <meta charset="UTF-8" /> # 换行标签 & ...
- Syncfusion HTMLUI研究一
HTMLUI可以加载HTML页面,并且相比WebKit等占用资源特别少 <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitiona ...