[Noi2014]魔法森林( 动态mst lct)
以前一直觉得lct特别难写,自从学了丽洁姐的lct之后,觉得lct居然能这么短,这个主程序能40行左右解决~~~~
这道嘛~~虽说能用spfa解决,但还是写下lct吧
把边按a值排序后一条一条插入并维护bmx值就行了= =
不过还是得膜拜云神在考场上a了这道题,我花了好久,自己还是太弱
NOIP就要来了,在这里祝自己RP++吧
毕竟升高二后,每天几乎都是怀着梦想睡觉的
希望这一年能不辜负我的愿望吧
向着最高点,冲吧!!!
CODE:
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<map>
using namespace std;
#define maxn 101000
struct node{
node *ch[2],*p,*mx;bool rec;int val;
node(int v=0):val(v) {
rec=false;
p=ch[0]=ch[1]=NULL;
mx=this;
}
inline bool d();
inline bool isroot();
inline void update();
inline void relax();
}t[maxn];
inline bool node::d(){return p->ch[1]==this;}
inline bool node::isroot(){
return (p==NULL)||(p->ch[1]!=this&&p->ch[0]!=this);
}
inline void node::relax(){
if (!this) return;
if (!rec) return ;
if (ch[1]) ch[1]->rec^=1;
if (ch[0]) ch[0]->rec^=1;
swap(ch[1],ch[0]);
rec=0;
}
inline void node::update(){
if (!this) return ;
mx=this;
if (ch[0]&&ch[0]->mx->val > mx->val ) mx=ch[0]->mx;
if (ch[1]&&ch[1]->mx->val > mx->val ) mx=ch[1]->mx;
}
inline void rotate(node*u) {
node *v=u->p;
v->relax();u->relax();
bool d=u->d();
if (!v->isroot()) v->p->ch[v->d()]=u;
u->p=v->p;
if (u->ch[d^1]) u->ch[d^1]->p=v;
v->ch[d]=u->ch[d^1];
u->ch[d^1]=v;
v->p=u;
v->update(),u->update();
}
inline void splay(node *u) {
u->relax();
while (!u->isroot()){
if (u->p->isroot()) rotate(u);
else if (u->d()!=u->p->d()) rotate(u),rotate(u);
else rotate(u->p),rotate(u);
}
u->update();
}
inline node* expose(node *u) {
node *v;
for (v=NULL;u;v=u,u=u->p){
splay(u);
u->ch[1]=v;
u->update();
}
return v;
}
inline void evert(node *u) {
expose(u)->rec^=1;
splay(u);
}
inline void link(node *u,node *v) {
evert(u);
u->p=v;
expose(u);
}
inline void cut(node *u,node *v) {
evert(u);
expose(v);
splay(v);
v->ch[0]=u->p=NULL;
v->update();u->update();
}
inline int query(node *u,node *v) {
evert(u);
expose(v);
splay(v);
return v->mx->val;
}
int fa[maxn];
inline int find(int x) {
if (fa[x]!=x) fa[x]=find(fa[x]);
return fa[x];
}
int x[maxn],y[maxn],a[maxn],b[maxn],id[maxn];
bool cmp(int p,int q) {return a[p]<a[q];}
map<node*,int> ma;
int n,m;
int main(){
scanf("%d%d",&n,&m);
for (int i=1;i<=n;i++) fa[i]=i;
for (int i=1;i<=m;i++) {
scanf("%d%d%d%d",x+i,y+i,a+i,b+i);
id[i]=i;
}
sort(id+1,id+1+m,cmp);
int ans=0x7fffffff;
for (int i=1;i<=m;i++) {
if (x[id[i]]==y[id[i]]) continue;
node *u=&t[x[id[i]]],*v=&t[y[id[i]]];
if (find(x[id[i]])!=find(y[id[i]])) {
node *e=new node(b[id[i]]);
ma[e]=id[i];
link(u,e);
link (v,e);
fa[find(x[id[i]])]=find(y[id[i]]);
}else {
evert(u);
expose(v);
splay(v);
if (v->mx->val>b[id[i]]){
node *e=v->mx;
splay(e);
node *_u=&t[x[ma[e]]],*_v=&t[y[ma[e]]];
cut(e,_u);
cut(e,_v);
delete e;
e=new node(b[id[i]]);
ma[e]=id[i];
link (u,e);
link (v,e);
}
}
if (find(1)==find(n)) ans=min(ans,a[id[i]]+query(&t[1],&t[n]));
}
if (ans==0x7fffffff) ans=-1;
printf("%d\n",ans);
return 0;
}
[Noi2014]魔法森林( 动态mst lct)的更多相关文章
- bzoj 3669: [Noi2014]魔法森林 动态树
3669: [Noi2014]魔法森林 Time Limit: 30 Sec Memory Limit: 512 MBSubmit: 363 Solved: 202[Submit][Status] ...
- BZOJ3669: [Noi2014]魔法森林(瓶颈生成树 LCT)
Time Limit: 30 Sec Memory Limit: 512 MBSubmit: 3558 Solved: 2283[Submit][Status][Discuss] Descript ...
- [bzoj] 3669 NOI2014 魔法森林 || LCT
原题 copy一篇题解:原链接 将边按照a排序,然后从小到大枚举,加入图中. 在图中用lct维护一棵两点之间b最大值尽量小的生成树. 当加入一条边(u, v)时: 如果(u, v)不连通,则直接加入这 ...
- BZOJ 3669 [Noi2014]魔法森林(贪心+LCT)
[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=3669 [题目大意] 给出一张图,每条边上有两个值ai和bi,现在需要求出一条1到N的路 ...
- BZOJ 3669: [Noi2014]魔法森林( LCT )
排序搞掉一维, 然后就用LCT维护加边MST. O(NlogN) ------------------------------------------------------------------- ...
- P2387 [NOI2014]魔法森林(LCT)
P2387 [NOI2014]魔法森林 LCT边权维护经典题 咋维护呢?边化为点,边权变点权. 本题中我们把边对关键字A进行排序,动态维护关键字B的最小生成树 加边后出现环咋办? splay维护最大边 ...
- loj2245 [NOI2014]魔法森林 LCT
[NOI2014]魔法森林 链接 loj 思路 a排序,b做动态最小生成树. 把边拆成点就可以了. uoj98.也许lct复杂度写假了..越卡常,越慢 代码 #include <bits/std ...
- bzoj 3669: [Noi2014]魔法森林 (LCT)
链接:https://www.lydsy.com/JudgeOnline/problem.php?id=3669 题面: 3669: [Noi2014]魔法森林 Time Limit: 30 Sec ...
- [NOI2014]魔法森林 LCT
题面 [NOI2014]魔法森林 题解 一条路径的代价为路径上的\(max(a[i]) + max(b[i])\),因为一条边同时有$a[i], b[i]$2种权值,直接处理不好同时兼顾到,所以我们考 ...
随机推荐
- [题解]bzoj 1861 Book 书架 - Splay
1861: [Zjoi2006]Book 书架 Time Limit: 4 Sec Memory Limit: 64 MBSubmit: 1396 Solved: 803[Submit][Stat ...
- 如何在Eclipse下安装myeclipse插件
来自http://www.blogjava.net/show911/archive/2008/04/27/86284.html 下载myeclipse插件 支持eclipse3.1.x, 具体安装步骤 ...
- PHP使用GD2库画图,图像无法输出解决方法
在CS6编辑环境下,使用php的GD2库绘制图像,在浏览器上输出可能会出现图像无法输出的情况,目前发现两个解决方法:方法1:用记事本写,再另存为utf-8 no bom的格式. 方法2:在代码中添加o ...
- iOS 之 获取View所在控制器
1. UIResponder UIViewController *uvc; UIResponder* nextResponder = [self.superview.superview.supervi ...
- Voilin 之 握弓
握弓要像拿杯子,手要圆:整个手型是左倾.
- Alamofire 4.0 迁移指南
Alamofire 4.0 是 Alamofire 最新的一个大版本更新, 一个基于 Swift 的 iOS, tvOS, macOS, watchOS 的 HTTP 网络库. 作为一个大版本更新, ...
- react中常用的一些方法
React.createClass:创建一个ReactClass(组件类),参数是一个对象且必须带有 render 属性方法,该方法必须返回一个封闭的容器(容器内可以有其它不限结构的容器)或 null ...
- swift 定位 根据定位到的经纬度转换城市名
好久没写随笔了 最近这段时间项目有点紧 天天在加班 国庆 一天假都没放 我滴娃娃 好啦 牢骚就不发了 毕竟没有什么毛用 待我那天闲了专门写一篇吐槽的随笔
- C++编程之面向对象的三个基本特征
面向对象的三个基本特征是:封装.继承.多态. 封装 封装最好理解了.封装是面向对象的特征之一,是对象和类概念的主要特性. 封装,也就是把客观事物封装成抽象的类,并且类可以把自己的数据和方法只让可信的类 ...
- java中Integer比较需要注意的问题
java中Integer比较需要注意的问题 package com.srie.test; import java.util.HashMap; import java.util.Map; public ...