bzoj4668 冷战 (并查集按秩合并)

83 篇文章 0 订阅

bzoj4668 冷战

原题地址http://www.lydsy.com/JudgeOnline/problem.php?id=4668

题意:
1946 年 3 月 5 日,英国前首相温斯顿·丘吉尔在美国富尔顿发表“铁幕演说”,正式拉开了冷战序幕。
美国和苏联同为世界上的“超级大国”,为了争夺世界霸权,两国及其盟国展开了数十年的斗争。在这段时期,虽然分歧和冲突严重,但双方都尽力避免世界范围的大规模战争(第三次世界大战)爆发,其对抗通常通过局部代理战争、科技和军备竞赛、太空竞争、外交竞争等“冷”方式进行,即“相互遏制,不动武力”,因此称之为“冷战”。
Reddington 是美国的海军上将。由于战争局势十分紧张,因此他需要时刻关注着苏联的各个活动,避免使自己的国家陷入困境。苏联在全球拥有 N 个军工厂,但由于规划不当,一开始这些军工厂之间是不存在铁路的,为了使武器制造更快,苏联决定修建若干条道路使得某些军工厂联通。Reddington 得到了苏联的修建日程表,并且他需要时刻关注着某两个军工
厂是否联通,以及最早在修建哪条道路时会联通。具体而言,现在总共有M 个操作,操作分为两类:
• 0 u v,这次操作苏联会修建一条连接 u 号军工厂及 v 号军工厂的铁路,注意铁路都是双向的;
• 1 u v, Reddington 需要知道 u 号军工厂及 v 号军工厂最早在加入第几条条铁路后会联通,假如到这次操作都没有联通,则输出 0;
作为美国最强科学家, Reddington 需要你帮忙设计一个程序,能满足他的要求。
Input
第一行两个整数 N, M。
接下来 M 行,每行为 0 u v 或 1 u v 的形式。
数据是经过加密的,对于每次加边或询问,真正的 u, v 都等于读入的u, v 异或上上一次询问的答案。一开始这个值为 0。
1 ≤ N, M ≤ 500000,解密后的 u, v 满足1 ≤ u, v ≤ N, u不等于v
Output
对于每次 1 操作,输出 u, v 最早在加入哪条边后会联通,若到这个操
作时还没联通,则输出 0。
Sample Input
5 9
0 1 4
1 2 5
0 2 4
0 3 4
1 3 1
0 7 0
0 6 1
0 1 6
1 2 6
Sample Output
0
3
5
数据范围
1 ≤ N, M ≤ 500000,解密后的 u, v 满足1 ≤ u, v ≤ N, u不等于v

题解:
并查集按秩合并,点权为其与父亲连通的时刻,每次查询两点路径上最大点权。
我们再跳lca时和一般求lca一样,先把dep跳到相同,再一步一步往上跳。这里无需纠结怎么维护每个点的dep,由于树高logn,我们可以先跳一边查dep。
之前考试题还给了一种思路:

与普通并查集不同,这道题我们如果要合并两棵子树x,y,我们不在x,y之间连边,而是新建一个点a,分别将ax,ay 连边。
在没有路径压缩的情况下,很容易知道,最终形成一棵树,根节点为最后一次操作的新建节点。
求x,y 的合并时间,可以直接查找x,y 在树上的lca。
当然,为了快速根节点,在之前的解法上加路径压缩不会影响正确性,而且还可以减少一个log
这道题还可以通过并查集+整体二分解决,但是这个解法不在讨论范围内。

(发现之前的按秩合并写错了,应维护树高而不是size。)

代码:

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
const int N=500005;
int n,m,fa[N],h[N],t[N],last=0,cnt=0;
int getfa(int x)
{
    if(x==fa[x]) return x;
    else return getfa(fa[x]);
}
void merge(int x,int y,int val)
{
    if(x==y) return;
    if(h[x]<h[y]) swap(x,y);
    fa[y]=x; t[y]=val;
    if(h[x]==h[y]) h[x]++;
}
int getdep(int x)
{
    if(fa[x]==x) return 0;
    else return getdep(fa[x])+1;
}
int query(int x,int y)
{
    if(getfa(x)!=getfa(y)) return 0;
    int ans=0;
    int depx=getdep(x); int depy=getdep(y); 
    if(depx<depy) {swap(x,y); swap(depx,depy);}
    while(depx!=depy)
    {
        ans=max(t[x],ans);
        x=fa[x]; depx--;
    }
    while(x!=y)
    {
        ans=max(ans,max(t[x],t[y]));
        x=fa[x]; y=fa[y];
    }
    return ans;
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) fa[i]=i;
    for(int i=1;i<=m;i++)
    {
        int opt,x,y;
        scanf("%d%d%d",&opt,&x,&y);
        x=x^last; y=y^last;
        if(opt==0)
        merge(getfa(x),getfa(y),++cnt);
        else
        {
            last=query(x,y);
            printf("%d\n",last);
        }
    }

    return 0;
}
  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值