[BZOJ2959]长跑(LCT+并查集)

题目:

我是超链接

题解:

诶这个题连成一个环是不是还要tarjan缩个点啥的啊,加边的时候会使原本构不成环的部分构成新的环所以考虑在树上用LCT进行动态的缩点和维护操作。

如果加入某条边之前两个端点就联通,那么我们就将两个点之间的路径提出来,加入改边会使这条路径上的点可以互相到达,所以我们把这些点缩成一个代表节点,然后将权值都加给这个节点。并用另外一个并查集(find-fa)维护每个点所属的代表元素信息。

对于询问2,其实就变成了从LCT中提取缩点后两点所属的点之间的路径上的权值总和。

写LCT是需要注意,以前找的是f[x],现在应该是找find(f[x]),就是f[x]所属的缩点的代表元素。

这样总体来说我们要开两个并查集,一个用来维护缩点的信息(find-fa),另一个维护两条边是否相连(finds-fs)
然后就T了,常数怎么这么大,能不find就不要find,多用变量代替

代码:

#include <cstdio>
#include <iostream>
using namespace std;
const int N=150005;
int f[N],fa[N],fs[N],ch[N][2],size[N],sum[N],a[N],delta[N],stack[N];
int read()
{
    char ch=getchar(); int x=0;
    while (ch>'9' || ch<'0') ch=getchar();
    while (ch<='9' && ch>='0') x=x*10+ch-'0',ch=getchar();
    return x;
}
int find(int x)
{
    if (fa[x]!=x) fa[x]=find(fa[x]);
    return fa[x];
}
int finds(int x)
{
    if (fs[x]!=x) fs[x]=finds(fs[x]);
    return fs[x];
}
void updata(int now)
{
    sum[now]=size[now];
    if (ch[now][0]) sum[now]+=sum[ch[now][0]];
    if (ch[now][1]) sum[now]+=sum[ch[now][1]];
}
bool isroot(int x){int sb=find(f[x]); return ch[sb][0]!=x && ch[sb][1]!=x;}
int get(int x){return ch[find(f[x])][1]==x;}
void rotate(int x)
{
    int old=find(f[x]),oldf=find(f[old]),which=get(x);
    f[x]=oldf; if (!isroot(old)) ch[oldf][get(old)]=x;
    ch[old][which]=ch[x][which^1]; f[ch[x][which^1]]=old;
    ch[x][which^1]=old; f[old]=x;
    updata(old); updata(x);
}
void pushdown(int now)
{
    if (delta[now])
    {
        delta[ch[now][0]]^=1; delta[ch[now][1]]^=1;
        swap(ch[now][0],ch[now][1]);
        delta[now]=0;
    }
}
void splay(int x)
{
    int top=0,i;
    for (i=x;!isroot(i);i=find(f[i])) stack[++top]=i;
    stack[++top]=i;
    for (int i=top;i>=1;i--) pushdown(stack[i]);
    for (;!isroot(x);rotate(x))
    {
        int sb=find(f[x]);
        if (!isroot(sb)) rotate(get(x)==get(sb)?sb:x);
    }
}
void access(int x)
{
    int t=0;
    for (;x;t=x,x=find(f[x]))
    {
        splay(x);
        ch[x][1]=t;
        updata(x);
    }
}
void dfs(int x,int go)
{
    x=find(x);
    if (ch[x][0]) dfs(ch[x][0],go);
    if (ch[x][1]) dfs(ch[x][1],go);
    if (x!=go){fa[x]=go; size[go]+=size[x];}
}
void reverse(int x){access(x);splay(x);delta[x]^=1;}
void un(int x,int y){reverse(x); access(y); splay(y); dfs(y,y); updata(y);}
void Link(int x,int y){reverse(x); f[x]=y; splay(x);}
int qurry(int x,int y){reverse(x); access(y); splay(y);return sum[y];}
int main()
{
    int n,m;n=read();m=read();
    for (int i=1;i<=n;i++) fa[i]=i,fs[i]=i;
    for (int i=1;i<=n;i++) size[i]=read(),a[i]=size[i];
    for (int i=1;i<=m;i++)
    {
        int id,x,y,xx;id=read();x=read();y=read();xx=x;
        x=find(x); 
        if (id==1)
        {
            y=find(y);if (x==y) continue;
            int f1=finds(x),f2=finds(y);
            if (f1==f2) un(x,y);
            else Link(x,y),fs[f1]=f2;
        }else if (id==2)
        {
            access(x); splay(x);
            size[x]=size[x]+y-a[xx];a[xx]=y;
            updata(x);
        }else
        {
            y=find(y);
            if (finds(x)!=finds(y)) printf("-1\n");else printf("%d\n",qurry(x,y));
        }
    }
}
  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值