3674: 可持久化并查集加强版
Time Limit: 15 Sec Memory Limit: 256 MBSubmit: 3289 Solved: 1229
[Submit][Status][Discuss]
Description
Description:
自从zkysb出了可持久化并查集后……
hzwer:乱写能AC,暴力踩标程
KuribohG:我不路径压缩就过了!
ndsf:暴力就可以轻松虐!
zky:……
n个集合 m个操作
操作:
1 a b 合并a,b所在集合
2 k 回到第k次操作之后的状态(查询算作操作)
3 a b 询问a,b是否属于同一集合,是则输出1否则输出0
请注意本题采用强制在线,所给的a,b,k均经过加密,加密方法为x = x xor lastans,lastans的初始值为0
0<n,m<=2*10^5
Sample Input
5 6
1 1 2
3 1 2
2 1
3 0 3
2 1
3 1 2
1 1 2
3 1 2
2 1
3 0 3
2 1
3 1 2
Sample Output
1
0
1
0
1
可持久化并查集=可持久化数组=可持久化线段树
可持久化线段树维护当前状态下集合的关系和秩
秩:以该元素为代表元的所有元素中的最大深度
对于并查集,就算不压缩路径
按照秩去启发式合并
效率也能保持在O(logn)以内(别问我为啥)
#include<cmath>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<complex>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<string>
#include<bitset>
#include<queue>
#include<map>
#include<set>
using namespace std;
inline int read()
{
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch<='9'&&ch>='0'){x=10*x+ch-'0';ch=getchar();}
return x*f;
}
inline void print(int x)
{if(x<0)putchar('-'),x=-x;if(x>=10)print(x/10);putchar(x%10+'0');}
const int N=200100,M=5000100;
struct seg_tree{int ls,rs,fa,d;}tr[M];
int n,m,root[N],ans;
int cnt=0;
void build(int &k,int l,int r)
{
k=++cnt;
if(l==r){tr[k].fa=l;tr[k].d=1;return ;}
int mid=l+r>>1;
build(tr[k].ls,l,mid);build(tr[k].rs,mid+1,r);
}
inline int findpos(int rt,int x)
{
int t=root[rt],l=1,r=n;
while(l<r)
{
int mid=l+r>>1;
if(x<=mid)r=mid,t=tr[t].ls;
else l=mid+1,t=tr[t].rs;
}
return t;
}
int findfapos(int rt,int x)
{
int t=findpos(rt,x);
if(tr[t].fa==x)return t;
else return findfapos(rt,tr[t].fa);
}
void modify(int &k,int fr,int l,int r,int x,int fa)
{
k=++cnt;if(l==r){tr[k].fa=fa;tr[k].d=tr[fr].d;return ;}
tr[k].ls=tr[fr].ls;tr[k].rs=tr[fr].rs;
int mid=l+r>>1;
if(x<=mid)modify(tr[k].ls,tr[fr].ls,l,mid,x,fa);
else modify(tr[k].rs,tr[fr].rs,mid+1,r,x,fa);
}
void add(int &k,int fr,int l,int r,int x)
{
k=++cnt;
if(l==r){tr[k].fa=x;tr[k].d=tr[fr].d+1;return ;}
tr[k].ls=tr[fr].ls;tr[k].rs=tr[fr].rs;
int mid=l+r>>1;
if(x<=mid)add(tr[k].ls,tr[fr].ls,l,mid,x);
else add(tr[k].rs,tr[fr].rs,mid+1,r,x);
}
inline void merger(int rt,int x,int y)
{
int fax=findfapos(rt,x),fay=findfapos(rt,y);
if(fax==fay)return ;
if(tr[fax].d>tr[fay].d)swap(fax,fay);
modify(root[rt],root[rt-1],1,n,tr[fax].fa,tr[fay].fa);
if(tr[fax].d==tr[fay].d)add(root[rt],root[rt],1,n,tr[fay].fa);
}
inline void query(int rt,int x,int y)
{
int fax=findfapos(rt,x),fay=findfapos(rt,y);
if(fax==fay)puts("1"),ans=1;else puts("0"),ans=0;
}
int main()
{
n=read();m=read();
build(root[0],1,n);
for(int i=1,opt,x,y;i<=m;++i)
{
opt=read();x=read()^ans;
root[i]=root[i-1];
if(opt==1)y=read()^ans,merger(i,x,y);
else if(opt==2)root[i]=root[x];
else y=read()^ans,query(i,x,y);
}
return 0;
}
/*
5 6
1 1 2
3 1 2
2 1
3 0 3
2 1
3 1 2
1
0
1
*/