NOI2014 洛谷2387 bzoj3669 魔法森林 LCT

题目描述：

为了得到书法大家的真传，小 E 同学下定决心去拜访住在魔法森林中的隐 士。魔法森林可以被看成一个包含 n 个节点 m 条边的无向图，节点标号为 1,2,3,…,n，边标号为 1,2,3,…,m。初始时小 E 同学在 1 号节点，隐士则住在 n 号节点。小 E 需要通过这一片魔法森林，才能够拜访到隐士。

魔法森林中居住了一些妖怪。每当有人经过一条边的时候，这条边上的妖怪 就会对其发起攻击。幸运的是，在 1 号节点住着两种守护精灵：A 型守护精灵与 B 型守护精灵。小 E 可以借助它们的力量，达到自己的目的。

只要小 E 带上足够多的守护精灵，妖怪们就不会发起攻击了。具体来说，无 向图中的每一条边 ei 包含两个权值 ai 与 bi 。若身上携带的 A 型守护精灵个数不 少于 ai ，且 B 型守护精灵个数不少于 bi ，这条边上的妖怪就不会对通过这条边 的人发起攻击。当且仅当通过这片魔法森林的过程中没有任意一条边的妖怪向 小 E 发起攻击，他才能成功找到隐士。

由于携带守护精灵是一件非常麻烦的事，小 E 想要知道，要能够成功拜访到 隐士，最少需要携带守护精灵的总个数。守护精灵的总个数为 A 型守护精灵的 个数与 B 型守护精灵的个数之和。

题解：

用LCT维护动态最小生成树。首先我们按每条边的A精灵排序，用LCT维护需要多少B精灵。我们需要判断每条边的两端的点是否已经连通，维护连通性应该写并查集会速度更快，我写的是LCT维护连通性。我们如果两端点x,y不连通，那么我们把它们link一下，连通的话就判断当前边的B精灵是否比原来生成树上x,y之间的最大的B精灵数量小，如果大的话什么都不干，如果小的话把原来那条大的边cut掉，然后link这条新的边。

特别来说明一下，这里的link和cut操作与一般的lct写法并不一样。cut不是很难理解，主要讲讲link。

lct上的点表示的含义和一般的lct有所不同，有的表示图上的边，有的表示图上的点，link时分别把树上这条边的点与表示这条边两端的点的点相连。我为了区别树上表示边的点和表示点的点，把表示点的点x进行了一个x|p的操作。

 

之后我们每处理完一条边，就判断一下当前1号点和n号点的连通性，如果连通就更新答案。

 

具体操作见代码吧。

 

#include <bits/stdc++.h>
using namespace std;

int n,m,f[200001],c[200001][2],mx[200001],rev[200001],st[200001],ans=999999999;
int p=131072;
struct node
{
    int u,v,x,y;
}e[400001];
int cmp(node a,node b)
{
    return a.x<b.x;
}
int nroot(int x)
{
    return c[f[x]][1]==x||c[f[x]][0]==x;
}
void pushup(int x)
{
    mx[x]=x;
    if(e[mx[x]].y<e[mx[c[x][0]]].y)
    mx[x]=mx[c[x][0]];
    if(e[mx[x]].y<e[mx[c[x][1]]].y)
    mx[x]=mx[c[x][1]];
}
void pushdown(int x)
{
    if(rev[x])
    {
        swap(c[x][0],c[x][1]);
        rev[c[x][0]]^=1;
        rev[c[x][1]]^=1;
        rev[x]=0;
    }
}
void rotate(int x)
{
    int y=f[x],z=f[y],k=c[y][1]==x,w=c[x][!k];
    if(nroot(y))
    c[z][c[z][1]==y]=x;
    c[x][!k]=y;
    c[y][k]=w;
    if(w)
    f[w]=y;
    f[x]=z;
    f[y]=x;
    pushup(y);
} 
void splay(int x)
{
    int y=x,z=0;
    st[++z]=y;
    while(nroot(y))
    {
        y=f[y];
        st[++z]=y;
    }
    while(z)
    pushdown(st[z--]);
    while(nroot(x))
    {
        y=f[x],z=f[y];
        if(nroot(y))
        {
            if(c[z][0]==y ^ c[y][0]==x)
            rotate(y);
            else
            rotate(x);
        }
        rotate(x);
    }
    pushup(x);
}
void access(int x)
{
    int y=0;
    while(x!=0)
    {
        splay(x);
        c[x][1]=y;
        pushup(x);
        if(y!=0)
        f[y]=x;
        y=x;
        x=f[x];
    }
}
void makeroot(int x)
{
    access(x);
    splay(x);
    rev[x]^=1;
}
int findroot(int x)
{
    access(x);
    splay(x);
    while(c[x][0])
    x=c[x][0];
    return x;
}
void link(int x)
{
    int y=e[x].u,z=e[x].v;
    makeroot(z);
    f[z]=x;
    f[x]=y;
}
void cut(int x)
{
    access(e[x].v);
    splay(x);
    c[x][0]=c[x][1]=f[c[x][0]]=f[c[x][1]]=0;
    pushup(x);
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d%d%d",&e[i].u,&e[i].v,&e[i].x,&e[i].y);
        e[i].u|=p;
        e[i].v|=p;
    }
    sort(e+1,e+m+1,cmp);
    for(int i=1;i<=m;i++)
    {
        int y,z;
        y=e[i].u,z=e[i].v;
        if(y==z)
        continue;
        makeroot(y);
        if(y!=findroot(z))
        link(i);
        else if(e[i].y<e[mx[z]].y)
        {
            cut(mx[z]);
            link(i);
        }
        makeroot(1|p);
        if((1|p)==findroot(n|p))
        ans=min(ans,e[i].x+e[mx[n|p]].y);
    }
    if(ans==999999999)//始终没法连通就输出-1
    printf("-1\n");
    else
    printf("%d\n",ans);
    return 0;
}