这个题目是一个比较裸的树剖题,很好写。

http://acm.hdu.edu.cn/showproblem.php?pid=3966

  1. #include <cstdio>
  2. #include <cstring>
  3. #include <algorithm>
  4. #include <iostream>
  5. #include <algorithm>
  6. #include <cstdlib>
  7. #include <vector>
  8. #include <stack>
  9. #include <map>
  10. #include <string>
  11. #define inf 0x3f3f3f3f
  12. #define inf64 0x3f3f3f3f3f3f3f3f
  13. using namespace std;
  14. typedef long long ll;
  15. const int maxn = 5e4 + 10;
  16.  
  17. int f[maxn];//f 保存u的父亲节点
  18. int dep[maxn];//dep保存节点u 的深度
  19. int siz[maxn];//siz保存以u为根的子节点的个数
  20. int son[maxn];//son 保存u的重儿子
  21. int rk[maxn];//rk当前dfs序在树中所对应的节点
  22. int top[maxn];// top保存当前结点所在链的顶端结点
  23. int id[maxn];//dfs的执行顺序
  24.  
  25. int a[maxn];
  26. int n;
  27. ll sum[maxn * 4], lazy[maxn * 4];
  28. //------------------线段树部分---------------//
  29. void push_up(int id) {
  30. sum[id] = sum[id << 1] + sum[id << 1 | 1];
  31. // printf("sum[%d]=%d sum[%d]=%d\n", id << 1, sum[id << 1], id << 1 | 1, sum[id << 1 | 1]);
  32. // printf("sum[%d]=%d\n", id, sum[id]);
  33. }
  34.  
  35. void build(int id, int l, int r) {
  36. lazy[id] = 0;
  37. if (l == r) {
  38. sum[id] = a[rk[l]];
  39. // printf("id=%d sum=%d\n", id, sum[id]);
  40. return;
  41. }
  42. int mid = (l + r) >> 1;
  43. build(id << 1, l, mid);
  44. build(id << 1 | 1, mid + 1, r);
  45. push_up(id);
  46. }
  47.  
  48. void push_down(int id, int len1, int len2) {
  49. if (lazy[id] == 0) return;
  50. sum[id << 1] += lazy[id] * len1;
  51. lazy[id << 1] += lazy[id];
  52.  
  53. sum[id << 1 | 1] += lazy[id] * len2;
  54. lazy[id << 1 | 1] += lazy[id];
  55.  
  56. lazy[id] = 0;
  57. }
  58.  
  59. void update(int id, int l, int r, int x, int y, int val) {
  60. // printf("id=%d l=%d r=%d x=%d y=%d val=%d\n", id, l, r, x, y, val);
  61. if (x <= l && y >= r) {
  62. // printf("id=%d sum=%d\n", id, sum[id]);
  63. sum[id] += val * (r - l + 1);
  64. lazy[id] += val;
  65. // printf("%d\n", sum[id]);
  66. return;
  67. }
  68. int mid = (l + r) >> 1;
  69. push_down(id, mid - l + 1, r - mid);
  70. if (x <= mid) update(id << 1, l, mid, x, y, val);
  71. if (y > mid) update(id << 1 | 1, mid + 1, r, x, y, val);
  72. push_up(id);
  73. }
  74.  
  75. ll query(int id, int l, int r, int x, int y) {
  76. if (x <= l && y >= r) return sum[id];
  77. int mid = (l + r) >> 1;
  78. ll ans = 0;
  79. push_down(id, mid - l + 1, r - mid);
  80. if (x <= mid) ans = (ans + query(id << 1, l, mid, x, y));
  81. if (y > mid) ans = (ans + query(id << 1 | 1, mid + 1, r, x, y));
  82. return ans;
  83. }
  84.  
  85. //------------------------树链剖分-------------------//
  86. // int f[maxn];//f 保存u的父亲节点
  87. // int dep[maxn];//dep保存节点u 的深度
  88. // int siz[maxn];//siz保存以u为根的子节点的个数
  89. // int son[maxn];//son 保存u的重儿子
  90. // int rk[maxn];//rk当前dfs序在树中所对应的节点
  91. // int top[maxn];// top保存当前结点所在链的顶端结点
  92. // int id[maxn];//dfs的执行顺序
  93. struct node {
  94. int v, nxt;
  95. node(int v = 0, int nxt = 0) :v(v), nxt(nxt) {}
  96. }ex[maxn];
  97. int head[maxn], cnt = 0, tot;
  98. void init() {
  99. cnt = 0, tot = 0;
  100. memset(son, 0, sizeof(son));
  101. memset(head, -1, sizeof(head));
  102. }
  103. void add(int u, int v) {
  104. ex[cnt] = node(v, head[u]);
  105. head[u] = cnt++;
  106. ex[cnt] = node(u, head[v]);
  107. head[v] = cnt++;
  108. }
  109.  
  110. void dfs1(int u, int fa, int depth) {
  111. f[u] = fa; dep[u] = depth; siz[u] = 1;
  112. for (int i = head[u]; i != -1; i = ex[i].nxt) {
  113. int v = ex[i].v;
  114. if (v == fa) continue;
  115. dfs1(v, u, depth + 1);
  116. siz[u] += siz[v];
  117. if (siz[v] > siz[son[u]]) son[u] = v;
  118. }
  119. }
  120.  
  121. void dfs2(int u, int t) {
  122. top[u] = t;
  123. id[u] = ++tot;//标记dfs序
  124. rk[tot] = u;//序号tot对应的结点u
  125. if (!son[u]) return;
  126. dfs2(son[u], t);
  127. /*我们选择优先进入重儿子来保证一条重链上各个节点dfs序连续,
  128. 一个点和它的重儿子处于同一条重链,所以重儿子所在重链的顶端还是t*/
  129. for (int i = head[u]; i != -1; i = ex[i].nxt) {
  130. int v = ex[i].v;
  131. if (v != son[u] && v != f[u]) dfs2(v, v);//一个点位于轻链底端,那么它的top必然是它本身
  132. }
  133. }
  134.  
  135. void update2(int x, int y, int z)//修改x到y路径的值
  136. {
  137. while (top[x] != top[y])//不在同一条链上
  138. {
  139. if (dep[top[x]] < dep[top[y]]) swap(x, y);//x为深度大的链
  140. update(1, 1, n, id[top[x]], id[x], z);//x为深度大的链
  141. x = f[top[x]];//深度大的向上跳
  142. }
  143. if (dep[x] > dep[y]) swap(x, y); //这里x和y在同一条链
  144. update(1, 1, n, id[x], id[y], z); //x和y这条链的更新
  145. }
  146.  
  147. ll query2(int x, int y) {
  148. ll ret = 0;
  149. while (top[x] != top[y]) {
  150. if (dep[top[x]] < dep[top[y]]) swap(x, y);
  151. ret = (ret + query(1, 1, n, id[top[x]], id[x]));
  152. x = f[top[x]];
  153. }
  154. if (dep[x] > dep[y]) swap(x, y);
  155. ret = (ret + query(1, 1, n, id[x], id[y]));
  156. return ret;
  157. }
  158.  
  159. //------------------树链剖分结束-------------------//
  160.  
  161. int main() {
  162. int m, p;
  163. while (scanf("%d%d%d", &n, &m, &p) != EOF) {
  164. init();
  165. for (int i = 1; i <= n; i++) scanf("%d", &a[i]);
  166. for (int i = 1; i < n; i++) {
  167. int u, v;
  168. scanf("%d%d", &u, &v);
  169. add(u, v);
  170. }
  171. dfs1(1, -1, 1), dfs2(1, 1);
  172. build(1, 1, n);
  173. while (p--) {
  174. char s[10];
  175. int l, r, k;
  176. scanf("%s", s);
  177. if (s[0] == 'I') {
  178. scanf("%d%d%d", &l, &r, &k);
  179. update2(l, r, k);
  180. }
  181. else if (s[0] == 'D') {
  182. scanf("%d%d%d", &l, &r, &k);
  183. update2(l, r, -k);
  184. }
  185. else {
  186. scanf("%d", &l);
  187. ll ans = query(1, 1, n, id[l], id[l]);
  188. printf("%lld\n", ans);
  189. }
  190. }
  191. }
  192. return 0;
  193. }

  

A - Aragorn's Story HDU - 3966 树剖裸题的更多相关文章

  1. Aragorn's Story HDU - 3966 -树剖模板

    HDU - 3966 思路 :树链剖分就是可以把一个路径上的点映射成几段连续的区间上.这样对于连续的区间可以用线段树维护, 对于每一段连续的区间都可以通过top [ ]数组很快的找到这段连续区间的头. ...

  2. 树剖裸题——BZOJ1036 树的统计

    #include<cstring> #include<cmath> #include<algorithm> #include<cstdio> #defi ...

  3. HDU 3966 (树链剖分+线段树)

    Problem Aragorn's Story (HDU 3966) 题目大意 给定一颗树,有点权. 要求支持两种操作,将一条路径上的所有点权值增加或减少ai,询问某点的权值. 解题分析 树链剖分模板 ...

  4. AC日记——Aragorn's Story HDU 3966

    Aragorn's Story Time Limit: 10000/3000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) ...

  5. 洛谷树剖模板题 P3384 | 树链剖分

    原题链接 对于以u为根的子树,后代节点的dfn显然比他的dfn大,我们可以记录一下回溯到u的dfn,显然这两个dfn构成了一个连续区间,代表u及u的子树 剩下的就和树剖一样了 #include< ...

  6. 树剖想法题——BZOJ3626

    本来是打算作为树剖练习的最后一题的,结果一直WA. 本来以为是自己写的太丑. 最后发现5w的数据 我开了10w的数组 然而有一个数组要×2 哦,好棒棒. #include<cstring> ...

  7. hdu 3966(树链剖分+线段树区间更新)

    传送门:Problem 3966 https://www.cnblogs.com/violet-acmer/p/9711441.html 学习资料: [1]线段树区间更新:https://blog.c ...

  8. HDU 3966 /// 树链剖分+树状数组

    题意: http://acm.hdu.edu.cn/showproblem.php?pid=3966 给一棵树,并给定各个点权的值,然后有3种操作: I x y z : 把x到y的路径上的所有点权值加 ...

  9. hdu 3966 树链剖分

    思路:树链剖分入门题,我这门入得好苦啊,程序很快写出来了,可是在LCA过程中把update函数里的左右边界位置写反了,一直RE到死. #pragma comment(linker, "/ST ...

随机推荐

  1. spring 管理事务配置时,结果 报错: No Hibernate Session bound to thread, and configuration does not allow creation of non-transactional one here这个异常

    java.lang.IllegalStateException: No Hibernate Session bound to thread, and configuration does not al ...

  2. mysql 的CURDATE() 与 NOW() 的区别

    SELECT CURDATE() 查询出的是当前天的开始时间点,比如今天是 2015.02.03号,那不管我在今天什么时间点查询,结果都是今天的凌晨,即今天的开始的那个时间点,因为它只具体到年月日,没 ...

  3. MAC 上brew 更新 出错

    在MAC上brew update的时候出现报错:Error: /usr/local must be writable! 错误,在该文章中也给出解决办法(sudo chown -R $(whoami) ...

  4. ES5与ES6 this 指向详细解析(箭头函数)

    首先要明白箭头函数的作用: 箭头函数除了让函数的书写变得很简洁,可读性很好外:最大的优点是解决了this执行环境所造成的一些问题.比如:解决了匿名函数this指向的问题(匿名函数的执行环境具有全局性) ...

  5. 用threejs 实现3D物体在浏览器展示

    用threejs 实现3D物体在浏览器展示,通过鼠标平移,缩放,键盘箭头按钮左右移动等功能展示. <!DOCTYPE html> <html> <head> < ...

  6. Julia的基本知识

    知识来源 1.变量.整数和浮点数 Julia和Matllab挺像的,基本的变量,数值定义都差不多,所以就没必要记录了. 2.数学运算 3.函数

  7. 一口气带你踩完五个 List 的大坑,真的是处处坑啊!

    List 可谓是我们经常使用的集合类之一,几乎所有业务代码都离不开 List.既然天天在用,那就没准就会踩中这几个 List 常见坑. 今天我们就来总结这些常见的坑在哪里,捞自己一手,防止后续同学再继 ...

  8. [HarekazeCTF2019] web

    在 buuoj 上看到的这个比赛题目,期间平台关了,就拿了 Dockerfile 本地做了,web 题目感觉还不错 encode_and_encode [100] 打开靶机,前两个页面都是 html ...

  9. jmeter json path espressions学习

    jsonpath表达式可以使用点或者括号来取值 管方详细介绍:https://goessner.net/articles/JsonPath/ $ :表示根对象 @:表示当前对象 . 或者[]:表示子运 ...

  10. 【mybatis xml】数据层框架应用--Mybatis 基于XML映射文件实现数据的CRUD

    使用MyBatis框架进行持久层开发 MyBatis是支持普通SQL查询,存储过程和高级映射的优秀持久层框架. MyBatis消除了几乎所有的JDBC代码和参数的手工设置以及对结果集的检索. MyBa ...