LCT（维护边双）--bzoj1969: [Ahoi2005]LANE 航线规划

传送门

用$LCT$维护边双，删边不好操作，可以离线后变成加边。

如果两个点$u,v$已经联通，再加一条$(u,v)$的边就会让它变成边双，所以可以用并查集维护联通性。

当一条链变成边双就把它在$LCT$上缩成一个点，具体的可以直接在$splay$上$dfs$把所有点的父亲都修改成一个点，而做其他操作如$splay$和$access$时，用到的$f[x]$实际上也是$find(f[x])$。

注意这里两个并查集是不同的，一个是维护$LCT$上的代表点，一个是维护原树联通性的，加边的时候首先要把$u,v$变成$LCT$上的，再用维护联通性的并查集判断是不是要变成边双.

代码如下：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<map>
#define N 30005
#define M 100005
#define ls ch[x][0]
#define rs ch[x][1]
using namespace std;

inline int rd(){
	int x=0,f=1;char c=' ';
	while(c<'0' || c>'9') f=c=='-'?-1:1,c=getchar();
	while(c<='9' && c>='0') x=x*10+c-'0',c=getchar();
	return x*f;
}

int n,m,ans[40005];
int f[N],fa[N],ch[N][2],rev[N],sum[N],fa2[N];
pair<int,int> pr;
map< pair<int,int>,int > mp;

inline int get(int x){return ch[f[x]][1]==x;}
inline int isrt(int x){return ch[f[x]][0]!=x && ch[f[x]][1]!=x;}
inline void update(int x){sum[x]=sum[ls]+sum[rs]+1;}
inline void rever(int x){rev[x]^=1;swap(ls,rs);}
inline int find(int x){
	if(x==fa[x]) return x; return fa[x]=find(fa[x]);
}
inline int find2(int x){
	if(x==fa2[x]) return x; return fa2[x]=find2(fa2[x]);
}
inline void pushdown(int x){
	if(rev[x]){
		if(ls) rever(ls);
		if(rs) rever(rs);
		rev[x]=0;
	}
}
inline void pushup(int x){if(!isrt(x)) pushup(find(f[x]));pushdown(x);f[x]=find(f[x]);}

inline void rotate(int x){
	int old=f[x],oldf=f[old],wh=get(x);
	if(!isrt(old)) ch[oldf][get(old)]=x;
	ch[old][wh]=ch[x][wh^1]; f[ch[x][wh^1]]=old;
	ch[x][wh^1]=old; f[x]=oldf; f[old]=x;
	update(old); update(x);
}
inline void splay(int x){
	pushup(x);
	for(;!isrt(x);rotate(x))
		if(!isrt(f[x])) rotate(get(x)==get(f[x])?f[x]:x);
}

inline void access(int x){
	for(int y=0;x;y=x,x=find(f[x])) splay(x),rs=y,update(x); 
}
inline void makeroot(int x) {access(x);splay(x);rever(x);}
inline void split(int x,int y) {makeroot(x);access(y);splay(y);}
inline void link(int x,int y) {makeroot(x);f[x]=y;}
inline void cut(int x,int y) {split(x,y);f[x]=ch[y][0]=0;update(y);}
inline int findroot(int x) {access(x);splay(x); while(ls) pushdown(x),x=ls; return x;}

void dfs(int x,int y){
	fa[find(x)]=find(y);
	pushdown(x);
	if(ls) dfs(ls,y);
	if(rs) dfs(rs,y);
	ch[x][0]=ch[x][1]=0;
	return ;
}

struct EDGE{int to,fr;}edge[M];
struct qwq{int tp,x,y;}op[40005];

int main(){
	n=rd(); m=rd();
	for(int i=1;i<=n;i++) fa[i]=i,fa2[i]=i;
	for(int i=1;i<=m;i++) edge[i].fr=rd(),edge[i].to=rd();
	int tot=0;
	while(1){
		op[++tot].tp=rd(); if(op[tot].tp==-1) {--tot;break;}
		op[tot].x=rd(),op[tot].y=rd();
		if(op[tot].tp==0) mp[make_pair(op[tot].x,op[tot].y)]=mp[make_pair(op[tot].y,op[tot].x)]=1;
	}
	for(int i=1;i<=m;i++){
		if(mp.find(make_pair(edge[i].fr,edge[i].to))!=mp.end()) continue;
		int x=find(edge[i].fr),y=find(edge[i].to);
		if(x==y) continue;
		int u=find2(x),v=find2(y);
		if(u!=v) fa2[u]=v,link(x,y);
		else{
			split(x,y); dfs(y,y); update(y);
		}
	}
	for(int i=tot;i;i--){
		int fx=find(op[i].x),fy=find(op[i].y);
		if(op[i].tp==0){
			if(fx==fy) continue;
			int u=find2(fx),v=find2(fy);
			if(u!=v) fa2[u]=v,link(fx,fy);
			else split(fx,fy),dfs(fy,fy),update(fy);
		}
		else{
			split(fx,fy);
			ans[i]=sum[fy]-1;
		}
	}
	for(int i=1;i<=tot;i++) if(op[i].tp==1) printf("%d\n",ans[i]);
	return 0;
}