题目链接

无优化版本(170行)：

  1. /*
  2. 首先树剖可以维护树上的链Sum、Max
  3. 可以对每个宗教建一棵线段树，那这题就很好做了
  4. 不过10^5需要动态开点 (不明白为什么nlogn不需要回收就可以 不是每个Insert加log个节点?)
  5. 操作修改完更改原数列！盲人。。
  6. */
  7. #include<cstdio>
  8. #include<cctype>
  9. #include<algorithm>
  10. #define gc() getchar()
  11. #define now node[rt]
  12. #define lson node[rt].ls
  13. #define rson node[rt].rs
  14. const int N=1e5+5,MAXN=N*19;
  16. int n,q,cnt,H[N],Enum,to[N<<1],nxt[N<<1],fa[N],dep[N],id[N],son[N],sz[N],top[N];
  17. int idx,pos[N],W[N],rel[N],root[N],pool[MAXN];
  18. struct Seg_Tree
  19. {
  20. 	struct Node
  21. 	{
  22. 		int sum,maxn,ls,rs,l,r;
  23. 		inline void Init() {sum=maxn=l=r=ls=rs=0;}
  24. 	}node[MAXN];
  25. 	inline int new_Node() {return pool[++idx];}
  26. 	inline void del_Node(int rt) {now.Init(), pool[idx--]=rt;}
  27. 	inline void PushUp(int rt)
  28. 	{
  29. 		now.sum = node[lson].sum + node[rson].sum,
  30. 		now.maxn= std::max(node[lson].maxn, node[rson].maxn);
  31. 	}
  32. 	void Insert(int l,int r,int &rt,int x)
  33. 	{
  34. 		if(!rt) rt=new_Node(), now.l=l, now.r=r;
  35. 		if(l==r) now.sum=now.maxn=W[x];
  36. 		else
  37. 		{
  38. 			int m=l+r>>1;
  39. 			if(id[x]<=m) /*lson?0:lson=new_Node(),*/ Insert(l,m,lson,x);//不想打括号--
  40. 			else /*rson?0:rson=new_Node(),*/ Insert(m+1,r,rson,x);
  41. 			PushUp(rt);
  42. 		}
  43. 	}
  44. 	void Delete(int rt,int x,int &son)//可以用Insert将其值变为0来删除 但是就不能重复利用了
  45. 	{
  46. 		if(now.l==now.r) son=0, del_Node(rt);//父亲的这个儿子也要删!
  47. 		else
  48. 		{
  49. 			int m=now.l+now.r>>1;
  50. 			if(id[x]<=m) Delete(lson,x,lson);
  51. 			else Delete(rson,x,rson);
  52. 			if(now.sum==W[x]) del_Node(rt), son=0;
  53. 			else PushUp(rt);
  54. 		}
  55. 	}
  56. 	void Modify(int rt,int p,int v)
  57. 	{
  58. 		if(now.l==now.r) now.maxn=now.sum=v;
  59. 		else
  60. 		{
  61. 			int m=now.l+now.r>>1;
  62. 			if(p<=m) Modify(lson,p,v);
  63. 			else Modify(rson,p,v);
  64. 			PushUp(rt);
  65. 		}
  66. 	}
  67. 	int Query_Max(int rt,int L,int R)
  68. 	{
  69. 		if(!rt) return 0;
  70. 		if(L<=now.l && now.r<=R) return now.maxn;
  71. 		int m=now.l+now.r>>1;
  72. 		if(L<=m)
  73. 			if(m<R) return std::max(Query_Max(lson,L,R),Query_Max(rson,L,R));
  74. 			else return Query_Max(lson,L,R);
  75. 		return Query_Max(rson,L,R);
  76. 	}
  77. 	int Query_Sum(int rt,int L,int R)
  78. 	{
  79. 		if(!rt) return 0;
  80. 		if(L<=now.l && now.r<=R) return now.sum;
  81. 		int m=now.l+now.r>>1;
  82. 		if(L<=m)
  83. 			if(m<R) return Query_Sum(lson,L,R)+Query_Sum(rson,L,R);
  84. 			else return Query_Sum(lson,L,R);
  85. 		return Query_Sum(rson,L,R);
  86. 	}
  87. }t;
  88. #undef now
  89. inline int read()
  90. {
  91. 	int now=0,f=1;register char c=gc();
  92. 	for(;!isdigit(c);c=gc()) if(c=='-') f=-1;
  93. 	for(;isdigit(c);now=now*10+c-'0',c=gc());
  94. 	return now*f;
  95. }
  96. inline void AddEdge(int u,int v)
  97. {
  98. 	to[++Enum]=v, nxt[Enum]=H[u], H[u]=Enum;
  99. 	to[++Enum]=u, nxt[Enum]=H[v], H[v]=Enum;
  100. }
  101. void DFS1(int x)
  102. {
  103. 	sz[x]=1; int mx=0;
  104. 	for(int v,i=H[x]; i; i=nxt[i])
  105. 		if((v=to[i])!=fa[x])
  106. 		{
  107. 			fa[v]=x, dep[v]=dep[x]+1, DFS1(v), sz[x]+=sz[v];
  108. 			if(mx<sz[v]) mx=sz[v], son[x]=v;
  109. 		}
  110. }
  111. void DFS2(int x,int tp)
  112. {
  113. 	id[x]=++cnt, top[x]=tp;
  114. 	if(son[x])
  115. 	{
  116. 		DFS2(son[x],tp);
  117. 		for(int i=H[x]; i; i=nxt[i])
  118. 			if(to[i]!=fa[x] && to[i]!=son[x])
  119. 				DFS2(to[i], to[i]);
  120. 	}
  121. }
  122. void Query_Chain_Sum(int x,int y)
  123. {
  124. 	long long res=0;int rt=root[rel[x]];
  125. 	while(top[x]!=top[y])
  126. 	{
  127. 		if(dep[top[x]]<dep[top[y]]) std::swap(x,y);
  128. 		res += t.Query_Sum(rt,id[top[x]],id[x]);
  129. 		x=fa[top[x]];
  130. 	}
  131. 	if(dep[x]>dep[y]) std::swap(x,y);
  132. 	res += t.Query_Sum(rt,id[x],id[y]);
  133. 	printf("%lld\n",res);
  134. }
  135. void Query_Chain_Max(int x,int y)
  136. {
  137. 	int res=0, rt=root[rel[x]];
  138. 	while(top[x]!=top[y])
  139. 	{
  140. 		if(dep[top[x]]<dep[top[y]]) std::swap(x,y);
  141. 		res = std::max(res, t.Query_Max(rt,id[top[x]],id[x]));
  142. 		x=fa[top[x]];
  143. 	}
  144. 	if(dep[x]>dep[y]) std::swap(x,y);
  145. 	res = std::max(res, t.Query_Max(rt,id[x],id[y]));
  146. 	printf("%d\n",res);
  147. }
  149. int main()
  150. {
  151. #ifndef ONLINE_JUDGE
  152. 	freopen("3531.in","r",stdin);
  153. 	freopen("my.out","w",stdout);
  154. #endif
  156. 	n=read(),q=read();
  157. 	for(int i=1; i<MAXN; ++i) pool[i]=i;
  158. 	for(int i=1; i<=n; ++i) W[i]=read(), rel[i]=read();
  159. 	for(int u,v,i=1; i<n; ++i) u=read(),v=read(),AddEdge(u,v);
  160. 	DFS1(1), DFS2(1,1);
  161. 	for(int i=1; i<=n; ++i) t.Insert(1,n,root[rel[i]],i);
  162. 	char opt[5]; int x,y;
  163. 	while(q--)
  164. 	{
  165. 		scanf("%s",opt), x=read(), y=read();
  166. 		switch(opt[1])
  167. 		{
  168. 			case 'C': t.Delete(root[rel[x]],x,root[rel[x]]), rel[x]=y, t.Insert(1,n,root[y],x);
  169. 					  break;
  170. 			case 'W': t.Modify(root[rel[x]],id[x],y), W[x]=y;//!
  171. 					  break;
  172. 			case 'S': Query_Chain_Sum(x,y);
  173. 					  break;
  174. 			case 'M': Query_Chain_Max(x,y);
  175. 					  break;
  176. 		}
  177. 	}
  178. 	return 0;
  179. }

优化(短多了)：

  1. //不知道为什么只需nlogn个节点 不需要回收
  2. //+各种函数简化
  3. #include<cstdio>
  4. #include<cctype>
  5. #include<algorithm>
  6. #define gc() getchar()
  7. #define now node[rt]
  8. #define lson node[rt].ls
  9. #define rson node[rt].rs
  10. const int N=1e5+5,MAXN=N*19;
  12. int n,q,cnt,H[N],Enum,to[N<<1],nxt[N<<1],fa[N],dep[N],id[N],son[N],sz[N],top[N];
  13. int idx,pos[N],W[N],rel[N],root[N],pool[MAXN];
  14. struct Seg_Tree
  15. {
  16. 	struct Node
  17. 	{
  18. 		int sum,maxn,ls,rs;
  19. 		inline void Init() {sum=maxn=ls=rs=0;}
  20. 	}node[MAXN];
  21. 	inline int new_Node() {return ++idx;}
  22. //	inline int new_Node() {return pool[++idx];}
  23. //	inline void del_Node(int rt) {now.Init(), pool[idx--]=rt;}
  24. 	inline void PushUp(int rt)
  25. 	{
  26. 		now.sum = node[lson].sum + node[rson].sum,
  27. 		now.maxn= std::max(node[lson].maxn, node[rson].maxn);
  28. 	}
  29. 	void Update(int l,int r,int &rt,int p,int v)
  30. 	{
  31. 		if(!rt) rt=new_Node();
  32. 		if(l==r) now.sum=now.maxn=v;
  33. 		else
  34. 		{
  35. 			int m=l+r>>1;
  36. 			if(p<=m) Update(l,m,lson,p,v);
  37. 			else Update(m+1,r,rson,p,v);
  38. 			PushUp(rt);
  39. 		}
  40. 	}
  41. 	int Query(int l,int r,int rt,int L,int R,bool opt)
  42. 	{
  43. 		if(!rt) return 0;
  44. 		if(L<=l && r<=R) return opt?now.sum:now.maxn;
  45. 		int m=l+r>>1;
  46. 		if(L<=m)
  47. 			if(m<R) return opt?Query(l,m,lson,L,R,opt)+Query(m+1,r,rson,L,R,opt)
  48. 							  :std::max(Query(l,m,lson,L,R,opt),Query(m+1,r,rson,L,R,opt));
  49. 			else return Query(l,m,lson,L,R,opt);
  50. 		return Query(m+1,r,rson,L,R,opt);
  51. 	}
  52. }t;
  53. #undef now
  54. inline int read()
  55. {
  56. 	int now=0,f=1;register char c=gc();
  57. 	for(;!isdigit(c);c=gc()) if(c=='-') f=-1;
  58. 	for(;isdigit(c);now=now*10+c-'0',c=gc());
  59. 	return now*f;
  60. }
  61. inline void AddEdge(int u,int v)
  62. {
  63. 	to[++Enum]=v, nxt[Enum]=H[u], H[u]=Enum;
  64. 	to[++Enum]=u, nxt[Enum]=H[v], H[v]=Enum;
  65. }
  66. void DFS1(int x)
  67. {
  68. 	sz[x]=1; int mx=-1;
  69. 	for(int v,i=H[x]; i; i=nxt[i])
  70. 		if((v=to[i])!=fa[x])
  71. 		{
  72. 			fa[v]=x, dep[v]=dep[x]+1, DFS1(v), sz[x]+=sz[v];
  73. 			if(mx<sz[v]) mx=sz[v], son[x]=v;
  74. 		}
  75. }
  76. void DFS2(int x,int tp)
  77. {
  78. 	id[x]=++cnt, top[x]=tp;
  79. 	if(son[x])
  80. 	{
  81. 		DFS2(son[x],tp);
  82. 		for(int i=H[x]; i; i=nxt[i])
  83. 			if(to[i]!=fa[x] && to[i]!=son[x])
  84. 				DFS2(to[i], to[i]);
  85. 	}
  86. }
  87. void Query_Chain(int x,int y,bool opt)
  88. {
  89. 	int res=0;int rt=root[rel[x]];
  90. 	while(top[x]!=top[y])
  91. 	{
  92. 		if(dep[top[x]]<dep[top[y]]) std::swap(x,y);
  93. 		if(opt) res += t.Query(1,n,rt,id[top[x]],id[x],1);
  94. 		else res = std::max(res, t.Query(1,n,rt,id[top[x]],id[x],0));
  95. 		x=fa[top[x]];
  96. 	}
  97. 	if(id[x]>id[y]) std::swap(x,y);
  98. 	if(opt) res += t.Query(1,n,rt,id[x],id[y],1);
  99. 	else res = std::max(res, t.Query(1,n,rt,id[x],id[y],0));
  100. 	printf("%d\n",res);
  101. }
  103. int main()
  104. {
  105. #ifndef ONLINE_JUDGE
  106. 	freopen("35312.in","r",stdin);
  107. 	freopen("my.out","w",stdout);
  108. #endif
  110. 	n=read(),q=read();
  111. 	for(int i=1; i<MAXN; ++i) pool[i]=i;
  112. 	for(int i=1; i<=n; ++i) W[i]=read(), rel[i]=read();
  113. 	for(int u,v,i=1; i<n; ++i) u=read(),v=read(),AddEdge(u,v);
  114. 	DFS1(1), DFS2(1,1);
  115. 	for(int i=1; i<=n; ++i) t.Update(1,n,root[rel[i]],id[i],W[i]);
  116. 	char opt[5]; int x,y;
  117. 	while(q--)
  118. 	{
  119. 		scanf("%s",opt), x=read(), y=read();
  120. 		switch(opt[1])
  121. 		{
  122. 			case 'C': t.Update(1,n,root[rel[x]],id[x],0), rel[x]=y, t.Update(1,n,root[y],id[x],W[x]);
  123. 					  break;
  124. 			case 'W': t.Update(1,n,root[rel[x]],id[x],y),W[x]=y/*!*/; break;
  125. 			case 'S': Query_Chain(x,y,1); break;
  126. 			case 'M': Query_Chain(x,y,0); break;
  127. 		}
  128. 	}
  129. 	return 0;
  130. }

BZOJ.3531.旅行(树链剖分 动态开点)的更多相关文章

  1. BZOJ 3531 [Sdoi2014]旅行 树链剖分+动态开点线段树

    题意 S国有N个城市,编号从1到N.城市间用N-1条双向道路连接,满足从一个城市出发可以到达其它所有城市.每个城市信仰不同的宗教,如飞天面条神教.隐形独角兽教.绝地教都是常见的信仰. 为了方便,我们用 ...

  2. 洛谷P3313 [SDOI2014]旅行(树链剖分 动态开节点线段树)

    题意 题目链接 Sol 树链剖分板子 + 动态开节点线段树板子 #include<bits/stdc++.h> #define Pair pair<int, int> #def ...

  3. 【BZOJ3531】[Sdoi2014]旅行 树链剖分+动态开点线段树

    [BZOJ3531][Sdoi2014]旅行 Description S国有N个城市,编号从1到N.城市间用N-1条双向道路连接,满足从一个城市出发可以到达其它所有城市.每个城市信仰不同的宗教,如飞天 ...

  4. 【bzoj4999】This Problem Is Too Simple！ 树链剖分+动态开点线段树

    题目描述 给您一颗树,每个节点有个初始值. 现在支持以下两种操作: 1. C i x(0<=x<2^31) 表示将i节点的值改为x. 2. Q i j x(0<=x<2^31) ...

  5. [bzoj 3531][SDOI2014]旅行（树链剖分+动态开点线段树）

    题目:http://www.lydsy.com:808/JudgeOnline/problem.php?id=3531 分析: 对于每个颜色(颜色<=10^5)都建立一颗线段树 什么!那么不是M ...

  6. 2018.07.24 bzoj3531: [Sdoi2014]旅行（树链剖分+动态开点）

    传送门 树链剖分. 如何维护? 如果颜色少直接每种颜色一颗线段树走人. 但这题颜色数量不大于1e5" role="presentation" style="po ...

  7. BZOJ 2157 旅行(树链剖分码农题)

    写了5KB,1发AC... 题意:给出一颗树,支持5种操作. 1.修改某条边的权值.2.将u到v的经过的边的权值取负.3.求u到v的经过的边的权值总和.4.求u到v的经过的边的权值最大值.5.求u到v ...

  8. bzoj3531: [Sdoi2014]旅行 （树链剖分 && 动态开点线段树）

    感觉动态开点线段树空间复杂度好优秀呀 树剖裸题 把每个宗教都开一颗线段树就可以了 但是我一直TLE 然后调了一个小时 为什么呢 因为我 #define max(x, y) (x > y ? x ...

  9. BZOJ 3531（树链剖分+线段树）

    Problem 旅行 (BZOJ 3531) 题目大意 给定一颗树,树上的每个点有两个权值(x,y). 要求维护4种操作: 操作1:更改某个点的权值x. 操作2:更改某个点的权值y. 操作3:求a-- ...

随机推荐

  1. 【vim】缩写 :ab [缩写] [要替换的文字]

    一个很可能是最令人印象深刻的窍门是你可以在 Vim 中定义缩写,它可以实时地把你输入的东西替换为另外的东西.语法格式如下: :ab [缩写] [要替换的文字] 一个通用的例子是: :ab asap a ...

  2. ARMV8 datasheet学习笔记3：AArch64应用级体系结构之Memory Type and Attributes

    1.前言 2. Memory类型和属性 memory分为normal memory和device memory,两种类型的Memory有各自的属性,除了下面介绍的几种属性外,还有其他一些杂项属性 2. ...

  3. 百度地图的Icon

    在百度地图的类说明中,查看对Icon的构建: 定制IconOptions 看下面的差别 发现在IconOptions没有imageSize属性 而在实际测试中,代码如下 <script type ...

  4. RestTemplate的使用

    1.postForObject :传入一个业务对象,返回是一个String 调用方: BaseUser baseUser=new BaseUser(); baseUser.setUserid(user ...

  5. python读取两个csv文件数据，进行查找匹配出现次数

    现有需求 表1 表2 需要拿表1中的编码去表2中的门票编码列匹配,统计出现的次数,由于表2编码列是区域间,而且列不是固定的,代码如下 #encoding:utf-8 ##导入两个CSV进行比对 imp ...

  6. python 全栈开发，Day64(视图,触发器,函数,存储过程,事务)

    昨日内容回顾 pymysql:属于python的一个模块 pip3 install pymysql conn = pymysql.connect(...,charset = 'uft8') 创建游标 ...

  7. 步步为营-71-asp.net的简单练习(图片处理)

    1 原有图片添加水印 1.1 封装一个类,用于获取文件路径 using System; using System.Collections.Generic; using System.IO; using ...

  8. JQuery编写自己的插件（七）

    一:jQuery插件的编写基础1.插件的种类编写插件的目的是给一系列已经方法或函数做一个封装,以便在其他地方重复使用,方便后期维护和提高开发效率.常见的种类有以下三种:封装对象方法的插件

  9. hdu 1728 迷宫 给定最大转弯次数 (BFS)

    给出起点 终点 以及转弯次数 在<=转弯次数的条件 能否走到终点 Sample Input25 5...** // .可走 *不可走*.**...........*....1 1 1 1 3 / ...

  10. zjoi2017 仙人掌

    题解: 好难的dp啊...看题解看了好久才看懂 http://blog.csdn.net/akak__ii/article/details/65935711 如果一开始的图就不是仙人掌,答案显然为0, ...