传送门

WC这个题真是丧心病狂啊,就是想学习一下怎么处理边权,给我来了这么一个破题!

ORZ hzwer 临摹黄学长代码233 但还是复杂的一匹

理一下思路吧

题目大意:给定一个无向图,多次删除图中的某一条边,求两点间路径最大值的最小值

求两点间的路径最大值的最小值的话,可以求最小生成树,那么这个值就是最小生成树上两点间路径上的最大值

但是题目要求是删除边,LCT维护最小生成树不支持删边操作,那么就离线处理,倒着加边,用LCT维护。

就是这个离线处理是最恶心的。

来说说如何处理边权,把边也抽象成点,那么这个点就和原来边所连的两个点相连,其中边抽象成的点点权为边权,所连接的两个点点权为 0

给边从小到大排序,那么边所抽象成的点的序号即为 边的序号 + n (n 为原来点的个数)

然后再将边按照它所连的两个点为关键字排序,二分查找确定所删除的边的序号。

再按照从小到大的顺序再排回来,跑 kruskal,依次添加没有被删除的边

最后从后往前处理询问,添加一条边的时候先判断两端点是否连通(当然此题不必,因为题目中说:任何时候我们考虑的水管网络都是连通的,即从任一结点A必有至少一条水管路径通往任一结点B。。所以肯定连通)

然后找出两点间路径边权的最大的,判断最大的边权是否大于要加入的边的边权,如果大于,就删除这条边,连接新边,否则就不加入(维护最小生成树)

具体细节看代码

——代码

  1. #include <cstdio>
  2. #include <cstring>
  3. #include <iostream>
  4. #include <algorithm>
  5. #define N 1500005
  6. #define get(x) (son[f[x]][1] == (x))
  7. #define swap(x, y) ((x) ^= (y) ^= (x) ^= (y))
  8. #define isroot(x) (son[f[x]][0] ^ (x) && son[f[x]][1] ^ (x))
  9.  
  10. int n, m, Q, tot;
  11. int mx[N], rev[N], fa[N], f[N], val[N], son[N][], s[N];
  12.  
  13. struct node
  14. {
  15. int x, y, z, id;
  16. bool d;
  17. }e[N];
  18.  
  19. struct ask
  20. {
  21. int f, x, y, ans, id;
  22. }q[N];
  23.  
  24. inline int read()
  25. {
  26. int x = , f = ;
  27. char ch = getchar();
  28. for(; !isdigit(ch); ch = getchar()) if(ch == '-') f = -;
  29. for(; isdigit(ch); ch = getchar()) x = (x << ) + (x << ) + ch - '';
  30. return x * f;
  31. }
  32.  
  33. inline bool cmp1(node x, node y)
  34. {
  35. return x.z < y.z;
  36. }
  37.  
  38. inline bool cmp2(node x, node y)
  39. {
  40. return x.x < y.x || (x.x == y.x && x.y < y.y);
  41. }
  42.  
  43. inline bool cmp3(node x, node y)
  44. {
  45. return x.id < y.id;
  46. }
  47.  
  48. inline int getf(int x)
  49. {
  50. return x == fa[x] ? x : fa[x] = getf(fa[x]);
  51. }
  52.  
  53. inline int find(int x, int y)
  54. {
  55. int l = , r = m, mid;
  56. while(l <= r)
  57. {
  58. mid = (l + r) >> ;
  59. if(e[mid].x < x || (e[mid].x == x && e[mid].y < y)) l = mid + ;
  60. else if(e[mid].x == x && e[mid].y == y) return mid;
  61. else r = mid - ;
  62. }
  63. }
  64.  
  65. inline void update(int x)
  66. {
  67. if(x)
  68. {
  69. mx[x] = x;
  70. if(son[x][] && val[mx[son[x][]]] > val[mx[x]]) mx[x] = mx[son[x][]];
  71. if(son[x][] && val[mx[son[x][]]] > val[mx[x]]) mx[x] = mx[son[x][]];
  72. }
  73. }
  74.  
  75. inline void pushdown(int x)
  76. {
  77. if(x && rev[x])
  78. {
  79. swap(son[x][], son[x][]);
  80. if(son[x][]) rev[son[x][]] ^= ;
  81. if(son[x][]) rev[son[x][]] ^= ;
  82. rev[x] = ;
  83. }
  84. }
  85.  
  86. inline void rotate(int x)
  87. {
  88. int old = f[x], oldf = f[old], wh = get(x);
  89.  
  90. if(!isroot(old))
  91. son[oldf][get(old)] = x;
  92. f[x] = oldf;
  93.  
  94. son[old][wh] = son[x][wh ^ ];
  95. f[son[old][wh]] = old;
  96.  
  97. son[x][wh ^ ] = old;
  98. f[old] = x;
  99.  
  100. update(old);
  101. update(x);
  102. }
  103.  
  104. inline void splay(int x)
  105. {
  106. int i, fat, t = ;
  107. s[++t] = x;
  108. for(i = x; !isroot(i); i = f[i]) s[++t] = f[i];
  109. for(i = t; i; i--) pushdown(s[i]);
  110. for(; !isroot(x); rotate(x))
  111. if(!isroot(fat = f[x]))
  112. rotate(get(x) ^ get(fat) ? x : fat);
  113. }
  114.  
  115. inline void access(int x)
  116. {
  117. for(int t = ; x; t = x, x = f[x]) splay(x), son[x][] = t, update(x);
  118. }
  119.  
  120. inline void reverse(int x)
  121. {
  122. access(x);
  123. splay(x);
  124. rev[x] ^= ;
  125. }
  126.  
  127. inline void cut(int x, int y)
  128. {
  129. reverse(x);
  130. access(y);
  131. splay(y);
  132. son[y][] = f[x] = ;
  133. update(y);
  134. }
  135.  
  136. inline void link(int x, int y)
  137. {
  138. reverse(x);
  139. f[x] = y;
  140. access(x);
  141. }
  142.  
  143. inline int query(int x, int y)
  144. {
  145. reverse(x);
  146. access(y);
  147. splay(y);
  148. return mx[y];
  149. }
  150.  
  151. int main()
  152. {
  153. int i, j, k, t, x, y, fx, fy;
  154. n = read();
  155. m = read();
  156. Q = read();
  157. for(i = ; i <= n; i++) fa[i] = i;
  158. for(i = ; i <= m; i++)
  159. {
  160. e[i].x = read();
  161. e[i].y = read();
  162. e[i].z = read();
  163. if(e[i].x > e[i].y) swap(e[i].x, e[i].y);
  164. }
  165. std::sort(e + , e + m + , cmp1);
  166. for(i = ; i <= m; i++)
  167. {
  168. e[i].id = i;
  169. val[n + i] = e[i].z;
  170. mx[n + i] = n + i;
  171. }
  172. std::sort(e + , e + m + , cmp2);
  173. for(i = ; i <= Q; i++)
  174. {
  175. q[i].f = read();
  176. q[i].x = read();
  177. q[i].y = read();
  178. if(q[i].x > q[i].y) swap(q[i].x, q[i].y);
  179. if(q[i].f == )
  180. {
  181. t = find(q[i].x, q[i].y);
  182. e[t].d = ;
  183. q[i].id = e[t].id;
  184. }
  185. }
  186. std::sort(e + , e + m + , cmp3);
  187. for(i = ; i <= m; i++)
  188. if(!e[i].d)
  189. {
  190. x = e[i].x;
  191. y = e[i].y;
  192. fx = getf(x);
  193. fy = getf(y);
  194. if(fx ^ fy)
  195. {
  196. fa[fx] = fy;
  197. link(x, i + n);
  198. link(y, i + n);
  199. tot++;
  200. if(tot == n - ) break;
  201. }
  202. }
  203. for(i = Q; i; i--)
  204. {
  205. if(q[i].f == ) q[i].ans = val[query(q[i].x, q[i].y)];
  206. else
  207. {
  208. x = q[i].x;
  209. y = q[i].y;
  210. k = q[i].id;
  211. t = query(x, y);
  212. if(e[k].z < val[t])
  213. {
  214. cut(e[t - n].x, t);
  215. cut(e[t - n].y, t);
  216. link(x, k + n);
  217. link(y, k + n);
  218. }
  219. }
  220. }
  221. for(i = ; i <= Q; i++)
  222. if(q[i].f == )
  223. printf("%d\n", q[i].ans);
  224. return ;
  225. }

[BZOJ2594] [Wc2006]水管局长数据加强版(LCT + kruskal + 离线)的更多相关文章

  1. BZOJ2594 [Wc2006]水管局长数据加强版 LCT kruskal

    欢迎访问~原文出处——博客园-zhouzhendong 去博客园看该题解 题目传送门 - BZOJ2594 题意概括 N个点的图,M条带权边.(N<=100000,M<=1000000) ...

  2. [bzoj2594][Wc2006]水管局长数据加强版 (lct)

    论蒟蒻的自我修养T_T.. 和noi2014魔法森林基本一样...然而数据范围大得sxbk...UPD:这题如果用lct判联通的话可能会被卡到O(mlogm)..所以最好还是用并查集吧 一开始数组开太 ...

  3. BZOJ 2594: [Wc2006]水管局长数据加强版 [LCT kruskal]

    2594: [Wc2006]水管局长数据加强版 Time Limit: 25 Sec  Memory Limit: 128 MBSubmit: 2917  Solved: 918[Submit][St ...

  4. [bzoj2594][Wc2006]水管局长数据加强版——lct+离线

    Brief Description 您有一个无向带权图,您需要支持两种操作. 询问两个点之间的最大权最小路径. 删除一条边. Algorithm Design 我们首先提出一个猜想:最优路径一定在原图 ...

  5. bzoj2594 [Wc2006]水管局长数据加强版——LCT

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=2594 时间倒序一下,就是 魔法森林 那道题: 有个不解的地方,是 access 里面关于 p ...

  6. BZOJ2594: [Wc2006]水管局长数据加强版

    题解: 裸LCT+离线+二分+MST... 代码:(几乎摘抄自hzwer) #include<cstdio> #include<cstdlib> #include<cma ...

  7. BZOJ 2594: [Wc2006]水管局长数据加强版( LCT )

    离线然后就是维护加边的动态MST, Link cut tree秒掉..不过我写+调了好久...时间复杂度O(NlogN + MlogM) ------------------------------- ...

  8. BZOJ2594 [Wc2006]水管局长数据加强版 【LCT维护最小生成树】

    题目 SC省MY市有着庞大的地下水管网络,嘟嘟是MY市的水管局长(就是管水管的啦),嘟嘟作为水管局长的工作就是:每天供水公司可能要将一定量的水从x处送往y处,嘟嘟需要为供水公司找到一条从A至B的水管的 ...

  9. BZOJ 2594: [Wc2006]水管局长数据加强版(kruskal + LCT)

    Description SC省MY市有着庞大的地下水管网络,嘟嘟是MY市的水管局长(就是管水管的啦),嘟嘟作为水管局长的工作就是:每天供水公司可能要将一定量的水从x处送往y处,嘟嘟需要为供水公司找到一 ...

随机推荐

  1. hdu 1043 Eight

    欸我一直以为双向bfs是搜完一半再搜另一半呢,妹想到是两个一起搜 然后队列里放的结构体里不能直接存答案,所以做一个邻接表一样的东西,直接指向需要的字符即可 记录状态用康托展开来hash 以及居然是多组 ...

  2. bzoj 1635: [Usaco2007 Jan]Tallest Cow 最高的牛【差分】

    s[i]为差分后的"i这头牛前有几头比它高",计算答案的时候加成前缀和,假设第一头最高减一下即可 用map记录一下被加过的区间,避免重复 #include<iostream& ...

  3. VUE学习之计算属性computed

    计算属性:computed 先看一下官网的说法 模板内的表达式非常便利,但是设计它们的初衷是用于简单运算的.在模板中放入太多的逻辑会让模板过重且难以维护.例如: <div id="ex ...

  4. 【工具】---- json-server基本使用

    一.概念 在开发过程中,前端通常需要等待后端开发完接口后,再调用接口渲染相应的数据,这会影响开发效率.而json-server的作用就是为了解决前后端并行开发的痛点,在本地模拟后端接口用来测试前端效果 ...

  5. mysql之replace into与 insert into duplicat key for update

    mysql实际应用中,我们在插入数据的时候,经常遇到主键冲突的情况,这是因为库中已经存在相同主键的数据,这时,我们只能更新数据:在判断是更新数据还是插入数据,我们还需要在此之前做一些必要的判断:在my ...

  6. C++帮助文档(自己写的)

    以下所有记录几乎都是摘抄自<C++ primer 5th 中文> auto   类型说明符   P61 特点: 1.     定义的变量必须有初始值 2.     通过初始值来推算变量的类 ...

  7. Java设置全局热键——第三方包jintellitype实现

    Java原生API并不支持为应用程序设置全局热键.要实现全局热键,需要用JNI方式实现,这就涉及到编写C/C++代码,这对于大多数不熟悉C/C++的javaer来说,有点困难.不过幸好,国外有人已经实 ...

  8. 求助:可以使用任何编程工具做成一个控件或组件,使得在VB中能调用并得到摄像头的参数及图片。

    请看下网址上的这个问题,看是否有解决的方式http://www.educity.cn/wenda/338634.html

  9. npm——nrm

    nrm 是镜像地址,方便国内下载 npm i nrm -g nrm ls nrm use taobao // 切换地址

  10. linux环境下为php7装phpredis扩展

    phpredis在php7.php5下都有不同的版本,装岔了可能会编译报错,所以在安装之前请先看下自己的php是啥版本. 我的redis装的是redis3.2.3版本. 用phpinfo()查看安装的 ...