#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#define maxn 100005
#define maxm 400005
using namespace std; int n,m,ans,vis[maxn],head,tot,now[maxn],prep[maxm],son_[maxm],fa[maxn],son[maxn][],size[maxn],add[maxn],list[maxn],tail;
bool rev[maxn]; bool which(int x){
return son[fa[x]][]==x;
} bool isroot(int x){
return son[fa[x]][]!=x&&son[fa[x]][]!=x;
} void insert_(int x,int y){
size[x]+=y,add[x]+=y;
} void pushdown(int x){
if (!x) return;
if (rev[x]){
rev[x]^=,swap(son[x][],son[x][]);
if (son[x][]) rev[son[x][]]^=;
if (son[x][]) rev[son[x][]]^=;
}
if (add[x]){
if (son[x][]) insert_(son[x][],add[x]);
if (son[x][]) insert_(son[x][],add[x]);
add[x]=;
}
} void relax(int x){
if (!isroot(x)) relax(fa[x]);
pushdown(x);
} void rotata(int x){
int y=fa[x],d=which(x),dd=which(y);
if (!isroot(y)) son[fa[y]][dd]=x; fa[x]=fa[y];
fa[son[x][d^]]=y,son[y][d]=son[x][d^];
fa[y]=x,son[x][d^]=y;
} void splay(int x){
relax(x);
while (!isroot(x)){
if (isroot(fa[x])) rotata(x);
else if (which(x)==which(fa[x])) rotata(fa[x]),rotata(x);
else rotata(x),rotata(x);
}
} void access(int x){
for (int p=;x;x=fa[x]){
splay(x);
son[x][]=p;
p=x;
}
} int succ(int x){
splay(x);
int y=son[x][];
for (pushdown(y);son[y][];y=son[y][],pushdown(y));
return y;
} int find_root(int x){
access(x);
splay(x);
while (son[x][]) x=son[x][];
return x;
} void insert(int x,int y){
tot++,prep[tot]=now[x],now[x]=tot,son_[tot]=y;
} void bfs(int x,int y){
int t;
head=,tail=,list[]=y,size[y]=,fa[y]=son[y][]=son[y][]=,vis[y]=m,add[y]=rev[y]=;
while (head<tail){
head++,t=list[head];
for (int i=now[t],so=son_[i];i;i=prep[i],so=son_[i]){
if (vis[so]!=m){
vis[so]=m,size[so]=,fa[so]=son[so][]=son[so][]=add[so]=rev[so]=,fa[so]=t;
list[++tail]=so;
}
}
}
for (int i=tail;i>=;i--){
t=list[i];
if (fa[t]) size[fa[t]]+=size[t];
}
} void make_root(int x){
access(x);
splay(x);
rev[x]^=;
} void merge(int x,int y){
int rx,ry,t;
rx=find_root(x),ry=find_root(y);
if (rx==ry) return;
if (size[rx]<size[ry]) swap(x,y),swap(rx,ry);
ans^=ry;
bfs(x,y),insert(x,y),insert(y,x),fa[y]=x;
access(x),splay(x),insert_(x,size[y]);
// access(ry),splay(rx),relax(ry);
access(ry);
for (;;){
t=succ(rx);
if (!t){
make_root(rx);
break;
}
if (size[t]*>size[rx]||(size[t]*==size[rx]&&t<rx)){
// son[t][0]=0;
swap(size[t],size[rx]),size[rx]=size[t]-size[rx];
ans^=rx,ans^=t,rx=t;
}else{
make_root(rx);
break;
}
}
} int main(){
freopen("capital.in","r",stdin);
freopen("capital.out","w",stdout);
char st[];
int u,v;
memset(vis,-,sizeof(vis));
scanf("%d%d",&n,&m),ans=;
tot=,memset(now,,sizeof(now));
memset(rev,,sizeof(rev));
memset(size,,sizeof(size));
for (int i=;i<=n;i++) add[i]=fa[i]=son[i][]=son[i][]=,size[i]=;
for (int i=;i<=n;i++) ans^=i;
while (m--){
scanf("%s",st+);
if (st[]=='X') printf("%d\n",ans);
else if (st[]=='Q') scanf("%d",&u),printf("%d\n",find_root(u));
else scanf("%d%d",&u,&v),merge(u,v);
}
return ;
}

题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=3510

题意:一个森林,要求维护每棵树的重心(多个重心取编号最小的那个),支持三个操作:

1.          连接两个节点(保证两个节点不在同一棵树中)。

2.          询问某个节点所在树的重心。

3.          询问所有树的重心编号的异或和。

关于重心有几个性质:

  1. 以重心为根,任意一个子树的节点个数<=总结点数/2,否则重心可向该子树方向转移。
  2. 设树中节点x的权值为该点到树中每个点的距离之和,则重心的权值最小。
  3. 两个树合并为一颗新树后,新树的重心一定在原来的两个重心之间唯一的路径上。

做法:考虑用lct动态维护树的重心,以原树中的重心做为原树根,

由于只有合并操作(合并节点x和y,sizex>sizey),可以考虑启发式合并,即把节点数少的那棵树(y所在树)暴力重建(一次BFS),接到另一棵树(x所在树)上,节点y所在树上的信息可在bfs过程中维护好,节点x所在的树的size信息会有什么变化呢,我们发现节点x到原x树重心的路径上size[i]会加上y树的总结点数,这就是用lct的缘故。接下来我们考虑如何维护新树的重心,并作为新树的根呢?根据重心的性质3,我们考虑暴力转移重心,由于是启发式合并,据说复杂度可靠,我们将原来的两个重心的路径放到同一棵splay中来,从x树中重心开始暴力转移,如果节点p的后继节点q,若(sizeq*2>sizep)or(sizeq*2=sizep&&q<p),则将重心从p转移至q点,最后将make_root(新重心),操作1得以解决。

对于操作2,用lct求该节点所在树的根即可,因为树根为重心。

对于操作3, 用ans维护即可,每次更新重心更新时更新ans。

Lct动态维护树的重心(Lct+启发式合并+暴力转移重心)。

Bzoj3510首都的更多相关文章

  1. BZOJ3510 首都(LCT)

    即动态维护树的重心.考虑合并后的新重心一定在两棵树的重心的连线上.于是对每个点维护其子树大小,合并时在这条链的splay上二分即可.至于如何维护子树大小,见https://blog.csdn.net/ ...

  2. BZOJ3510 首都

    题目描述 在X星球上有N个国家,每个国家占据着X星球的一座城市.由于国家之间是敌对关系,所以不同国家的两个城市是不会有公路相连的. X星球上战乱频发,如果A国打败了B国,那么B国将永远从这个星球消失, ...

  3. BZOJ3510首都(LCT)

    Description 在X星球上有N个国家,每个国家占据着X星球的一座城市.由于国家之间是敌对关系,所以不同国家的两个城市是不会有公路相连的. X星球上战乱频发,如果A国打败了B国,那么B国将永远从 ...

  4. bzoj3510 首都 LCT 维护子树信息+树的重心

    题目传送门 https://lydsy.com/JudgeOnline/problem.php?id=3510 题解 首先每一个连通块的首都根据定义,显然就是直径. 然后考虑直径的几个性质: 定义:删 ...

  5. Some Conclusions.

    目录 DP 四边形不等式 数论 & 数学 数据结构 树链剖分 左偏树的性质及\(O(n)\)的构造 图论 树 二分图 竞赛图 平面图 双连通分量 字符串 后缀自动机 复杂度分析 没什么好写的. ...

  6. LCT维护子树信息

    有些题目,在要求支持link-cut之外,还会在线询问某个子树的信息.LCT可以通过维护虚边信息完成这个操作. 对于LCT上每个节点,维护两个两sz和si,后者维护该点所有虚儿子的信息,前者维护该点的 ...

  7. 【BZOJ3510】首都 LCT维护子树信息+启发式合并

    [BZOJ3510]首都 Description 在X星球上有N个国家,每个国家占据着X星球的一座城市.由于国家之间是敌对关系,所以不同国家的两个城市是不会有公路相连的. X星球上战乱频发,如果A国打 ...

  8. 洛谷P4299 首都(BZOJ3510)(LCT,树的重心,二分查找)

    Update:原来的洛谷U21715已成坑qwq 已经被某位管理员巨佬放进公共题库啦!又可以多一个AC记录啦! 洛谷题目传送门 其实也可以到这里交啦 思路分析 动态维护树的重心 题目中说到国家的首都会 ...

  9. 【bzoj3510】首都 LCT维护子树信息(+启发式合并)

    题目描述 在X星球上有N个国家,每个国家占据着X星球的一座城市.由于国家之间是敌对关系,所以不同国家的两个城市是不会有公路相连的. X星球上战乱频发,如果A国打败了B国,那么B国将永远从这个星球消失, ...

随机推荐

  1. Android手机浏览器访问本地网络相关问题

    为了测试开发的手机网站,常常需要使手机直接访问本地网络. 在这个过程中碰到几个问题,记下来供以后参考 1. 在本地主机运行apache后,使用localhost和127.0.0.1可以访问页面,但使用 ...

  2. BZOJ 1087 【SCOI2005】 互不侵犯King

    Description 在N×N的棋盘里面放K个国王,使他们互不攻击,共有多少种摆放方案.国王能攻击到它上下左右,以及左上左下右上右下八个方向上附近的各一个格子,共8个格子. Input 只有一行,包 ...

  3. [转]各种有用的PHP开源库精心收集

    FROM : http://my.oschina.net/caroltc/blog/324024 1.html2ps and html2pdf    下载地址: http://www.tufat.co ...

  4. SQLServer 分布式查询MySQL

    这学期开了分布式数据库这门课,开始编程实现,今天调试了一早上,写下此配置文件方便查询. 本文实现的是SQLServer2008 Express 链式添加MySql-5.6.10数据库,进行远程操作. ...

  5. 分布式中使用Redis实现Session共享(一)

    上一篇介绍了如何使用nginx+iis部署一个简单的分布式系统,文章结尾留下了几个问题,其中一个是"如何解决多站点下Session共享".这篇文章将会介绍如何使用Redis,下一篇 ...

  6. 单例模式的两种实现方式对比:DCL (double check idiom)双重检查 和 lazy initialization holder class(静态内部类)

    首先这两种方式都是延迟初始化机制,就是当要用到的时候再去初始化. 但是Effective Java书中说过:除非绝对必要,否则就不要这么做. 1. DCL (double checked lockin ...

  7. lecture5-对象识别与卷积神经网络

    Hinton第五课 突然不知道object recognition 该翻译成对象识别好,还是目标识别好,还是物体识别好,但是鉴于范围性,还是翻译成对象识别吧.这一课附带了两个论文<Convolu ...

  8. Python2.7-异常和工具

    来自<python学习手册第四版>第七部分,而且本书发布的时候3.1还未发布,所以针对本书的一些知识会有些滞后于python的版本,具体更多细节可以参考python的标准手册. 一.异常基 ...

  9. 是什么时候开始学习gulp了

    转自:http://www.ydcss.com/archives/18 简介: gulp是前端开发过程中对代码进行构建的工具,是自动化项目的构建利器:她不仅能对网站资源进行优化,而且在开发过程中很多重 ...

  10. Chrome 监听 console 打开

    这个算是 Chrome only 其他的我没测试,也不想测试.因为我的控制台脚本仅仅在 Chrome 下加载. 如果你需要全平台,那么这肯定不是你需要的结果. 需求 其实我很早就想折腾这个了,但是,, ...