bzoj3669(NOI2014)魔法森林

将边按a排序,一条条加入,若当前的x,y未联通则直接加入这条边,

否则看当前边的b是否小于x,y,路径上最大的b,小于则拆掉该边,加入当前边。

显然可以用lct维护,splay维护b的最大值

把边拆成点权为b的点,点则是点权为0

Cut时找到改点,splay到根,切掉左右儿子

Link时加入新点,把x,y newroot,虚父亲指向新点

//Achen
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<vector>
#include<queue>
#include<ctime>
#include<cmath>
const int N=200000+7; 
typedef long long LL;
using namespace std;

template<typename T> void read(T &x) {
    char ch=getchar(); x=0; T f=1;
    while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
    if(ch=='-') f=-1,ch=getchar();
    for(;ch>='0'&&ch<='9';ch=getchar()) x=x*10+ch-'0'; x*=f;
}

int n,m,ecnt;
struct edge{
    int x,y,a,b;
    edge(){}
    edge(int x,int y,int a,int b):x(x),y(y),a(a),b(b){} 
    friend bool operator <(const edge&A,const edge&B) {
        return A.a<B.a||(A.a==B.a&&A.b<B.b);
    }
}e[N];

int tot,ch[N][2],p[N],v[N],bb[N],flip[N];
#define lc ch[x][0]
#define rc ch[x][1]
void update(int x) {v[x]=max(max(bb[x],v[lc]),v[rc]);}

int isroot(int x) {return (ch[p[x]][0]!=x&&ch[p[x]][1]!=x);}

void down(int x) {
    if(!flip[x]) return;
    swap(lc,rc); 
    flip[lc]^=1;
    flip[rc]^=1;
    flip[x]^=1; 
}

void rotate(int x) {
    int y=p[x],z=p[y],l=(x==ch[y][1]),r=l^1;
    if(!isroot(y)) ch[z][y==ch[z][1]]=x; p[x]=z;
    ch[y][l]=ch[x][r]; p[ch[x][r]]=y;
    ch[x][r]=y; p[y]=x;
    update(y); update(x);
}

void splay(int x) {
    static int g[N],top=0,tp;
    for(tp=x;!isroot(tp);tp=p[tp]) g[++top]=tp;
    g[++top]=tp;
    while(top) down(g[top--]);
    for(;!isroot(x);rotate(x)) {
        int y=p[x],z=p[y];
        if(!isroot(y))
            ((x==ch[y][1])^(y==ch[z][1]))?rotate(x):rotate(y); 
    } 
}

void access(int x) {
    for(int t=0;x;x=p[t=x]) {
        splay(x);
        ch[x][1]=t;
        update(x);
    }
}

int findroot(int x) {
    access(x);
    splay(x);
    while(lc) x=lc;
    return x;
}

int newroot(int x) {
    access(x);
    splay(x);
    flip[x]^=1;
}

int qry(int x,int y) {
    newroot(x);
    access(y);
    splay(y);
    return v[y];
} 

void lik(int x,int y,int w) {
    int z=++tot;
    bb[z]=v[z]=w;
    newroot(x);
    p[x]=z;
    newroot(y);
    p[y]=z;
}

int find(int x,int y) {
    if(bb[x]==y) return x;
    if(v[lc]==y) return find(lc,y);
    else return find(rc,y);
}

int cut(int x) {
    splay(x);
    p[ch[x][0]]=0;
    p[ch[x][1]]=0;
}

void work() {
    int ans=1e9+7;
    sort(e+1,e+m+1);
    for(int i=1;i<=m;i++) {
        int u=e[i].x,v=e[i].y,a=e[i].a,b=e[i].b;
        if(findroot(u)!=findroot(v)) lik(u,v,b);
        else {
            int tp=qry(u,v);
            if(b<tp) {
                int tpp=find(v,tp);
                cut(tpp);
                lik(u,v,b);
            }
        }
        if(findroot(1)==findroot(n)) {
            int tp=qry(1,n);
            ans=min(ans,a+tp);
        }
    } 
    if(ans==1e9+7) ans=-1;
    printf("%d\n",ans);
}

void init() {
    read(n); read(m);
    tot=n;
    for(int i=1;i<=m;i++) {
        int x,y,a,b;
        read(x); read(y); read(a); read(b);
        if(x!=y) e[++ecnt]=edge(x,y,a,b); 
    }
    m=ecnt;
}

int main() {
#ifdef DEBUG
    freopen(".in","r",stdin);
    freopen(".out","w",stdout);
#endif
    init();
    work();
    return 0;
}
/*
4 5
1 2 19 1
2 3 8 12
2 4 12 15
1 3 17 8
3 4 1 17

3 1
1 2 1 1
*/
View Code

 

posted @ 2017-12-11 07:28  啊宸  阅读(192)  评论(0编辑  收藏  举报