线段树分治 - 「离线可过」动态图连通性（LOJ121）

传送门

Analysis

线段树分治入门

所谓线段树分治就是先把操作分成logn个连续不相交的区间分别维护信息。
最后按线段树从上到下再从左到右的遍历方式一起统计答案。
这道题可以按时间建树，每次相当于在一段区间里增加边。
这样我们就避免了删边的操作

由于每次回溯的时候要撤销本次操作，我们用并查集就不可以路径压缩
又因为我们每次都是后加入的先撤销，所以并查集可做，用一个栈记录一下即可

总的时间复杂度 $m\ast logm\ast logn$
有m条边，每条边的存在时间在线段树上被划分为了log m个区间
每次合并的复杂度是O（logn），撤销是O（1）的

Code

#include<bits/stdc++.h>
#define in read()
#define re register
using namespace std;
inline int read(){
	char ch;int f=1,res=0;
	while((ch=getchar())<'0'||ch>'9') if(ch=='-') f=-1;
	while(ch>='0'&&ch<='9'){
		res=(res<<1)+(res<<3)+(ch^48);
		ch=getchar();
	}
	return f==1?res:-res;
}
const int N=5009,M=5e5+10;
int n,m;
int link[N][N];
struct query{
	int op,x,y;
	query(int op=0,int x=0,int y=0):op(op),x(x),y(y){}
}q[M];
#define lc (k<<1)
#define rc (k<<1|1)
struct segment{
	vector<int> v;
}tree[M<<2];
void modify(int k,int l,int r,int x,int y,int val){
	if(x<=l&&r<=y){
		tree[k].v.push_back(val);
		return;
	}
	int mid=l+r>>1;
	if(y<=mid) modify(lc,l,mid,x,y,val);
	else if(x>mid) modify(rc,mid+1,r,x,y,val);
	else modify(lc,l,mid,x,y,val),modify(rc,mid+1,r,x,y,val);
}
int fa[N],sze[N];
inline int getfa(int x){
	return x==fa[x]?x:getfa(fa[x]);
}
int top=0;
pair<int,int> s[M];
inline void merge(int x,int y){
	x=getfa(x);y=getfa(y);
	if(x==y) return;
	if(sze[x]<sze[y]) swap(x,y);
	fa[y]=x;
	s[++top]=make_pair(y,sze[y]);
	s[++top]=make_pair(x,sze[x]);
	sze[x]+=sze[y];
}
inline void delet(int cur){
	while(top>cur){
		fa[s[top].first]=s[top].first;
		sze[s[top].first]=s[top].second;
		top--;
	}
}
void dfs(int k,int l,int r){
	int cur=top;
	for(re int i=0;i<tree[k].v.size();++i) merge(q[tree[k].v[i]].x,q[tree[k].v[i]].y);
	if(l==r){
		if(q[l].op==2) printf("%c\n",getfa(q[l].x)==getfa(q[l].y)?'Y':'N');
	}
	else{
		int mid=l+r>>1;
		dfs(lc,l,mid);dfs(rc,mid+1,r);
	}
	delet(cur);
}
int main(){
	n=in;m=in;
	for(re int i=1;i<=n;++i) fa[i]=i,sze[i]=1;
	for(re int i=1;i<=m;++i){
		int op=in,x=in,y=in;
		if(x>y) swap(x,y);
		if(op==0) link[x][y]=i;
		if(op==1) modify(1,1,m,link[x][y],i,i),link[x][y]=0;
		q[i]=(query){op,x,y};
	}
	for(re int i=1;i<=n;++i)
		for(re int j=1;j<=n;++j)
			if(link[i][j]) modify(1,1,m,link[i][j],m,link[i][j]);
	dfs(1,1,m);	
	return 0;
}