BZOJ3669: [Noi2014]魔法森林

题目链接

学习了一下黄学长的做法：http://hzwer.com/3845.html
边集按a值从小到大排序，依次加边，若构成环则删去环上最大边。
LCT维护的是点，由于是动态，所以不方便像树链剖分一样直接将边转到点上，一个很好用的技巧就是将边建新点，然后就好做了。只有边建成的点有权值，权值为b[i]，LCT维护最大值的位置。当1和n连通时更新答案。（还是比较好实现的吧。。打完莫名其妙就胡对了）

（学长写的加了各种奇奇怪怪优化的最短路跑出了1100ms+。。我的LCT直接跑到了5700ms+。。）

【代码】

#include <cstdio>
#include <iostream>
#include <queue>
#include <vector>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <stack>
#define N 150005 
#define M 2001005
#define INF 1e9
#define mod 51061
using namespace std;
typedef unsigned int ll;
typedef pair<ll,ll> pa;

int read()
{
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
    return x*f;
}

int n,m,top,ans=INF;
int son[N][2],st[N],fa[N],mx[N],a[N];
bool rev[N];

class Edge{
    public:
        int x,y,A,B;
}e[100005];

bool cmp(Edge a,Edge b){
    return a.A<b.A;
}

void Input_Init()
{
    n=read();m=read();a[0]=-INF;
    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+1+m,cmp);
}

bool isroot(int x){
    return son[fa[x]][0]!=x&&son[fa[x]][1]!=x;
}

void pushup(int k)
{
    int l=son[k][0],r=son[k][1];
    mx[k]=k;
    if(a[mx[k]]<a[mx[l]]) mx[k]=mx[l];
    if(a[mx[k]]<a[mx[r]]) mx[k]=mx[r];
}

void pushdown(int k)
{
    if(!rev[k]) return;
    rev[k]=0,rev[son[k][0]]^=1,rev[son[k][1]]^=1;
    swap(son[k][1],son[k][0]);
}

void Rotate(int x)
{
    int y=fa[x],z=fa[y],l,r;
    l=son[y][1]==x;r=l^1;
    if(!isroot(y)) son[z][son[z][1]==y]=x;
    fa[x]=z;fa[y]=x;fa[son[x][r]]=y;
    son[y][l]=son[x][r];son[x][r]=y;
    pushup(y);pushup(x);
}

void Splay(int x)
{
    st[++top]=x;
    for(int i=x;!isroot(i);i=fa[i]) st[++top]=fa[i];
    while(top) pushdown(st[top--]);
    while(!isroot(x))
    {
        int y=fa[x],z=fa[y];
        if(!isroot(y)) {
            if(son[z][0]==y^son[y][0]==x) Rotate(x);
            else Rotate(y);
        }
        Rotate(x);
    }
}

void Access(int x){
    for(int t=0;x;t=x,x=fa[x])
        Splay(x),son[x][1]=t,pushup(x);
}

void Make_Root(int x){
    Access(x);Splay(x);rev[x]^=1;
}

int Find(int x)
{
    Access(x);Splay(x);
    while(son[x][0]) x=son[x][0];
    return x; 
}

void Link(int x,int y){
    Make_Root(x);fa[x]=y;
}

void Cut(int x,int y){
    Make_Root(x);Access(y);Splay(y);
    son[y][0]=fa[x]=0;pushup(y);
}

int Query(int x,int y){
    Make_Root(x);Access(y);Splay(y);
    return mx[y];
}

void Solve()
{
    for(int i=1;i<=m;i++)
    {
        static int x,y,A,B;
        x=e[i].x,y=e[i].y,A=e[i].A,B=e[i].B;
        if(Find(x)==Find(y))
        {
            int t=Query(x,y);
            if(a[t]>B)
            {
                Cut(t,e[t-n].x);
                Cut(t,e[t-n].y);
            }
            else {
                if(Find(1)==Find(n))
                    ans=min(ans,A+a[Query(1,n)]);
                continue;
            }
        }
        a[n+i]=B;mx[n+i]=n+i;
        Link(x,i+n);Link(y,i+n);
        if(Find(1)==Find(n)) ans=min(ans,A+a[Query(1,n)]);
    }
    if(ans==INF) printf("-1\n");
    else printf("%d\n",ans);
}

int main()
{
    Input_Init();
    Solve(); 
    return 0;
}