[WC2006]水管局长

原创于 2020-01-12 12:54:00 发布 · 212 阅读

0 ·

CC 4.0 BY-SA版权

最小生成树同时被 2 个专栏收录

10 篇文章

订阅专栏

link-cut-tree

8 篇文章

订阅专栏

本文介绍了一种算法，用于动态维护最小生成树，支持在图中进行边的删除和查询操作，通过逆序处理边的删除转换为边的添加，利用lctlctlct数据结构实现最小生成树的高效更新。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

一、题目

点此看题

二、解法

简述一下题意吧，本题就是给你 $n$ 个点 $m$ 条边，支持下面两种操作：

1、询问 $(x, y)$ 所有路径上最大值的最小值。
2、去掉边 $(x, y)$ ，保证边在给定的 $m$ 条边中，并且图任何时候联通。

考虑询问，如果我们能维护一棵最小生成树，那么 $(x, y)$ 在最小生成树上路径的最大值即是答案，但是本题需要删边，我们可以把删边逆序，就变成了加边，我们要做的就是动态维护这颗最小生成树。

这里我们使用 $l c t$ ，我们重新写一个最小生成树加边函数，如果边的两端不连通，直接加入；否则找到这条边构成的环上的最大值，和新加的这一条边作比较，如果新加的这条边大则不操作，否则我们删除环上的最大边之后再加上这条边， $l c t$ 中维护最大边权和最大边权的编号即可。

还有一个问题，就是 $l c t$ 处理不了边，我们可以把边看成点，连接 $(x, y)$ 就代表连接 $(x, i d), (i d, y)$ ， $i d$ 是边的编号，我们把边权放在这个点上。

据说本题还有一个数据加强版，蒟蒻的代码写的太丑了，过不了，代码仅供参考。

#include <cstdio>
#include <iostream>
#include <map>
using namespace std;
const int Q = 100005;
const int M = 110005;
int read()
{
    int x=0,flag=1;
    char c;
    while((c=getchar())<'0' || c>'9') if(c=='-') flag=-1;
    while(c>='0' && c<='9') x=(x<<3)+(x<<1)+(c^48),c=getchar();
    return x*flag;
}
int n,m,q,a[M],b[M],c[M],op[Q],x[Q],y[Q],ans[Q];
int par[M],ch[M][2],val[M],Max[M],pos[M],fl[M],st[M];
map<int,int> mp,id;
int nrt(int x)
{
    return ch[par[x]][0]==x || ch[par[x]][1]==x;
}
int Hash(int a,int b)
{
    if(a>=b) swap(a,b);
    return 1ll*a*n+b;
}
int chk(int x)
{
    return ch[par[x]][1]==x;
}
void push_up(int x)
{
    if(!x) return ;
    Max[x]=max(val[x],max(Max[ch[x][0]],Max[ch[x][1]]));
    if(Max[x]==val[x]) pos[x]=x;
    else if(Max[x]==Max[ch[x][0]]) pos[x]=pos[ch[x][0]];
    else if(Max[x]==Max[ch[x][1]]) pos[x]=pos[ch[x][1]];
}
void flip(int x)
{
    if(!x) return ;
    swap(ch[x][0],ch[x][1]);
    fl[x]^=1;
}
void push_down(int x)
{
    if(!x) return ;
    if(fl[x])
    {
        flip(ch[x][0]);
        flip(ch[x][1]);
        fl[x]=0;
    }
}
void rotate(int x)
{
    int y=par[x],z=par[y],k=chk(x),w=ch[x][k^1];
    ch[y][k]=w;
    par[w]=y;
    if(nrt(y)) ch[z][chk(y)]=x;
    par[x]=z;
    ch[x][k^1]=y;
    par[y]=x;
    push_up(y);
    push_up(x);
}
void splay(int x)
{
    int y=x,z=0;
    st[++z]=y;
    while(nrt(y)) st[++z]=y=par[y];
    while(z) push_down(st[z--]);
    while(nrt(x))
    {
        int y=par[x],z=par[y];
        if(nrt(y))
        {
            if(chk(y)==chk(x)) rotate(y);
            else rotate(x);
        }
        rotate(x);
    }
}
void access(int x)
{
    for(int y=0; x; x=par[y=x])
        splay(x),ch[x][1]=y,push_up(x);
}
void makeroot(int x)
{
    access(x);
    splay(x);
    flip(x);
}
int findroot(int x)
{
    access(x);
    splay(x);
    while(ch[x][0]) push_down(x),x=ch[x][0];
    splay(x);
    return x;
}
void split(int x,int y)
{
    makeroot(x);
    access(y);
    splay(y);
}
void link(int x,int y)
{
    makeroot(x);
    if(findroot(y)!=x) par[x]=y;
}
void cut(int x,int y)
{
    makeroot(x);
    if(findroot(y)==x && par[y]==x && !ch[y][0])
    {
        par[y]=ch[x][1]=0;
        push_up(x);
    }
}
void add(int i)
{
    val[i+n]=c[i];
    push_up(i+n);
    if(findroot(a[i])^findroot(b[i]))
    {
        link(a[i],i+n);
        link(i+n,b[i]);
        return ;
    }
    split(a[i],b[i]);
    if(Max[b[i]]<=c[i]) return ;
    int t=pos[b[i]];
    cut(a[t-n],pos[t]);
    cut(pos[t],b[t-n]);
    link(a[i],i+n);
    link(i+n,b[i]);
}
int main()
{
    Max[0]=-(1<<30);
    n=read();
    m=read();
    q=read();
    for(int i=1; i<=m; i++)
    {
        a[i]=read();
        b[i]=read();
        c[i]=read();
        id[Hash(a[i],b[i])]=i;
    }
    for(int i=1; i<=q; i++)
    {
        op[i]=read();
        x[i]=read();
        y[i]=read();
        if(op[i]==2)
            mp[Hash(x[i],y[i])]=1;
    }
    for(int i=1; i<=m; i++)
        if(!mp[Hash(a[i],b[i])])
            add(i);
    for(int i=q; i>0; i--)
    {
        if(op[i]==1)
        {
            split(x[i],y[i]);
            ans[i]=Max[y[i]];
        }
        else
            add(id[Hash(x[i],y[i])]);
    }
    for(int i=1; i<=q; i++)
        if(op[i]==1)
            printf("%d\n",ans[i]);
}