BZOJ 3123 SDOI2013 森林
首先对于查询操作就是裸的COT QAQ
在树上DFS建出主席树就可以了
对于连接操作,我们发现并没有删除
所以我们可以进行启发式合并,每次将小的树拍扁插入大的树里并重构即可
写完了之后第一个和第二个点迷のRE
然后又重新写了一遍就A了(并不知道为什么,难道第一遍写挫了?
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cstdlib>
using namespace std; const int maxn=100010;
int T,n,m,Q,u,v,k,lca,w;
int ans=0;
char ch;
int p[maxn];
int fa[maxn],dep[maxn],sz[maxn];
int anc[maxn][20];
int h[maxn],cnt=0;
struct edge{
int to,next;
}G[maxn<<1];
void add(int x,int y){++cnt;G[cnt].to=y;G[cnt].next=h[x];h[x]=cnt;}
int ufs(int x){return fa[x]==x?x:fa[x]=ufs(fa[x]);} int a[maxn];
struct val{
int v,id;
}c[maxn];
bool cmpv(const val &A,const val &B){return A.v<B.v;}
bool cmpid(const val &A,const val &B){return A.id<B.id;} int tot=0,rt[maxn];
struct Seg_Tree{
int L,R,v;
}t[11000010]; void read(int &num){
num=0;ch=getchar();
while(ch<'!')ch=getchar();
while(ch>='0'&&ch<='9')num=num*10+ch-'0',ch=getchar();
}
void build(int &o,int L,int R){
o=++tot;
if(L==R)return;
int mid=(L+R)>>1;
build(t[o].L,L,mid);
build(t[o].R,mid+1,R);
}
void modify(int &o,int now,int L,int R,int p){
t[++tot]=t[now];o=tot;
if(L==R){t[o].v++;return;}
int mid=(L+R)>>1;
if(p<=mid)modify(t[o].L,t[now].L,L,mid,p);
else modify(t[o].R,t[now].R,mid+1,R,p);
t[o].v=t[t[o].L].v+t[t[o].R].v;
}
void DFS(int u,int f){
anc[u][0]=f;
modify(rt[u],rt[f],1,m,c[u].v);
for(int i=1;i<=18;++i){
if(anc[u][i-1]!=0){
int a=anc[u][i-1];
anc[u][i]=anc[a][i-1];
}else anc[u][i]=0;
}
for(int i=h[u];i;i=G[i].next){
int v=G[i].to;
if(v==f)continue;
dep[v]=dep[u]+1;
DFS(v,u);
}return;
}
int LCA(int p,int q){
if(dep[p]<dep[q])swap(p,q);
int d=dep[p]-dep[q];
for(int i=18;i>=0;--i){
if(d>>i&1)p=anc[p][i];
}
if(p==q)return p;
for(int i=18;i>=0;--i){
if(anc[p][i]!=anc[q][i]&&anc[p][i]!=0)p=anc[p][i],q=anc[q][i];
}return anc[p][0];
}
int ask(int u,int v,int lca,int L,int R,int k){
if(L==R)return L;
int mid=(L+R)>>1;
int now=t[t[u].L].v+t[t[v].L].v-2*t[t[lca].L].v+(w<=mid&&w>=L);
if(k>now)return ask(t[u].R,t[v].R,t[lca].R,mid+1,R,k-now);
else return ask(t[u].L,t[v].L,t[lca].L,L,mid,k);
}
int main(){
read(T);read(n);read(m);read(Q);
for(int i=1;i<=n;++i)read(c[i].v),c[i].id=i,a[i]=c[i].v;
for(int i=1;i<=n;++i)fa[i]=i;
for(int i=1;i<=m;++i){
read(u);read(v);
add(u,v);add(v,u);
int d1=ufs(u),d2=ufs(v);
if(d1!=d2)fa[d1]=d2;
}m=0;
sort(c+1,c+n+1,cmpv);sort(a+1,a+n+1);
for(int i=1;i<=n;++i){
if(i==1||a[i]!=a[i-1])c[i].v=++m,p[m]=a[i];
else c[i].v=m;
}
sort(c+1,c+n+1,cmpid);
build(rt[0],1,m);
for(int i=1;i<=n;++i){
if(ufs(i)==i)DFS(i,0);
sz[ufs(i)]++;
}
while(Q--){
ch=getchar();
while(ch<'!')ch=getchar();
if(ch=='Q'){
scanf("%d%d%d",&u,&v,&k);
u^=ans;v^=ans;k^=ans;
lca=LCA(u,v);w=c[lca].v;
//cout<<u<<' '<<v<<' '<<k<<endl;
ans=p[ask(rt[u],rt[v],rt[lca],1,m,k)];
printf("%d\n",ans);
}else{
scanf("%d%d",&u,&v);
u^=ans;v^=ans;
//cout<<u<<' '<<v<<endl;
int d1=ufs(u),d2=ufs(v);
//cout<<d1<<' '<<d2<<endl;
if(sz[d1]<sz[d2])swap(u,v),swap(d1,d2);
sz[d1]+=sz[d2];
if(d1!=d2)fa[d2]=d1;//d2->d1
dep[v]=dep[u]+1;
DFS(v,u);add(u,v);add(v,u);
}
}return 0;
}
BZOJ 3123 SDOI2013 森林的更多相关文章
- BZOJ 3123: [Sdoi2013]森林 [主席树启发式合并]
3123: [Sdoi2013]森林 题意:一个森林,加边,询问路径上k小值.保证任意时刻是森林 LCT没法搞,树上kth肯定要用树上主席树 加边?启发式合并就好了,小的树dfs重建一下 注意 测试点 ...
- bzoj 3123: [Sdoi2013]森林(45分暴力)
3123: [Sdoi2013]森林 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 4184 Solved: 1235[Submit][Status ...
- Bzoj 3123: [Sdoi2013]森林(主席树+启发式合并)
3123: [Sdoi2013]森林 Time Limit: 20 Sec Memory Limit: 512 MB Description Input 第一行包含一个正整数testcase,表示当前 ...
- ●BZOJ 3123 [Sdoi2013]森林
题链: http://www.lydsy.com/JudgeOnline/problem.php?id=3123 题解: 主席树,在线,启发式合并 简单版(只有询问操作):[2588: Spoj 10 ...
- BZOJ 3123 [SDOI2013] 森林 - 启发式合并 主席树
Description 给你一片森林, 支持两个操作: 查询$x$到$y$的$K$大值, 连接两棵树中的两个点 Solution 对每个节点$x$动态开权值线段树, 表示从$x$到根节点路径上权值出 ...
- bzoj 3123 [Sdoi2013]森林(主席树,lca,启发式合并)
Description Input 第一行包含一个正整数testcase,表示当前测试数据的测试点编号.保证1≤testcase≤20. 第二行包含三个整数N,M,T,分别表示节点数.初始边数.操作数 ...
- bzoj 3123 [Sdoi2013]森林(主席树+启发式合并+LCA)
Description Input 第一行包含一个正整数testcase,表示当前测试数据的测试点编号.保证1≤testcase≤20. 第二行包含三个整数N,M,T,分别表示节点数.初始边数.操作数 ...
- 3123: [Sdoi2013]森林
3123: [Sdoi2013]森林 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 3336 Solved: 978[Submit][Status] ...
- 【BZOJ】3123: [Sdoi2013]森林
题解 ------------------ 我莫不是一个智障吧 我把testdata的编号 当成数据组数读进来 我简直有毒 以为哪里写错了自闭了好久 实际上这题很简单,只要愉悦地开个启发式合并,然后每 ...
随机推荐
- (转)设置Win7防火墙规则 顺畅访问局域网
在Windows 7系统的电脑上搭建WAMP环境后,发现在局域网中其他电脑不能访问.有朋友告诉小强,这可能是因为当时Windows 7自带的防火墙屏蔽了80端口,只需要重新设置规则就可以了. 点击Wi ...
- (转) DockPanel 右键增加关闭,除此之外全部关闭的功能
在项目中新建一个class文件,代码如下: using System; using System.Collections.Generic; using System.ComponentModel; u ...
- VS2013编译GLUI
vs自带的OpenGL为1.1版本,太老了. 1,编译glut https://www.opengl.org/resources/libraries/glut/glut37.zip 查看生成路径,可以 ...
- C++与Lua交互(四)
引言 通过前几篇,我们已经对Lua的C API有了一定的了解,如lua_push*.lua_is*.lua_to*等等.用C++调用Lua数据时,我们主要运用lua_getglobal与lua_pus ...
- 小黑的镇魂曲(HDU2155:贪心+dfs+奇葩解法)
题目:点这里 题目的意思跟所谓的是英雄就下100层一个意思……在T秒内能够下到地面,就可以了(还有一个板与板之间不能超过H高). 接触这题目是在昨晚的训练赛,当时拍拍地打了个贪心+dfs,果断跟我想的 ...
- vc 编译执行bat
转载:
- xamarin.ios 豆瓣电台视频教程
视频中提到的网址: http://www.sufeinet.com/thread-655-1-1.html https://github.com/akfish/fm-terminal/blob/dev ...
- 网站如何做到完全不需要jQuery
jQuery是现在最流行的JavaScript工具库. 据统计,目前全世界57.3%的网站使用它.也就是说,10个网站里面,有6个使用jQuery.如果只考察使用工具库的网站,这个比例就会上升到惊人的 ...
- React-router 要点
1.关于url中传参的问题 比如我想打开: /articles/detail/101 在url中要传一个参数 /articles/detail/:articleId 路由中:<Route pat ...
- jquery获得option的值(示例)
jquery获得option的值和对option的操作. jQuery获取Select元素,并选择的Text和Value: 复制代码代码如下: $("#select_id").ch ...