"BestCoder"杯中国大学生程序设计冠军赛 HDU 5221 Occupation

举报
Linux猿 发表于 2021/08/04 23:34:25 2021/08/04
【摘要】 题目链接~~> 做题感悟 :区域赛过后就没写过树剖 ,只记得思想,比赛时想到了,但是只打代码就打了半个多小时(真是醉了!!),然后就是不断的调试代码,悲催的是调了一个多小时也没调出来。。。。。 解题思路:这题只要想到某个节点的子树的所有节点编号都大于此节点的编号(在线段树中的位置的编号,树剖的时候编号都是连续的)且是连续的,那么我们只要深搜的时候记录一下子树最大的一...

题目链接~~>

做题感悟 :区域赛过后就没写过树剖 ,只记得思想,比赛时想到了,但是只打代码就打了半个多小时(真是醉了!!),然后就是不断的调试代码,悲催的是调了一个多小时也没调出来。。。。。

解题思路:这题只要想到某个节点的子树的所有节点编号都大于此节点的编号(在线段树中的位置的编号,树剖的时候编号都是连续的)且是连续的,那么我们只要深搜的时候记录一下子树最大的一个子节点的时间戳(在线段树中的编号),只要我们利用一个节点的进去的时间戳和出去的时间戳就可以对这个节点的子树操作了(CF上有一个处理子树类似的题),其它的两个操作就 so easy 了,线段树写的很挫,请见谅。。。

代码:


  
  1. #pragma comment(linker, "/STACK:102400000,102400000")
  2. #include <iostream>
  3. #include <string.h>
  4. #include <math.h>
  5. #include <queue>
  6. #include <algorithm>
  7. #include <stdlib.h>
  8. #include <map>
  9. #include <set>
  10. #include <stdio.h>
  11. #include <vector>
  12. #define L(x) (x*2)
  13. #define R(x) (x*2 + 1)
  14. using namespace std;
  15. const int INF = 0x3f3f3f3f ;
  16. const int MX = 100000 + 5 ;
  17. int n ,m ,num ,idx ,sum ;
  18. int dep[MX] ,top[MX] ,V[MX] ,ti[MX] ,to[MX] ,head[MX] ,father[MX] ,siz[MX] ,son[MX] ;
  19. struct Edge
  20. {
  21. int v ,next ;
  22. }E[MX*2] ;
  23. void add_edge(int u ,int v)
  24. {
  25. E[num].v = v ; E[num].next = head[u] ; head[u] = num++ ;
  26. E[num].v = u ; E[num].next = head[v] ; head[v] = num++ ;
  27. }
  28. void dfs_find(int u ,int fa)
  29. {
  30. dep[u] = dep[fa] + 1 ;
  31. father[u] = fa ;
  32. siz[u] = 1 ;
  33. son[u] = 0 ;
  34. for(int i = head[u] ;i != -1 ; i = E[i].next)
  35. {
  36. int v = E[i].v ;
  37. if(v == fa) continue ;
  38. dfs_find(v ,u) ;
  39. siz[u] += siz[v] ;
  40. if(siz[son[u]] < siz[v]) son[u] = v ;
  41. }
  42. }
  43. void dfs_time(int u ,int fa)
  44. {
  45. ti[u] = idx++ ;
  46. top[u] = fa ;
  47. if(son[u]) dfs_time(son[u] ,top[u]) ;
  48. for(int i = head[u] ;i != -1 ;i = E[i].next)
  49. {
  50. int v = E[i].v ;
  51. if(v == father[u] || v == son[u]) continue ;
  52. dfs_time(v ,v) ;
  53. }
  54. to[u] = idx-1 ;
  55. }
  56. void init()
  57. {
  58. sum = 0 ;
  59. num = 0 ;
  60. memset(head ,-1 ,sizeof(head)) ;
  61. }
  62. //线段树 ---------------------
  63. struct node
  64. {
  65. int le ,rt ,sum ,flag ;
  66. }T[MX*4] ;
  67. void build(int x ,int le ,int rt)
  68. {
  69. T[x].sum = 0 ;
  70. T[x].flag = true ;
  71. T[x].le = le ; T[x].rt = rt ;
  72. if(le == rt)
  73. return ;
  74. int Mid = (le + rt)>>1 ;
  75. build(L(x) ,le ,Mid) ;
  76. build(R(x),Mid+1 ,rt) ;
  77. }
  78. void push_up(int x)
  79. {
  80. T[x].sum = T[L(x)].sum + T[R(x)].sum ;
  81. T[x].flag = true ;
  82. }
  83. void push_down(int x)
  84. {
  85. if(!T[x].flag)
  86. {
  87. T[L(x)].sum = T[R(x)].sum = 0 ;
  88. T[L(x)].flag = T[R(x)].flag = false ;
  89. T[x].sum = 0 ;
  90. }
  91. }
  92. bool Findpos(int x ,int pos)
  93. {
  94. if(T[x].le == T[x].rt && T[x].le == pos)
  95. {
  96. if(!T[x].flag) T[x].sum = 0 ;
  97. return T[x].flag ;
  98. }
  99. push_down(x) ;
  100. int Mid = (T[x].le + T[x].rt)>>1 ;
  101. if(pos <= Mid) return Findpos(L(x) ,pos) ;
  102. else return Findpos(R(x) ,pos) ;
  103. push_up(x) ;
  104. }
  105. void insert(int x ,int pos ,int w) // 修改点
  106. {
  107. if(T[x].le == T[x].rt && T[x].le == pos)
  108. {
  109. T[x].sum = w ;
  110. T[x].flag = true ;
  111. return ;
  112. }
  113. push_down(x) ;
  114. int Mid = (T[x].le + T[x].rt)>>1 ;
  115. if(pos <= Mid) insert(L(x) ,pos ,w) ;
  116. else insert(R(x) ,pos ,w) ;
  117. push_up(x) ;
  118. }
  119. void update(int x ,int Le ,int Rt) // 区间清零
  120. {
  121. if(T[x].le == Le && T[x].rt == Rt)
  122. {
  123. T[x].flag = false ;
  124. T[x].sum = 0 ;
  125. return ;
  126. }
  127. push_down(x) ;
  128. int Mid = (T[x].le + T[x].rt)>>1 ;
  129. if(Le > Mid) update(R(x) ,Le ,Rt) ;
  130. else if(Rt <= Mid) update(L(x) ,Le ,Rt) ;
  131. else
  132. {
  133. update(L(x) ,Le ,Mid) ;
  134. update(R(x) ,Mid+1 ,Rt) ;
  135. }
  136. push_up(x) ;
  137. }
  138. int Query(int x ,int Le ,int Rt) // 区间求和
  139. {
  140. if(T[x].le == Le && T[x].rt == Rt)
  141. {
  142. if(!T[x].flag) T[x].sum = 0 ;
  143. int temp = T[x].sum ;
  144. T[x].sum = 0 ;
  145. T[x].flag = false ;
  146. return temp ;
  147. }
  148. push_down(x) ;
  149. int Mid = (T[x].le + T[x].rt)>>1 ;
  150. if(Le > Mid) return Query(R(x) ,Le ,Rt) ;
  151. else if(Rt <= Mid) return Query(L(x) ,Le ,Rt) ;
  152. else
  153. return Query(L(x) ,Le ,Mid) + Query(R(x) ,Mid+1 ,Rt) ;
  154. }
  155. void LCA(int u ,int v)
  156. {
  157. int temp ;
  158. while(top[u] != top[v])
  159. {
  160. if(dep[top[u]] < dep[top[v]])
  161. swap(u ,v) ;
  162. temp = Query(1 ,ti[top[u]] ,ti[u]) ;
  163. update(1 ,ti[top[u]] ,ti[u]) ;
  164. sum += temp ;
  165. u = father[top[u]] ;
  166. }
  167. if(dep[u] > dep[v])
  168. swap(u ,v) ;
  169. temp = Query(1 ,ti[u] ,ti[v]) ;
  170. update(1 ,ti[u] ,ti[v]) ;
  171. sum += temp ;
  172. }
  173. int main()
  174. {
  175. //freopen("input.txt" ,"r" ,stdin) ;
  176. int Tx ,u ,v ;
  177. scanf("%d" ,&Tx) ;
  178. while(Tx--)
  179. {
  180. init() ;
  181. scanf("%d" ,&n) ;
  182. for(int i = 1 ;i <= n ; ++i)
  183. scanf("%d" ,&V[i]) ;
  184. for(int i = 1 ;i < n ; ++i)
  185. {
  186. scanf("%d%d" ,&u ,&v) ;
  187. add_edge(u ,v) ;
  188. }
  189. dep[1] = siz[0] = 0 ;
  190. dfs_find(1 ,1) ;
  191. idx = 1 ;
  192. dfs_time(1 ,1) ;
  193. build(1 ,1 ,n) ;
  194. for(int i = 1 ;i <= n ; ++i)
  195. insert(1 ,ti[i] ,V[i]) ;
  196. scanf("%d" ,&m) ;
  197. int type ;
  198. for(int i = 0 ;i < m ; ++i)
  199. {
  200. scanf("%d%d" ,&type ,&u) ;
  201. if(type == 1)
  202. {
  203. scanf("%d" ,&v) ;
  204. LCA(u ,v) ;
  205. }
  206. else if(type == 2)
  207. {
  208. bool flag = Findpos(1 ,ti[u]) ;
  209. if(!flag) // 如果为 false 说明已经被占据
  210. {
  211. insert(1 ,ti[u] ,V[u]) ;
  212. sum -= V[u] ;
  213. }
  214. }
  215. else if(type == 3)
  216. {
  217. sum += Query(1 ,ti[u] ,to[u]) ;
  218. update(1 ,ti[u] ,to[u]) ;
  219. }
  220. cout<<sum<<endl ;
  221. }
  222. }
  223. return 0 ;
  224. }


文章来源: blog.csdn.net,作者:Linux猿,版权归原作者所有,如需转载,请联系作者。

原文链接:blog.csdn.net/nyist_zxp/article/details/45457685

【版权声明】本文为华为云社区用户转载文章,如果您发现本社区中有涉嫌抄袭的内容,欢迎发送邮件进行举报,并提供相关证据,一经查实,本社区将立刻删除涉嫌侵权内容,举报邮箱: cloudbbs@huaweicloud.com
  • 点赞
  • 收藏
  • 关注作者

评论(0

0/1000
抱歉,系统识别当前为高风险访问,暂不支持该操作

全部回复

上滑加载中

设置昵称

在此一键设置昵称,即可参与社区互动!

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。