[bzoj3669][LCT]魔法森林

Description

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

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

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

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

Input

第1行包含两个整数N,M，表示无向图共有N个节点，M条边。
接下来M行，第行包含4个正整数Xi,Yi,Ai,Bi，描述第i条无向边。其中Xi与Yi为该边两个端点的标号，Ai与Bi的含义如题所述。
注意数据中可能包含重边与自环。

Output

输出一行一个整数：如果小E可以成功拜访到隐士，输出小E最少需要携带的守护精灵的总个数；如果无论如何小E都无法拜访到隐士，输出“-1”（不含引号）。

Sample Input

【输入样例1】

4 5

1 2 19 1

2 3 8 12

2 4 12 15

1 3 17 8

3 4 1 17

【输入样例2】

3 1

1 2 1 1

Sample Output

【输出样例1】

32

【样例说明1】

如果小E走路径1→2→4，需要携带19+15=34个守护精灵；

如果小E走路径1→3→4，需要携带17+17=34个守护精灵；

如果小E走路径1→2→3→4，需要携带19+17=36个守护精灵；

如果小E走路径1→3→2→4，需要携带17+15=32个守护精灵。

综上所述，小E最少需要携带32个守护精灵。

【输出样例2】

-1

【样例说明2】

小E无法从1号节点到达3号节点，故输出-1。

HINT

2<=n<=50,000

0<=m<=100,000

1<=ai ,bi<=50,000

题解

LCT好题！！！
论省选前一天发现模板大问题的奇妙性
你先把边按a排序，然后LCT维护b的MST，看到1和n连通了就统计一下答案
然后。我就发现了我的模板一个大问题
我在Access的时候并没有更新每个点的答案
这就有个很大的问题了
如果x和y本来就在一个splay里，那么y旋上去时会删除右孩子
我的代码是只在splay的时候更新的，于是下一步统计答案的时候y的右孩子的影响并没有消除！
所以。。
一定要在access的时候也更新答案呀

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
int bel[110000];
struct edge{int x,y,a,b;}w[110000];
bool cmp(edge n1,edge n2){return n1.a<n2.a;}
int s[210000];
struct node
{
    int f,c,maxpos,son[2];
    bool fz;node(){fz=false;}
}tr[211000];int n,m,cnt;
void rev(int now)
{
    int lc=tr[now].son[0],rc=tr[now].son[1];
    swap(tr[now].son[0],tr[now].son[1]);
    tr[lc].fz^=1;tr[rc].fz^=1;
    tr[now].fz=false;
}
void upd(int now)
{
    tr[now].maxpos=now;
    int lc=tr[now].son[0],rc=tr[now].son[1];
    if(s[tr[lc].maxpos]>s[tr[now].maxpos])tr[now].maxpos=tr[lc].maxpos;
    if(s[tr[rc].maxpos]>s[tr[now].maxpos])tr[now].maxpos=tr[rc].maxpos;
}
void rotate(int x,int w)
{
    int f=tr[x].f,ff=tr[f].f;
    int R,r;
    R=f;r=tr[x].son[w];
    tr[R].son[1-w]=r;
    if(r!=0)tr[r].f=R;
    R=ff;r=x;
    if(tr[R].son[0]==f)tr[R].son[0]=r;
    else if(tr[R].son[1]==f)tr[R].son[1]=r;
    tr[r].f=R;
    R=x;r=f;
    tr[R].son[w]=r;
    tr[r].f=R;
    upd(f);upd(x);
}
int tmp[211000];
void splay(int x,int rt)
{
    int i=x,s=0;
    while(tr[i].f!=rt && (tr[tr[i].f].son[0]==i || tr[tr[i].f].son[1]==i))
    {
        tmp[++s]=i;
        i=tr[i].f;
    }
    tmp[++s]=i;
    while(s)
    {
        i=tmp[s--];
        if(tr[i].fz)rev(i);
    }
    while(tr[x].f!=rt && (tr[tr[x].f].son[0]==x || tr[tr[x].f].son[1]==x))
    {
        int f=tr[x].f,ff=tr[f].f;
        if(ff==rt || (tr[ff].son[0]!=f && tr[ff].son[1]!=f))
        {
            if(tr[f].son[0]==x)rotate(x,1);
            else rotate(x,0);
        }
        else
        {
            if(tr[ff].son[0]==f && tr[f].son[0]==x)rotate(f,1),rotate(x,1);
            else if(tr[ff].son[1]==f && tr[f].son[0]==x)rotate(x,1),rotate(x,0);
            else if(tr[ff].son[1]==f && tr[f].son[1]==x)rotate(f,0),rotate(x,0);
            else rotate(x,0),rotate(x,1);
        }
    }
}
void access(int x)
{
    int y=0;
    while(x!=0)
    {
        splay(x,0);
        tr[x].son[1]=y;
        if(y!=0)tr[y].f=x;
        upd(x);
        y=x;x=tr[x].f;
    }
}
void markroot(int x){access(x);splay(x,0);tr[x].fz^=1;}
void link(int x,int y){markroot(x);tr[x].f=y;access(x);}
void cut(int x,int y)
{
    markroot(x);access(y);splay(y,0);
    tr[tr[y].son[0]].f=0;tr[y].son[0]=0;
}
int fa[51000];
int findfa(int x){return fa[x]!=x?fa[x]=findfa(fa[x]):fa[x];}
int sol(int x,int y)
{
    markroot(x);access(y);splay(y,0);
    return tr[y].maxpos;
}
int main()
{
    //freopen("forest4.in","r",stdin);
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)scanf("%d%d%d%d",&w[i].x,&w[i].y,&w[i].a,&w[i].b);
    sort(w+1,w+1+m,cmp);
    for(int i=1;i<=n;i++)fa[i]=i;cnt=n;
    int ans=999999999;
    for(int i=1;i<=m;i++)
    {
        s[i+n]=w[i].b;
        int p=findfa(w[i].x),q=findfa(w[i].y);
        if(w[i].x==w[i].y)continue;
        int x=w[i].x,y=w[i].y;
        if(p==q)
        {
            int u=sol(x,y);
            if(s[u]<=w[i].b)continue;
            cut(u,w[u-n].x);cut(u,w[u-n].y);
        }
        else fa[p]=q;
        link(x,i+n);link(i+n,y);
        if(findfa(1)==findfa(n))ans=min(ans,w[i].a+s[sol(1,n)]);
    }
    if(ans==999999999)printf("-1\n");
    else printf("%d\n",ans);
    return 0;
}