以前一直觉得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)的更多相关文章

  1. bzoj 3669: [Noi2014]魔法森林 动态树

    3669: [Noi2014]魔法森林 Time Limit: 30 Sec  Memory Limit: 512 MBSubmit: 363  Solved: 202[Submit][Status] ...

  2. BZOJ3669: [Noi2014]魔法森林(瓶颈生成树 LCT)

    Time Limit: 30 Sec  Memory Limit: 512 MBSubmit: 3558  Solved: 2283[Submit][Status][Discuss] Descript ...

  3. [bzoj] 3669 NOI2014 魔法森林 || LCT

    原题 copy一篇题解:原链接 将边按照a排序,然后从小到大枚举,加入图中. 在图中用lct维护一棵两点之间b最大值尽量小的生成树. 当加入一条边(u, v)时: 如果(u, v)不连通,则直接加入这 ...

  4. BZOJ 3669 [Noi2014]魔法森林(贪心+LCT)

    [题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=3669 [题目大意] 给出一张图,每条边上有两个值ai和bi,现在需要求出一条1到N的路 ...

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

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

  6. P2387 [NOI2014]魔法森林(LCT)

    P2387 [NOI2014]魔法森林 LCT边权维护经典题 咋维护呢?边化为点,边权变点权. 本题中我们把边对关键字A进行排序,动态维护关键字B的最小生成树 加边后出现环咋办? splay维护最大边 ...

  7. loj2245 [NOI2014]魔法森林 LCT

    [NOI2014]魔法森林 链接 loj 思路 a排序,b做动态最小生成树. 把边拆成点就可以了. uoj98.也许lct复杂度写假了..越卡常,越慢 代码 #include <bits/std ...

  8. bzoj 3669: [Noi2014]魔法森林 (LCT)

    链接:https://www.lydsy.com/JudgeOnline/problem.php?id=3669 题面: 3669: [Noi2014]魔法森林 Time Limit: 30 Sec  ...

  9. [NOI2014]魔法森林 LCT

    题面 [NOI2014]魔法森林 题解 一条路径的代价为路径上的\(max(a[i]) + max(b[i])\),因为一条边同时有$a[i], b[i]$2种权值,直接处理不好同时兼顾到,所以我们考 ...

随机推荐

  1. [题解]bzoj 1861 Book 书架 - Splay

    1861: [Zjoi2006]Book 书架 Time Limit: 4 Sec  Memory Limit: 64 MBSubmit: 1396  Solved: 803[Submit][Stat ...

  2. 如何在Eclipse下安装myeclipse插件

    来自http://www.blogjava.net/show911/archive/2008/04/27/86284.html 下载myeclipse插件 支持eclipse3.1.x, 具体安装步骤 ...

  3. PHP使用GD2库画图,图像无法输出解决方法

    在CS6编辑环境下,使用php的GD2库绘制图像,在浏览器上输出可能会出现图像无法输出的情况,目前发现两个解决方法:方法1:用记事本写,再另存为utf-8 no bom的格式. 方法2:在代码中添加o ...

  4. iOS 之 获取View所在控制器

    1. UIResponder UIViewController *uvc; UIResponder* nextResponder = [self.superview.superview.supervi ...

  5. Voilin 之 握弓

    握弓要像拿杯子,手要圆:整个手型是左倾.

  6. Alamofire 4.0 迁移指南

    Alamofire 4.0 是 Alamofire 最新的一个大版本更新, 一个基于 Swift 的 iOS, tvOS, macOS, watchOS 的 HTTP 网络库. 作为一个大版本更新, ...

  7. react中常用的一些方法

    React.createClass:创建一个ReactClass(组件类),参数是一个对象且必须带有 render 属性方法,该方法必须返回一个封闭的容器(容器内可以有其它不限结构的容器)或 null ...

  8. swift 定位 根据定位到的经纬度转换城市名

    好久没写随笔了   最近这段时间项目有点紧  天天在加班  国庆 一天假都没放  我滴娃娃   好啦  牢骚就不发了  毕竟没有什么毛用    待我那天闲了专门写一篇吐槽的随笔  

  9. C++编程之面向对象的三个基本特征

    面向对象的三个基本特征是:封装.继承.多态. 封装 封装最好理解了.封装是面向对象的特征之一,是对象和类概念的主要特性. 封装,也就是把客观事物封装成抽象的类,并且类可以把自己的数据和方法只让可信的类 ...

  10. java中Integer比较需要注意的问题

    java中Integer比较需要注意的问题 package com.srie.test; import java.util.HashMap; import java.util.Map; public ...