题目链接

LCT(良心总结)

  1. #include <cstdio>
  2. #include <cctype>
  3. #include <algorithm>
  4. #define gc() getchar()
  5. const int N=3e5+5;
  7. inline int read()
  8. {
  9. 	int now=0;register char c=gc();
  10. 	for(;!isdigit(c);c=gc());
  11. 	for(;isdigit(c);now=now*10+c-'0',c=gc());
  12. 	return now;
  13. }
  14. namespace LCT
  15. {
  16. 	#define lson son[x][0]
  17. 	#define rson son[x][1]
  19. 	int fa[N],val[N],sum[N],son[N][2],sk[N];
  20. 	bool tag[N];
  21. 	inline void Update(int x){
  22. 		sum[x]=sum[lson]^sum[rson]^val[x];
  23. 	}
  24. 	inline void Rev(int x){
  25. 		std::swap(lson,rson), tag[x]^=1;
  26. 	}
  27. 	inline void PushDown(int x){
  28. 		if(tag[x]) Rev(lson),Rev(rson),tag[x]=0;
  29. 	}
  30. 	inline bool n_root(int x){
  31. 		return son[fa[x]][0]==x||son[fa[x]][1]==x;
  32. 	}
  33. 	void Rotate(int x)
  34. 	{
  35. 		int a=fa[x],b=fa[a],l=son[a][1]==x,r=l^1;
  36. 		if(n_root(a)) son[b][son[b][1]==a]=x;
  37. 		if(son[x][r]) fa[son[x][r]]=a;
  38. 		fa[a]=x, fa[x]=b;
  39. 		son[a][l]=son[x][r], son[x][r]=a;
  40. 		Update(a), Update(x);
  41. 	}
  42. 	void Splay(int x)
  43. 	{
  44. 		int t=1,a=x,b; sk[1]=x;
  45. 		while(n_root(a)) sk[++t]=a=fa[a];//别写成while(fa[x])!
  46. 		while(t) PushDown(sk[t--]);
  47. 		while(n_root(x))
  48. 		{
  49. 			a=fa[x], b=fa[a];
  50. 			if(n_root(a)) Rotate(son[a][1]==x^son[b][1]==a?x:a);
  51. 			Rotate(x);
  52. 		}
  53. 		Update(x);//还是加上吧
  54. 	}
  55. 	void Access(int x){
  56. 		for(int pre=0; x; x=fa[pre=x])//pre=x老漏。。
  57. 			Splay(x), rson=pre, Update(x);
  58. 	}
  59. 	void Make_root(int x){
  60. 		Access(x), Splay(x), Rev(x);
  61. 	}
  62. 	void Split(int x,int y){
  63. 		Make_root(x), Access(y), Splay(y);
  64. 	}
  65. 	int Find_root(int x)
  66. 	{
  67. 		Access(x), Splay(x);
  68. 		while(lson) x=lson;
  69. //		Splay(x);//需要保证复杂度？
  70. 		return x;
  71. 	}
  72. //	int Get_fa(int x){
  73. //		while(fa[x]) x=fa[x];
  74. //		return x;
  75. //	}
  76. //	void Link(int x,int y){
  77. //		if(Get_fa(x)!=Get_fa(y)) Make_root(x),fa[x]=y;//Get_fa()有点慢啊
  78. //	}
  79. 	void Link(int x,int y)
  80. 	{
  81. 		Make_root(x);
  82. 		if(Find_root(y)!=x) fa[x]=y;//连的是轻边！不要加rson=y;
  83. 	}
  84. 	void Cut(int x,int y)
  85. 	{
  86. 		Make_root(x);
  87. 		if(Find_root(y)==x&&fa[x]==y&&!rson)
  88. 			fa[x]=son[y][0]=0, Update(y);
  89. //		if(Find_root(y)==x&&fa[y]==x&&!son[y][0])//若Find_root()中把根又转回去了
  90. //			fa[y]=rson=0, Update(x);
  91. 	}
  92. }
  94. int main()
  95. {
  96. 	int n=read(),m=read();
  97. 	for(int i=1; i<=n; ++i) LCT::val[i]=read();
  98. 	int opt,x,y;
  99. 	while(m--)
  100. 	{
  101. 		opt=read(),x=read(),y=read();
  102. 		if(!opt) LCT::Split(x,y), printf("%d\n",LCT::sum[y]);
  103. 		else if(opt==1) LCT::Link(x,y);
  104. 		else if(opt==2) LCT::Cut(x,y);
  105. 		else LCT::Splay(x), LCT::val[x]=y;/*sum[y]^=val[x]^y, val[x]=y*/ //这步在Split()(查询时)的Splay()中 最后还会Update()
  106. 	}
  107. 	return 0;
  108. }

第二次(2018.4.5)：

  1. #include <cstdio>
  2. #include <cctype>
  3. #include <algorithm>
  4. #define gc() getchar()
  5. const int N=3e5+5;
  7. inline int read()
  8. {
  9. 	int now=0;register char c=gc();
  10. 	for(;!isdigit(c);c=gc());
  11. 	for(;isdigit(c);now=now*10+c-'0',c=gc());
  12. 	return now;
  13. }
  14. namespace LCT
  15. {
  16. 	#define lson son[x][0]
  17. 	#define rson son[x][1]
  19. 	int fa[N],son[N][2],sum[N],val[N],sk[N];
  20. 	bool rev[N];
  21. 	inline void Update(int x){
  22. 		sum[x]=sum[lson]^sum[rson]^val[x];
  23. 	}
  24. 	inline bool n_root(int x){
  25. 		return son[fa[x]][0]==x||son[fa[x]][1]==x;
  26. 	}
  27. 	inline void Rev(int x){
  28. 		std::swap(lson,rson), rev[x]^=1;
  29. 	}
  30. 	inline void PushDown(int x){
  31. 		if(rev[x]) Rev(lson),Rev(rson),rev[x]=0;
  32. 	}
  33. 	void Rotate(int x)
  34. 	{
  35. 		int a=fa[x],b=fa[a],l=son[a][1]==x,r=l^1;
  36. 		if(n_root(a)) son[b][son[b][1]==a]=x;
  37. 		if(son[x][r]) fa[son[x][r]]=a;
  38. 		fa[a]=x, fa[x]=b, son[a][l]=son[x][r], son[x][r]=a;
  39. 		Update(a);
  40. 	}
  41. 	void Splay(int x)
  42. 	{
  43. 		int t=1,a=x; sk[1]=x;
  44. 		while(n_root(a)) sk[++t]=a=fa[a];
  45. 		while(t) PushDown(sk[t--]);
  46. 		while(n_root(x))
  47. 		{
  48. 			if(n_root(a=fa[x])) Rotate(son[a][1]==x^son[fa[a]][1]==a?x:a);
  49. 			Rotate(x);
  50. 		}
  51. 		Update(x);
  52. 	}
  53. 	void Access(int x){
  54. 		for(int pre=0; x; x=fa[pre=x])
  55. 			Splay(x), rson=pre, Update(x);
  56. 	}
  57. 	void Make_root(int x){
  58. 		Access(x), Splay(x), Rev(x);
  59. 	}
  60. 	void Split(int x,int y){
  61. 		Make_root(x), Access(y), Splay(y);
  62. 	}
  63. 	int Find_root(int x)
  64. 	{
  65. 		Access(x), Splay(x);
  66. 		while(lson) x=lson;
  67. 		return x;
  68. 	}
  69. 	void Link(int x,int y)
  70. 	{
  71. 		Make_root(x);
  72. 		if(Find_root(y)!=x) fa[x]=y;
  73. 	}
  74. 	void Cut(int x,int y)
  75. 	{
  76. 		Make_root(x);
  77. 		if(Find_root(y)==x&&fa[x]==y&&!rson)
  78. 			fa[x]=son[y][0]=0, Update(y);
  79. 	}
  80. }
  81. using namespace LCT;
  83. int main()
  84. {
  85. 	int n=read(),m=read(),opt,x,y;
  86. 	for(int i=1; i<=n; ++i) val[i]=read();
  87. 	while(m--)
  88. 		switch(opt=read(),x=read(),y=read(),opt)
  89. 		{
  90. 			case 0: Split(x,y),printf("%d\n",sum[y]); break;
  91. 			case 1: Link(x,y); break;
  92. 			case 2: Cut(x,y); break;
  93. 			case 3: Splay(x), val[x]=y; break;
  94. 		}
  95. 	return 0;
  96. }

第三次（2018.6.28）：竟然快了那么多，好像没什么太大变化啊。。

  1. //592ms	3.76MB
  2. #include <cstdio>
  3. #include <cctype>
  4. #include <algorithm>
  5. //#define gc() getchar()
  6. #define MAXIN 150000
  7. #define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
  8. const int N=3e5+5;
  10. char IN[MAXIN],*SS=IN,*TT=IN;
  11. namespace LCT
  12. {
  13. 	#define lson son[x][0]
  14. 	#define rson son[x][1]
  15. 	int fa[N],son[N][2],val[N],sum[N],sk[N];
  16. 	bool rev[N];
  18. 	inline void Update(int x){
  19. 		sum[x]=sum[lson]^sum[rson]^val[x];
  20. 	}
  21. 	inline void Rev(int x){
  22. 		std::swap(lson,rson), rev[x]^=1;
  23. 	}
  24. 	inline void PushDown(int x){
  25. 		if(rev[x]) Rev(lson), Rev(rson), rev[x]=0;
  26. 	}
  27. 	inline bool n_root(int x){
  28. 		return son[fa[x]][0]==x||son[fa[x]][1]==x;
  29. 	}
  30. 	void Rotate(int x)
  31. 	{
  32. 		int a=fa[x],b=fa[a],l=son[a][1]==x,r=l^1;
  33. 		if(n_root(a)) son[b][son[b][1]==a]=x;
  34. 		if(son[x][r]) fa[son[x][r]]=a;
  35. 		fa[x]=b, fa[a]=x, son[a][l]=son[x][r], son[x][r]=a;
  36. 		Update(a);
  37. 	}
  38. 	void Splay(int x)
  39. 	{
  40. 		int a=x,t=1; sk[1]=x;
  41. 		while(n_root(a)) sk[++t]=a=fa[a];
  42. 		while(t) PushDown(sk[t--]);
  43. 		while(n_root(x))
  44. 		{
  45. 			if(n_root(a=fa[x])) Rotate(son[a][1]==x^son[fa[a]][1]==a?x:a);
  46. 			Rotate(x);
  47. 		}
  48. 		Update(x);
  49. 	}
  50. 	void Access(int x){
  51. 		for(int pre=0; x; x=fa[pre=x])
  52. 			Splay(x), rson=pre, Update(x);
  53. 	}
  54. 	void Make_root(int x){
  55. 		Access(x), Splay(x), Rev(x);
  56. 	}
  57. 	void Split(int x,int y){
  58. 		Make_root(x), Access(y), Splay(y);
  59. 	}
  60. 	int Find_root(int x)
  61. 	{
  62. 		Access(x), Splay(x);
  63. 		while(lson) x=lson;
  64. 		return x;
  65. 	}
  66. 	void Link(int x,int y)
  67. 	{
  68. 		Make_root(x);
  69. 		if(Find_root(y)!=x) fa[x]=y;
  70. 	}
  71. 	void Cut(int x,int y)
  72. 	{
  73. 		Make_root(x);
  74. 		if(Find_root(y)==x&&!rson&&fa[x]==y) son[y][0]=fa[x]=0, Update(y);
  75. 	}
  76. }
  77. using namespace LCT;
  79. inline int read()
  80. {
  81. 	int now=0;register char c=gc();
  82. 	for(;!isdigit(c);c=gc());
  83. 	for(;isdigit(c);now=now*10+c-'0',c=gc());
  84. 	return now;
  85. }
  87. int main()
  88. {
  89. 	int n=read(), m=read();
  90. 	for(int i=1; i<=n; ++i) val[i]=read();
  91. 	int opt,x,y;
  92. 	while(m--){
  93. 		switch(opt=read(),x=read(),y=read(),opt){
  94. 			case 0: Split(x,y), printf("%d\n",sum[y]); break;
  95. 			case 1: Link(x,y); break;
  96. 			case 2: Cut(x,y); break;
  97. 			case 3: Splay(x), val[x]=y; break;
  98. 		}
  99. 	}
  100. 	return 0;
  101. }

19.4.5

  1. //520ms	2984KB
  2. #include <cstdio>
  3. #include <cctype>
  4. #include <algorithm>
  5. #define gc() getchar()
  6. #define MAXIN 300000
  7. //#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
  8. typedef long long LL;
  9. const int N=3e5+5;
  11. char IN[MAXIN],*SS=IN,*TT=IN;
  12. struct LCT
  13. {
  14. 	#define ls son[x][0]
  15. 	#define rs son[x][1]
  16. 	int fa[N],son[N][2],sum[N],val[N],sk[N];
  17. 	bool rev[N];
  18. 	inline void Update(int x)
  19. 	{
  20. 		sum[x]=sum[ls]^sum[rs]^val[x];
  21. 	}
  22. 	inline bool n_root(int x)
  23. 	{
  24. 		return son[fa[x]][0]==x||son[fa[x]][1]==x;
  25. 	}
  26. 	inline void Rev(int x)
  27. 	{
  28. 		std::swap(ls,rs), rev[x]^=1;
  29. 	}
  30. 	inline void PushDown(int x)
  31. 	{
  32. 		if(rev[x]) Rev(ls), Rev(rs), rev[x]=0;
  33. 	}
  34. 	void Rotate(int x)
  35. 	{
  36. 		int a=fa[x],b=fa[a],l=son[a][1]==x,r=l^1;
  37. 		if(n_root(a)) son[b][son[b][1]==a]=x;
  38. 		if(son[x][r]) fa[son[x][r]]=a;
  39. 		fa[x]=b, fa[a]=x, son[a][l]=son[x][r], son[x][r]=a;
  40. 		Update(a);
  41. 	}
  42. 	void Splay(int x)
  43. 	{
  44. 		int t=1,a=x; sk[1]=a;
  45. 		while(n_root(a)) sk[++t]=a=fa[a];
  46. 		while(t) PushDown(sk[t--]);
  47. 		while(n_root(x))
  48. 		{
  49. 			if(n_root(a=fa[x])) Rotate(son[a][0]==x^son[fa[a]][0]==a?x:a);
  50. 			Rotate(x);
  51. 		}
  52. 		Update(x);
  53. 	}
  54. 	void Access(int x)
  55. 	{
  56. 		for(int pre=0; x; x=fa[pre=x])
  57. 			Splay(x), rs=pre, Update(x);
  58. 	}
  59. 	void MakeRoot(int x)
  60. 	{
  61. 		Access(x), Splay(x), Rev(x);
  62. 	}
  63. 	void Split(int x,int y)
  64. 	{
  65. 		MakeRoot(x), Access(y), Splay(y);
  66. 	}
  67. 	int FindRoot(int x)
  68. 	{
  69. 		Access(x), Splay(x);
  70. 		while(ls) x=ls;
  71. 		return x;
  72. 	}
  73. 	void Link(int x,int y)
  74. 	{
  75. 		MakeRoot(x);
  76. 		if(FindRoot(y)!=x) fa[x]=y;
  77. 	}
  78. 	void Cut(int x,int y)
  79. 	{
  80. 		MakeRoot(x);
  81. 		if(FindRoot(y)==x&&fa[x]==y&&!rs) fa[x]=son[y][0]=0, Update(y);
  82. 	}
  83. }T;
  85. inline int read()
  86. {
  87. 	int now=0;register char c=gc();
  88. 	for(;!isdigit(c);c=gc());
  89. 	for(;isdigit(c);now=now*10+c-48,c=gc());
  90. 	return now;
  91. }
  93. int main()
  94. {
  95. 	int n=read(),q=read();
  96. 	for(int i=1; i<=n; ++i) T.val[i]=read();
  97. 	for(int x,y; q--; )
  98. 		switch(read())
  99. 		{
  100. 			case 0: x=read(),y=read(),T.Split(x,y),printf("%d\n",T.sum[y]); break;
  101. 			case 1: x=read(),y=read(),T.Link(x,y); break;
  102. 			case 2: x=read(),y=read(),T.Cut(x,y); break;
  103. 			case 3: T.Splay(x=read()),T.val[x]=read(); break;
  104. 		}
  106. 	return 0;
  107. }

洛谷.3690.[模板]Link Cut Tree(动态树)的更多相关文章

  1. LCT总结——概念篇+洛谷P3690[模板]Link Cut Tree(动态树)（LCT，Splay）

    为了优化体验(其实是强迫症),蒟蒻把总结拆成了两篇,方便不同学习阶段的Dalao们切换. LCT总结--应用篇戳这里 概念.性质简述 首先介绍一下链剖分的概念(感谢laofu的讲课) 链剖分,是指一类 ...

  2. 洛谷P3690 [模板] Link Cut Tree [LCT]

    题目传送门 Link Cut Tree 题目背景 动态树 题目描述 给定n个点以及每个点的权值,要你处理接下来的m个操作.操作有4种.操作从0到3编号.点从1到n编号. 0:后接两个整数(x,y),代 ...

  3. Link Cut Tree 动态树 小结

    动态树有些类似 树链剖分+并查集 的思想,是用splay维护的 lct的根是动态的,"轻重链"也是动态的,所以并没有真正的轻重链 动态树的操作核心是把你要把 修改/询问/... 等 ...

  4. 洛谷P3690 Link Cut Tree (动态树)

    干脆整个LCT模板吧. 缺个链上修改和子树操作,链上修改的话join(u,v)然后把v splay到树根再打个标记就好. 至于子树操作...以后有空的话再学(咕咕咕警告) #include<bi ...

  5. LCT(link cut tree) 动态树

    模板参考:https://blog.csdn.net/saramanda/article/details/55253627 综合各位大大博客后整理的模板: #include<iostream&g ...

  6. 洛谷P2633 Count on a tree(主席树上树)

    题目描述 给定一棵N个节点的树,每个点有一个权值,对于M个询问(u,v,k),你需要回答u xor lastans和v这两个节点间第K小的点权.其中lastans是上一个询问的答案,初始为0,即第一个 ...

  7. 模板Link Cut Tree （动态树）

    题目描述 给定N个点以及每个点的权值,要你处理接下来的M个操作.操作有4种.操作从0到3编号.点从1到N编号. 0:后接两个整数(x,y),代表询问从x到y的路径上的点的权值的xor和.保证x到y是联 ...

  8. 洛谷 P2633 Count on a tree 主席树

    在一棵树上,我们要求点 $(u,v)$ 之间路径的第$k$大数. 对于点 $i$  ,建立 $i$  到根节点的一棵前缀主席树. 简单容斥后不难得出结果为$sumv[u]+sumv[v]−sumv[l ...

  9. 洛谷P2633 Count on a tree 主席树

    传送门:主席树 解题报告: 传送门! umm这题我还麻油开始做 所以 先瞎扯一波我的想法,如果错了我就当反面教材解释这种典型错误,对了我就不管了QwQ 就直接dfs,在dfs的过程中建树 然后就直接查 ...

随机推荐

  1. python3之Django多数据库

    1.定义数据库 在django项目中, 一个工程中存在多个APP应用很常见:有时候希望不同的APP连接不同的数据库,这个时候需要建立多个数据库连接.在Django的setting中使用DATABASE ...

  2. 【转】OpenCV对图片中的RotatedRect进行填充

    函数名:full_rotated_rect 函数参数: image输入图像,rect希望在图像中填充的RotatedRect,color填充的颜色 主要的思路是:先找到RotatedRect的四个顶点 ...

  3. 028_rync和inotify实现实时备份

    一.服务节点安装inotify-tools. 确保系统后以下输出=> [root@xxxx]# ll /proc/sys/fs/inotify/ total 0 -rw-r--r-- 1 roo ...

  4. nginx报错：403 Forbidden 并且访问首页index.php是下载文件的状态

    nginx报错:403 Forbidden 并且访问首页index.php是下载文件的状态,不能正常解析php 系统有其他两个站访问是正常的 看日志没有看到明显的错误 搜索了下: 答案如下: php的 ...

  5. Centos6安装FreeSWITCH 1.5时./configure问题解决记录

    系统:Centos 6.4 64位: FreeSWITCH版本:1.5 具体的安装过程参考FreeSWITCH 官网wiki (也可以参考我的博客<Centos6安装FreeSWITCH> ...

  6. 初识numpy

    from numpy import *   导入numpy包 random可以生成随机数组 通过mat函数,将数组转换成矩阵,可以对矩阵进行求逆计算等.其中.I操作实现了矩阵求逆计算操作. 执行矩阵乘 ...

  7. JAVA Random 随机类

    nextInt 方法 得到一个随机整数, 可以指定范围 package object; import static net.util.Print.*; import java.util.Random; ...

  8. ECLIPSE 导入外部文件或源码包

    步骤: 点击Project->Properties->Libraries->Add External Class Folder.. ->选择你的文件路径->确定 注:如果 ...

  9. Windows安装使用Openssl

    1.什么是openssl? 2.下载安装 三方下载地址 备用64位和32位下载地址 选择32位或者64位合适的版本下载,例如Win64OpenSSL_Light-1_0_2h.exe: 设置环境变量, ...

  10. 学习Struts2经验总结

    一.struts 访问路径问题 1) Struts2的思想:主要围着“action”转,只要找到“action”它就知道自己该干嘛了. 首先配置struts.xml ,我们可以明白的看到,action ...