bzoj3669 [Noi2014]魔法森林（并查集+lct）

把所有边按a值从小到大排序。不断往里加边，维护连通性，如果成环，查询目前x,y路径上的b的最大值，如果我当前边的b比最大值还大，则肯定不往里加。否则替换掉最大的b值那条边，连通性不变。如果1和n已经联通，就更新答案，用目前1，n路径上的b的最大值+当前边的a值来更新。这样为什么是正确的呢？考虑你这次替换了一条边进去，分两种情况
1.这条边在目前1~n的路径上，则显然这条边的a值是最大的a值。
2.这条边不在目前1~n的路径上，则显然目前答案的最优值是之前更新过的了。
然后把边权变成点权来做，更方便些。。。

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 150010
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,m,fat[N],fa[N],c[N][2],v[N],mx[N],q[N],ans=inf;bool rev[N];
struct edge{
    int x,y,a,b;
}e[100010];
inline bool cmp(edge a,edge b){return a.a<b.a;}
inline int find(int x){return x==fat[x]?x:fat[x]=find(fat[x]);}
inline bool isroot(int x){
    return x!=c[fa[x]][0]&&x!=c[fa[x]][1];
}
inline void update(int p){
    int l=c[p][0],r=c[p][1];mx[p]=p;
    if(v[mx[l]]>v[mx[p]]) mx[p]=mx[l];
    if(v[mx[r]]>v[mx[p]]) mx[p]=mx[r];
}
inline void pushdown(int p){
    if(!rev[p]) return;rev[p]=0;
    int l=c[p][0],r=c[p][1];
    rev[l]^=1;rev[r]^=1;swap(c[p][0],c[p][1]);
}
inline void rotate(int x){
    int y=fa[x],z=fa[y],l=x==c[y][1],r=l^1;
    if(!isroot(y)) c[z][y==c[z][1]]=x;
    fa[c[x][r]]=y;fa[y]=x;fa[x]=z;
    c[y][l]=c[x][r];c[x][r]=y;update(y);update(x);
}
inline void splay(int x){
    int top=0;q[++top]=x;
    for(int xx=x;!isroot(xx);xx=fa[xx]) q[++top]=fa[xx];
    while(top) pushdown(q[top--]);
    while(!isroot(x)){
        int y=fa[x],z=fa[y];
        if(!isroot(y)){
            if(x==c[y][1]^y==c[z][1]) rotate(x);
            else rotate(y);
        }rotate(x);
    }
}
inline void access(int x){
    int y=0;
    while(x){
        splay(x);c[x][1]=y;update(x);y=x;x=fa[x];
    }
}
inline void makeroot(int x){
    access(x);splay(x);rev[x]^=1;
}
inline void link(int x,int y){
    makeroot(x);fa[x]=y;
}
inline void cut(int x,int y){
    makeroot(x);access(y);splay(y);c[y][0]=fa[x]=0;update(y);
}
inline int qmax(int x,int y){
    makeroot(x);access(y);splay(y);return mx[y];
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();
    for(int i=1;i<=n;++i) fat[i]=i;
    for(int i=1;i<=m;++i) e[i].x=read(),e[i].y=read(),e[i].a=read(),e[i].b=read();
    sort(e+1,e+m+1,cmp);
    for(int i=1;i<=m;++i){
        int x=e[i].x,y=e[i].y;v[n+i]=e[i].b;
        int xx=find(x),yy=find(y);
        if(xx!=yy) fat[xx]=yy,link(x,n+i),link(n+i,y);
        else{
            int pos=qmax(x,y);if(v[pos]>e[i].b){
                cut(e[pos-n].x,pos);cut(pos,e[pos-n].y);
                link(x,n+i);link(n+i,y);
            }
        }if(find(1)==find(n)) ans=min(ans,e[i].a+v[qmax(1,n)]);
    }if(ans==inf) puts("-1");
    else printf("%d\n",ans);
    return 0;
}