【并查集+LCT维护生成树】BZOJ3669(Noi2014)[魔法森林]题解

题目描述

有n个点和m条边，每条边有a,b两个权值。求一条1->n的路径，使得max(a)+max(b)最小。

解题报告

ps：调了两天终于调出来了……
这道题是经典的并查集+LCT维护生成树（好像最短路搞好点可以水过？）。我们先把边按照a从小到大排序，这样每次枚举的边就是当前最大边。每当新来一条边，我们先用并查集检查两个点是否已经相连，如果没有，那么直接加入。否则如果当前边的b小于两个点之间b的最大值，就干掉原先的最大边，换成当前边。
然而我感觉完全不会写，之后才知道这道题是边权，所以LCT需要采用这种连边方式（E是边，x和y是E的两个节点）：x->E,E->y。掌握了新姿势啊！

示例程序

#include<cstdio>
#include<algorithm>
using namespace std;
const int maxn=50000,maxm=100000,maxt=maxn+maxm,MAXINT=((1<<30)-1)*2+1;

int n,m,ans,father[maxn+5];
struct Edge
{
    int x,y,a,b;
    bool operator < (const Edge &c) const {return a<c.a;}
};
Edge e[maxm+5];
struct LCT //LCT结构体
{
    int son[maxt+5][2],father[maxt+5],val[maxt+5],ID[maxt+5];
    bool flip[maxt+5];
    void clear()
    {
        for (int i=0;i<=n+m;i++)
        {
            val[i]=-MAXINT;if (i>n) val[i]=e[i-n].b;ID[i]=i;
            son[i][0]=son[i][1]=father[i]=0;
            flip[i]=false;
        }
    }
    bool is_ro(int p) {return p!=son[father[p]][0]&&p!=son[father[p]][1];}
    int wt_son(int p) {return p==son[father[p]][1];}
    void Pushup(int p)
    {
        ID[p]=p;
        for (int i=0;i<=1;i++)
            if (val[ID[son[p][i]]]>val[ID[p]]) ID[p]=ID[son[p][i]];
    }
    void Add_flip(int p) {swap(son[p][0],son[p][1]);flip[p]^=1;}
    void Pushdown(int p)
    {
        if (!flip[p]) return;
        flip[p]^=1;Add_flip(son[p][0]);Add_flip(son[p][1]);
    }
    void Rotate(int p)
    {
        int d=wt_son(p),fa=father[p];
        if (!is_ro(fa)) son[father[fa]][wt_son(fa)]=p;
        son[fa][d]=son[p][d^1];father[son[p][d^1]]=fa;son[p][d^1]=fa;
        father[p]=father[fa];father[fa]=p;
        Pushup(fa);Pushup(p);
    }
    int top,stk[maxn+5];
    void Splay(int p)
    {
        top=0;stk[++top]=p;
        for (int i=p;!is_ro(i);i=father[i]) stk[++top]=father[i];
        while (top) Pushdown(stk[top--]);
        while (!is_ro(p))
        {
            int fa=father[p],d1=wt_son(p);
            if (!is_ro(fa))
            {
                int d2=wt_son(fa);
                if (d1==d2) Rotate(fa); else Rotate(p);
            }
            Rotate(p);
        }
    }
    void Access(int p)
    {
        int lst=0;
        while (p)
        {
            Splay(p);son[p][1]=lst;Pushup(p);
            lst=p;p=father[p];
        }
    }
    void make_ro(int p) {Access(p);Splay(p);Add_flip(p);}
    void Link(int x,int y) {make_ro(x);father[x]=y;}
    void getpath(int x,int y) {make_ro(y);Access(x);Splay(y);}
    int getMAX(int x,int y) {getpath(x,y);return val[ID[y]];}
    void Cut(int x,int y)
    {
        make_ro(x);Access(y);Splay(y);
        son[y][0]=father[x]=0;Pushup(y);
    }
};
LCT tr;

char readc()
{
    static char buf[100000],*l=buf,*r=buf;
    if (l==r) r=(l=buf)+fread(buf,1,100000,stdin);
    if (l==r) return EOF; else return *l++;
}
bool Eoln(char ch) {return ch==10||ch==13||ch==EOF;}
int readi(int &x)
{
    int tot=0,f=1;char ch=readc(),lst=ch;
    while ('9'<ch||ch<'0') {if (ch==EOF) return EOF;lst=ch;ch=readc();}
    if (lst=='-') f=-f;
    while ('0'<=ch&&ch<='9') tot=tot*10+ch-48,ch=readc();
    x=tot*f;
    return Eoln(ch);
}
int getfa(int x)
{
    if (father[x]==x) return x;
    return father[x]=getfa(father[x]);
}
void Merge(int i) //将i边加入图
{
    int x=e[i].x,y=e[i].y,fx=getfa(x),fy=getfa(y);
    if (fx!=fy) {father[fx]=fy;tr.Link(x,n+i);tr.Link(n+i,y);return;}
    if (tr.getMAX(x,y)>e[i].b) //换边
    {
        int who=tr.ID[y];
        tr.Cut(e[who-n].x,who);tr.Cut(e[who-n].y,who);
        tr.Link(x,n+i);tr.Link(n+i,y);
    }
}
int main()
{
    freopen("program.in","r",stdin);
    freopen("program.out","w",stdout);
    readi(n);readi(m);
    for (int i=1;i<=m;i++)
        readi(e[i].x),readi(e[i].y),readi(e[i].a),readi(e[i].b);
    sort(e+1,e+1+m);tr.clear();
    for (int i=1;i<=n;i++) father[i]=i;
    ans=MAXINT;
    for (int i=1;i<=m;i++)
    {
        Merge(i);if (getfa(1)!=getfa(n)) continue;
        ans=min(ans,e[i].a+tr.getMAX(1,n));
    }
    if (ans==MAXINT) ans=-1;printf("%d\n",ans);
    return 0;
}