[bzoj 3669--NOI2014]魔法森林

爲了得到書法大家的真傳,小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型守護精靈的個數之和。

這道題可以算是一道LCT的經典題目。我們從簡單開始考慮,如果只有一維,我們就可以用最小生成樹解決。那二維的基本思想其實也沒多大變化,我們先按第一維排序一下,使第一維影響消去,那第二維用LCT維護之。不斷加邊,如果當前兩點已聯通,則找出當前兩點之間的最大值,如果比當前的b值要大,則刪掉之前聯通兩點之間的邊,加上當前的邊,否則就不管當前的邊。那這道題就解決了。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
using namespace std;
struct node
{
    int f,sd,son[2];
    bool lazy;
    node(){f=son[0]=son[1]=0;lazy=false;}
}tr[201000];
int a[201000];
void update(int x)
{
    int lc=tr[x].son[0],rc=tr[x].son[1];
    tr[x].sd=x;
    if(a[tr[lc].sd]>a[tr[x].sd])tr[x].sd=tr[lc].sd;
    if(a[tr[rc].sd]>a[tr[x].sd])tr[x].sd=tr[rc].sd;
}
void reverse(int x)
{
    tr[x].lazy=false;
    swap(tr[x].son[0],tr[x].son[1]);
    int lc=tr[x].son[0],rc=tr[x].son[1];
    tr[lc].lazy=1-tr[lc].lazy;
    tr[rc].lazy=1-tr[rc].lazy;
}
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;
    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;
    update(f);update(x);
}
int tmp[201000];
void splay(int x,int rt)
{
    int s=x,ss=0;
    while(tr[s].f!=rt && (tr[tr[s].f].son[0]==s || tr[tr[s].f].son[1]==s))
    {
        tmp[++ss]=s;
        s=tr[s].f;
    }
    tmp[++ss]=s;
    while(ss)
    {
        s=tmp[ss--];
        if(tr[s].lazy==true)reverse(s);
    }
    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[f].son[0]==x && tr[ff].son[0]==f)rotate(f,1),rotate(x,1);
            if(tr[f].son[1]==x && tr[ff].son[1]==f)rotate(f,0),rotate(x,0);
            if(tr[f].son[0]==x && tr[ff].son[1]==f)rotate(x,1),rotate(x,0);
            if(tr[f].son[1]==x && tr[ff].son[0]==f)rotate(x,0),rotate(x,1);
        }
    }
}
void access(int x)
{
    int y=0;
    while(x)
    {
        splay(x,0);
        tr[x].son[1]=y;
        if(y!=0)tr[y].f=x;
        y=x;x=tr[x].f;
    }
}
void makeroot(int x)
{
    access(x);splay(x,0);
    tr[x].lazy=1-tr[x].lazy;
}
void Link(int x,int y)
{
    makeroot(x);
    tr[x].f=y;
    access(x);
}
void Cut(int x,int y)
{
    makeroot(x);
    access(y);splay(y,0);
    tr[tr[y].son[0]].f=0;tr[y].son[0]=0;
    update(y);
}
int solve(int x,int y)
{
    makeroot(x);
    access(y);splay(y,0);
    update(y);
    return tr[y].sd;
}
int fa[51000];
int findfa(int x)
{
    if(fa[x]==x)return x;
    else
    {
        fa[x]=findfa(fa[x]);
        return fa[x];
    }
}
struct edge
{
    int x,y,w1,w2;
}e[110000];
bool cmp(edge a,edge b)
{
    if(a.w1>b.w1)return false;
    if(a.w1<b.w1)return true;
    return 0;
}
int main()
{
    int n,m,ans=999999999;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)fa[i]=i;
    for(int i=1;i<=m;i++)scanf("%d%d%d%d",&e[i].x,&e[i].y,&e[i].w1,&e[i].w2);
    sort(e+1,e+m+1,cmp);
    for(int i=1;i<=m;i++)
    {
        int x=e[i].x,y=e[i].y;
        int fx=findfa(x),fy=findfa(y);
        if(fx==fy)
        {
            int wy=solve(x,y);
            if(a[wy]<=e[i].w2)continue;
            Cut(wy,e[wy-n].x);
            Cut(wy,e[wy-n].y);
        }
        else fa[fx]=fy;
        a[n+i]=e[i].w2;tr[n+i].sd=n+i;
        Link(x,n+i);Link(y,n+i);
        if(findfa(1)==findfa(n))ans=min(ans,e[i].w1+a[solve(1,n)]);
    }
    if(ans==999999999)printf("-1\n");
    else printf("%d\n",ans);
    return 0;
}
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章