BZOJ 4129

http://www.lydsy.com/JudgeOnline/problem.php?id=4129

还是树上带修莫队，

只不过要离散化一下，

对于离散化后的值域还要分块一下

每次询问可以进行分块查询

复杂度O(n√n)

#include<cstdio>
#include<algorithm>
#include<cmath>
#define VIS(now) for(register int e=las[now];e;e=nxt[e])
#define FOR(i,s,t) for(register int i=s;i<=t;++i)
#define up(l,r) for(register int j=l;j!=r;j=fa[j])
using std::swap;
using std::lower_bound;
using std::sort;
using std::unique;
const int N=200011;
int n,m,blo,dfn_num,t,l,r,x,y,tot,times,now,lca;
int las[N],nxt[N],to[N],a[N],b[N<<1];
int cnt[N<<1],g[N],bo[N],ans[N]; 
int xu[N],top[N],sz[N],dep[N],fa[N];
struct qs{
	int l,r,t,id;
	inline bool operator<(qs A)const{
		if(xu[l]/blo!=xu[A.l]/blo)return xu[l]<xu[A.l];
		if(xu[r]/blo!=xu[A.r]/blo)return xu[r]<xu[A.r];
		if(t/blo!=A.t/blo)return t<A.t;
		return id<A.id;
	}
}q[N];
struct cg{int pos,num;}c[N];
inline void add(int x,int y){
	nxt[++tot]=las[x];las[x]=tot;to[tot]=y;
}
inline void dfs1(int now){
	sz[now]=1;
	VIS(now)
		if(!dep[to[e]]){
			dep[to[e]]=dep[now]+1;
			fa[to[e]]=now;
			dfs1(to[e]);
			sz[now]+=sz[to[e]];
		}
}
inline void dfs2(int now,int chain){
	xu[now]=++dfn_num;
	top[now]=chain;
	register int i=0;
	VIS(now)if(sz[to[e]]>sz[i]&&to[e]!=fa[now])i=to[e];
	if(!i)return;
	dfs2(i,chain);
	VIS(now)if(to[e]!=i&&to[e]!=fa[now])dfs2(to[e],to[e]);
}
inline int LCA(int x,int y){
	for(;top[x]!=top[y];dep[top[x]]>dep[top[y]]?x=fa[top[x]]:y=fa[top[y]]);
	return dep[x]<dep[y]?x:y;
}
inline void disc_init(){
	sort(b+1,b+b[0]+1);
	b[0]=unique(b+1,b+b[0]+1)-b-1; 
	FOR(i,1,n)a[i]=lower_bound(b+1,b+b[0]+1,a[i])-b;
	FOR(i,1,times)c[i].num=lower_bound(b+1,b+b[0]+1,c[i].num)-b;
}
inline void add(int now){
	int col;
	bo[now]^=1;
	col=a[now];
	if(bo[now]){if(!cnt[col])++g[(col-1)/blo+1];++cnt[col];}
	else{--cnt[col];if(!cnt[col])--g[(col-1)/blo+1];}
} 
inline void change(int now){
	int pos=c[now].pos;
	int col;
	if(bo[pos]){
		col=a[pos];--cnt[col];if(!cnt[col])--g[(col-1)/blo+1];
		col=c[now].num;if(!cnt[col])++g[(col-1)/blo+1];++cnt[col];
	}
	swap(a[pos],c[now].num);
}
inline int query(){
	register int i;
	for(i=1;g[i]==blo;++i);
	for(register int j=(i-1)*blo+1;;++j)
		if(!cnt[j])return j;
}
int main(){
	scanf("%d%d",&n,&m);
	blo=pow(n,2.00/3.00);
	b[++b[0]]=0;
	FOR(i,1,n){
		scanf("%d",a+i);
		b[++b[0]]=a[i];b[++b[0]]=a[i]+1;
	}
	FOR(i,2,n){
		scanf("%d%d",&x,&y);
		add(x,y);add(y,x);
	}
	dep[1]=1;dfs1(1);dfs2(1,1);
	tot=0;
	FOR(i,1,m){
		scanf("%d%d%d",&t,&l,&r);
		if(t==1){
			if(xu[l]>xu[r])swap(l,r);
			q[++tot]=(qs){l,r,times,tot};		
		}
		else{
			c[++times]=(cg){l,r};
			b[++b[0]]=r;
			b[++b[0]]=r+1;
		}
	}
	disc_init();
	sort(q+1,q+tot+1);
	l=1;r=1;now=0;
	FOR(i,1,tot){
		lca=LCA(l,q[i].l);up(l,lca)add(j);up(q[i].l,lca)add(j);
		lca=LCA(r,q[i].r);up(r,lca)add(j);up(q[i].r,lca)add(j);
		lca=LCA(q[i].l,q[i].r);add(lca);
		while(now<q[i].t)change(++now);while(q[i].t<now)change(now--);
		ans[q[i].id]=query();add(lca);
		l=q[i].l;r=q[i].r;
	}
	FOR(i,1,tot)printf("%d
",b[ans[i]]);
	return 0;
}