BZOJ 3674: 可持久化并查集加强版/BZOJ 3673: 可持久化并查集 by zky 可持久化线段树

3674: 可持久化并查集加强版

Time Limit: 15 Sec  Memory Limit: 256 MB
Submit: 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

Sample Output

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
*/