LCT维护生成树

先按照a的权值把边排序,离线维护b的最小生成树。

将a排序后,依次动态加边,我们只需要关注b的值。要保证1-n花费最少,两点间的b值肯定是越小越好,所以我们可以考虑以b为关键字维护最小生成树。

对于新加的边b,如果1-n已经联通,需要更新答案

#include <bits/stdc++.h>
#define INF 0x3f3f3f3f
#define full(a, b) memset(a, b, sizeof a)
using namespace std;
typedef long long ll;
inline int lowbit(int x){ return x & (-x); }
inline int read(){int X = 0, w = 0; char ch = 0;while(!isdigit(ch)) { w |= ch == '-'; ch = getchar(); }while(isdigit(ch)) X = (X << 3) + (X << 1) + (ch ^ 48), ch = getchar();return w ? -X : X;
}
inline int gcd(int a, int b){ return a % b ? gcd(b, a % b) : b; }
inline int lcm(int a, int b){ return a / gcd(a, b) * b; }
template<typename T>
inline T max(T x, T y, T z){ return max(max(x, y), z); }
template<typename T>
inline T min(T x, T y, T z){ return min(min(x, y), z); }
template<typename A, typename B, typename C>
inline A fpow(A x, B p, C lyd){A ans = 1;for(; p; p >>= 1, x = 1LL * x * x % lyd)if(p & 1)ans = 1LL * x * ans % lyd;return ans;
}
const int N = 400005;
int n, m, tot, ans, mx[N], fa[N], w[N], ch[N][2], rev[N], id[N], st[N];
struct Edge {int v, u, a, b;bool operator < (const Edge &rhs) const {return a < rhs.a;}
} e[N];int newNode(int v){++tot;w[tot] = mx[tot] = v, id[tot] = tot;fa[tot] = ch[tot][0] = ch[tot][1] = 0;return tot;
}bool isRoot(int x){return ch[fa[x]][0] != x && ch[fa[x]][1] != x;
}void reverse(int x){rev[x] ^= 1, swap(ch[x][0], ch[x][1]);
}void push_up(int x){int l = ch[x][0], r = ch[x][1];mx[x] = w[x], id[x] = x;if(mx[l] > mx[x]) mx[x] = mx[l], id[x] = id[l];if(mx[r] > mx[x]) mx[x] = mx[r], id[x] = id[r];
}void push_down(int x){if(rev[x]){rev[x] ^= 1;if(ch[x][0]) reverse(ch[x][0]);if(ch[x][1]) reverse(ch[x][1]);}
}void rotate(int x){int y = fa[x], z = fa[y], p = (ch[y][1] == x) ^ 1;ch[y][p^1] = ch[x][p], fa[ch[x][p]] = y;if(!isRoot(y)) ch[z][ch[z][1] == y] = x;fa[x] = z, fa[y] = x, ch[x][p] = y;push_up(y), push_up(x);
}void splay(int x){int pos = 0; st[++pos] = x;for(int i = x; !isRoot(i); i = fa[i]) st[++pos] = fa[i];while(pos) push_down(st[pos--]);while(!isRoot(x)){int y = fa[x], z = fa[y];if(!isRoot(y)) (ch[y][0] == x) ^ (ch[z][0] == y) ? rotate(x) : rotate(y);rotate(x);}push_up(x);
}void access(int x){for(int p = 0; x; p = x, x = fa[x])splay(x), ch[x][1] = p, push_up(x);
}void makeRoot(int x){access(x), splay(x), reverse(x);
}void link(int x, int y){makeRoot(x);fa[x] = y;
}int findRoot(int x){access(x), splay(x);while(ch[x][0]) push_down(x), x = ch[x][0];splay(x);return x;
}void split(int x, int y){makeRoot(x), access(y), splay(y);
}bool isConnect(int x, int y){makeRoot(x);return findRoot(y) == x;
}int main(){ans = INF;n = read(), m = read();for(int i = 1; i <= n; i ++) newNode(0);for(int i = 0; i < m; i ++){e[i].u = read(), e[i].v = read(), e[i].a = read(), e[i].b = read();}sort(e, e + m);for(int i = 0; i < m; i ++){int u = e[i].u, v = e[i].v, t = newNode(e[i].b);if(!isConnect(u, v)) link(u, t), link(t, v);else{split(u, v);if(e[i].b > mx[v]) continue;int tmp = id[v]; splay(tmp);fa[ch[tmp][0]] = fa[ch[tmp][1]] = 0;ch[tmp][0] = ch[tmp][1] = 0;link(u, t), link(t, v);}if(isConnect(1, n)){split(1, n);ans = min(ans, mx[n] + e[i].a);}}printf(ans == INF ? "-1\n" : "%d\n", ans);return 0;
}

转载于:https://www.cnblogs.com/onionQAQ/p/10738436.html

BZOJ 3669 魔法森林相关推荐

  1. BZOJ 3669: [Noi2014]魔法森林( LCT )

    排序搞掉一维, 然后就用LCT维护加边MST. O(NlogN) ------------------------------------------------------------------- ...

  2. BZOJ 3669 . JZOJ 3754. 【NOI2014】魔法森林

    Description 为了得到书法大家的真传,小 E 同学下定决心去拜访住在魔法森林中的隐士.魔法森林可以被看成一个包含 n 个节点 m 条边的无向图,节点标号为1,2,3, - , n,边标号为 ...

  3. 【BZOJ】【3669】【NOI2014】魔法森林

    LCT动态维护MST LCT动态维护MST 我们可以枚举a,然后找从1到n的一条路径使得:这条路径上的b的最大值最小.这个路径肯定在MST上--所以枚举一遍所有的边,动态维护一个关于b值的MST即可. ...

  4. 【bzoj 3669】[Noi2014]魔法森林

    Description 为了得到书法大家的真传,小E同学下定决心去拜访住在魔法森林中的隐士.魔法森林可以被看成一个包含个N节点M条边的无向图,节点标号为1..N,边标号为1..M.初始时小E同学在号节 ...

  5. 图论 BZOJ 3669 [Noi2014]魔法森林

    Description 为了得到书法大家的真传,小E同学下定决心去拜访住在魔法森林中的隐士.魔法森林可以被看成一个包含个N节点M条边的无向图,节点标号为1..N,边标号为1..M.初始时小E同学在号节 ...

  6. [BZOJ]3669: [Noi2014]魔法森林 lct

    Description 为了得到书法大家的真传,小E同学下定决心去拜访住在魔法森林中的隐士.魔法森林可以被看成一个包含个N节点M条边的无向图,节点标号为1..N,边标号为1..M.初始时小E同学在号节 ...

  7. 【BZOJ3669】【NOI2014】魔法森林 (spfa动态队列加点算法)

    3669: [Noi2014]魔法森林 Time Limit: 30 Sec  Memory Limit: 512 MB Submit: 254  Solved: 140 [Submit][Statu ...

  8. 刷题总结——魔法森林(bzoj3669)

    题目: Description 为了得到书法大家的真传,小E同学下定决心去拜访住在魔法森林中的隐士.魔法森林可以被看成一个包含个N节点M条边的无向图,节点标号为1..N,边标号为1..M.初始时小E同 ...

  9. [LCT动态树] [NOI2014]魔法森林,[ZJOI2018]历史

    [NOI2014] 魔法森林 题目 按照aaa精灵从小到大排序 按顺序插入每一条边 加入第iii条边后的最小代价为a[i]a[i]a[i]加上从111到nnn的所有路径中最大bbb最小的路径代价 维护 ...

最新文章

  1. bash: vue: command not found
  2. JVM自动内存管理:对象判定和回收算法
  3. php面向对象分页,PHP基于面向对象封装的分页类示例
  4. 如何了解这个世界(社会)运行的本质
  5. php显示表格,php – 显示所有表格行
  6. Java开发和运行环境的搭建(详细教程)
  7. 得物App获得2020“年度新经济企业”奖项
  8. webpack 基本功能和原理
  9. 线性代数的本质-B站视频
  10. [bzoj5158][Tjoi2014]Alice and Bob
  11. Netscreen ×××配置(一)---基于策略的点到点×××设置
  12. 《程序员面试宝典》.pdf
  13. 将Windows的IE11降级为原先的IE版本
  14. 身份证文字信息及人脸图片采集
  15. 解决AndroidStudio报错问题:Missing essential plugin
  16. 图像特征与描述(2)
  17. Unity网格编程篇(二) 非常详细的Mesh编程入门文章
  18. 管理者的人品的重要性
  19. 支持向量机(SVM)算法原理
  20. MySQL:由于找不到VCRUNTIME140_1.dll,无法继续执行代码。重新安装程序可能会解决此问题

热门文章

  1. 常考数据结构与算法:最长公共子串
  2. oracle 创建带参数的视图
  3. 02-CSS基础与进阶-day10_2018-09-14-20-38-10
  4. python学习之老男孩python全栈第九期_day002作业
  5. 基于rman 全备+归档在线搭建DG
  6. Hadoop1.x之hdfs集群搭建
  7. 定时器工作原理及初值快速计算
  8. 只读域控制器RODC的安装
  9. centos中python报错 SyntaxError: Non-ASCII character ‘\xe7‘ in file car-training.py on line 7, but no en
  10. formZ Pro 9(3D绘图软件)中文版