[LOJ]#121. 「离线可过」动态图连通性 线段树分治+撤销并查集 | LCT

Description

这是一道被离线爆艹的模板题。
你要维护一张无向简单图。你被要求加入删除一条边及查询两个点是否连通。
0：加入一条边。保证它不存在。
1：删除一条边。保证它存在。
2：查询两个点是否联通。

题解：

方法一：
删边不好删，考虑避免删边操作。共有$M$个时间点，以时间为下标建线段树，线段树每个节点用一个$vector$维护当前时间段有哪些边存在。对于每条边求出其存在的时间段$[L,R]$，然后在这一段插入这条边。最后再$dfs$一次，到每个叶子节点时求答案。由于要回溯，所以要用可以撤销的并查集，不路径压缩，按秩合并保证复杂度。
这里的可撤销并查集并不是可以随意删除的，而是要保证后加入的边先被删除才可以使用。
方法二：
同样求出每条边存在的时间段$[L,R]$，然后用LCT维护以消失时间为边权的最大生成树即可。

代码（方法一）：

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define pa pair<int,int>
const int Maxn=5010;
const int Maxm=500010;
const int inf=2147483647;
int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9')x=(x<<3)+(x<<1)+(ch^48),ch=getchar();
	return x*f;
}
int n,m,tim[Maxn][Maxn],ans[Maxm];
struct Query{int x,y,t,id;}q[Maxm];int lq=0;
struct Edge{int x,y,l,r;}e[Maxm];int le=0;
int fa[Maxn],top=0,dep[Maxn],p1[Maxn],p2[Maxn];
int findfa(int x){return((fa[x]==x)?x:findfa(fa[x]));}
void Merge(int x,int y)
{
	int fx=findfa(x),fy=findfa(y);
	if(fx==fy)return;
	if(dep[fx]<=dep[fy])
	{
		fa[fx]=fy;
		top++;p1[top]=fx;p2[top]=dep[fy];
		dep[fy]=max(dep[fy],dep[fx]+1);
	}
	else
	{
		fa[fy]=fx;
		top++;p1[top]=fy;p2[top]=dep[fx];
		dep[fx]=max(dep[fx],dep[fy]+1);
	}
}
void Delete(int p)
{
	while(top>p)
	{
		dep[fa[p1[top]]]=p2[top];
		fa[p1[top]]=p1[top];
		top--;
	}
}
vector<int>h[Maxm<<1];
vector<Query>hh[Maxm<<1];
struct Seg{int l,r,lc,rc;}tr[Maxm<<1];
int tot=0,P[Maxm];
void build(int l,int r)
{
	int x=++tot;
	tr[x].l=l;tr[x].r=r;
	if(l==r){P[l]=x;return;}
	int mid=l+r>>1;
	tr[x].lc=tot+1,build(l,mid);
	tr[x].rc=tot+1,build(mid+1,r);
}
void insert(int x,int l,int r,int id)
{
	if(tr[x].l==l&&tr[x].r==r){h[x].push_back(id);return;}
	int lc=tr[x].lc,rc=tr[x].rc,mid=tr[x].l+tr[x].r>>1;
	if(r<=mid)insert(lc,l,r,id);
	else if(l>mid)insert(rc,l,r,id);
	else insert(lc,l,mid,id),insert(rc,mid+1,r,id);
}
void dfs(int x)
{
	int now=top;
	for(int i=0;i<h[x].size();i++)Merge(e[h[x][i]].x,e[h[x][i]].y);
	if(tr[x].l==tr[x].r)
	{
		for(int i=0;i<hh[x].size();i++)
		ans[hh[x][i].id]=((findfa(hh[x][i].x)==findfa(hh[x][i].y))?1:0);
	}
	if(tr[x].l!=tr[x].r)dfs(tr[x].lc),dfs(tr[x].rc);
	Delete(now);
}
int main()
{
	n=read(),m=read();
	for(int i=1;i<=m;i++)
	{
		int op=read(),x=read(),y=read();
		if(x>y)swap(x,y);
		if(!op)
		{
			tim[x][y]=++le;
			e[le].x=x;e[le].y=y;e[le].l=i;e[le].r=m;
		}
		else if(op==1)e[tim[x][y]].r=i;
		else q[++lq].x=x,q[lq].y=y,q[lq].t=i,q[lq].id=lq;
	}
	for(int i=1;i<=n;i++)fa[i]=i,dep[i]=1;
	build(1,m);
	for(int i=1;i<=le;i++)insert(1,e[i].l,e[i].r,i);
	for(int i=1;i<=lq;i++)hh[P[q[i].t]].push_back(q[i]);
	dfs(1);
	for(int i=1;i<=lq;i++)puts(ans[i]?"Y":"N");
}