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;
}