AHOI2005航线规划 bzoj1969 洛谷2542 LCT 缩点

题意：链接

题解： 

LCT中的每一条边表示一个双连通分量，可以像水管局长那样离线处理，然后每次加点。如果加上这个点之后没有形成环，那么直接像link那样加上就行。如果形成环，就在LCT上暴力缩点。据说可以证明复杂度是对的，但是我现在还不太明白这个过程是复杂度，有巨佬说暴力修改的次数不超过nlogn。

暴力缩点并不需要像我想象中那样去写tarjan，代码中有一种递归+并查集来实现缩点。这个缩点的道理是，如果原来的splay加上这条边形成了一个环，那么我们对makeroot(x)之后，x在原树的子树就到了它的右子树（c[x][1]），那么我们只需要把此时整个右子树加上x缩成一个点，即把它们在并查集中的父节点均变为x在并查集中的父节点。询问相当于子树大小减一。

然后是代码

#include <bits/stdc++.h>
using namespace std;

int n,m,f[100001],fa[100010],l[100010],r[100010],xx,ji;//fa是并查集的 
int sz[100010],c[100010][2],rev[100100],st[100100],ans[100100];
struct node
{
	int x,y,opt;
}a[200010],t[200010];
int cmp(node x,node y)
{
	if(x.x==y.x)
	return x.y<y.y;
	else
	return x.x<y.x;
}
void div(int le,int ri,int x)
{
	int mid;
	while(le<=ri)
	{
		mid=(le+ri)>>1;
		if(a[mid].y==x)
		{
			a[mid].opt=1;
			return;
		}
		else if(a[mid].y<x)
		le=mid+1;
		else
		ri=mid-1;
	}
}
int getr(int x)
{
	if(x==fa[x])
	return x;
	else
	{
		fa[x]=getr(fa[x]);
		return fa[x];
	}
} 
void pushup(int x)
{
	sz[x]=sz[c[x][0]]+sz[c[x][1]]+1;
}
int nroot(int x)
{
	return c[f[x]][0]==x||c[f[x]][1]==x;
}
void pushdown(int x)
{
	if(rev[x])
	{
		swap(c[x][0],c[x][1]);
		rev[c[x][0]]^=1;
		rev[c[x][1]]^=1;
		rev[x]=0;
	}
}
void rotate(int x)
{
	int y=f[x],z=f[y],k=c[y][1]==x,w=c[x][!k];
	if(nroot(y))
	c[z][c[z][1]==y]=x;
	c[x][!k]=y;
	c[y][k]=w;
	if(w)
	f[w]=y;
	f[y]=x;
	f[x]=z;
	pushup(y);
	pushup(x);
}
void splay(int x)
{
	int y=x,z=0;
	st[++z]=y;
	while(nroot(y))
	{
		y=f[y];
		st[++z]=y;
	}
	while(z)
	pushdown(st[z--]);
	while(nroot(x))
	{
		y=f[x],z=f[y];
		if(nroot(y))
		{
			if(c[z][0]==y ^ c[y][0]==x)
			rotate(y);
			else
			rotate(x);
		}
		rotate(x);
	}
	pushup(x);
}
void access(int x)
{
	int y=0;
	while(x!=0)
	{
		splay(x);
		c[x][1]=y;
		pushup(x);
		y=x;
		f[y]=getr(f[x]);
		x=f[y];
	}
}
void makeroot(int x)
{
	access(x);
	splay(x);
	rev[x]^=1;
}
int findroot(int x)
{
	access(x);
	splay(x);
	pushdown(x);
	while(c[x][0])
	{
		x=c[x][0];
		pushdown(x);
	}
	splay(x);
	return x;
}
void del(int x,int y)
{
	if(x)
	{
		fa[x]=y;
		del(c[x][0],y);
		del(c[x][1],y);
	}	
}
void merge(int x,int y)
{
	if(x==y)
	return;
	makeroot(x);
	if(findroot(y)!=x)
	{
		f[x]=y;
		return;
	}	
	del(c[x][1],x);//递归缩点 	
	c[x][1]=0;//删点 
	pushup(x);
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;++i)
	{
		scanf("%d%d",&a[i].x,&a[i].y);
		if(a[i].x>a[i].y)
		swap(a[i].x,a[i].y);
	}
	sort(a+1,a+m+1,cmp);
	l[a[1].x]=1;
	r[a[m].x]=m;
	for(int i=2;i<=m;++i)
	{
		if(a[i-1].x!=a[i].x)
		{
			r[a[i-1].x]=i-1;
			l[a[i].x]=i;
		}
	}
	for(int i=1;i<=n;++i)
	fa[i]=i;
	for(int i=1;xx!=-1;++i)
	{
		scanf("%d",&xx);
		if(xx==-1)
		break;
		ji=i;
		scanf("%d%d",&t[i].x,&t[i].y);
		if(t[i].x>t[i].y)
		swap(t[i].x,t[i].y);
		t[i].opt=xx;
		if(t[i].opt==0)
		{
			div(l[t[i].x],r[t[i].x],t[i].y);
		}
	}
	for(int i=1;i<=m;++i)
	{
		if(a[i].opt==0)
		{
			merge(getr(a[i].x),getr(a[i].y));
		}
	}
	int i;
	for(i=0;ji;--ji)
	{
		int fx=getr(t[ji].x),fy=getr(t[ji].y);
		if(t[ji].opt)
		{
			makeroot(fy);
			access(fx);
			splay(fx);
			ans[++i]=sz[fx]-1;
		}
		else
		merge(fx,fy);
	}
	while(i)
	printf("%d\n",ans[i--]);
	return 0;
}