[JZOJ3690] 【CF418D】Big Problems for Organizers

题目

题目大意

给你一棵树，然后有一堆询问，每次给出两个点。
问所有点到两个点中最近点的距离的最大值。

---

正解

本来打了倍增，然后爆了，也懒得调……
显然可以在两个点之间的路径的中点处割开，一边归一个点管。
有个比较显然的思路是DP，设(f_x)表示(x)子树内的最远点，(g_x)向父亲那边走的最远点。
然后就可以倍增搞，合并一下……
代码复杂度极高。

然后有个简单又自然的思路是直接打(LCT)。直接把中间那条边断掉，然后求最远点即可。
想法倒是简单自然，接下来的问题就是，如何求最远点？
每条链可以维护子树到链顶的距离。设最远距离为(len)，则合并的时候，就是左区间的(len)、右区间的(len)加上左区间的(siz)加(1)、中间的答案（包括虚边转移过来的）加左区间的(siz)。这三个东西取最大值。
至于虚边信息，由于维护的是最大值，不能像和一样加加减减，所以要用multiset维护（或者手写数据结构也可以哈）。
然后你就会发现一个bug——翻转的时候怎么办？
(len)的计算是不满足交换律的，所以我们要维护一个(len2)，表示(len)反过来的样子。合并的时候跟上面相反就是了。换句话来说，就是在重链管的整棵子树中，到达链底的最远距离。
翻转的时候直接将两个交换就行了。

---

代码

using namespace std;
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <set>
#define N 100010
#define INF 1000000000
int n,m;
struct Node{
	Node *fa,*c[2];
	bool is_root,rev;
	int siz,len1,len2,mx;
	multiset<int> s;
	inline void reverse(){
		swap(c[0],c[1]);
		swap(len1,len2);
		rev^=1;
	}
	inline void pushdown(){
		if (rev){
			c[0]->reverse();
			c[1]->reverse();
			rev=0;
		}
	}
	void push(){
		if (!is_root)
			fa->push();
		pushdown();
	}
	inline void update(){
		siz=c[0]->siz+c[1]->siz+1;
		len1=max(max(c[0]->len1,c[0]->siz+1+c[1]->len1),c[0]->siz+mx+1);
		len2=max(max(c[1]->len2,c[1]->siz+1+c[0]->len2),c[1]->siz+mx+1);
	}
	inline bool getson(){return fa->c[0]!=this;}
	inline void rotate(){
		Node *y=fa,*z=y->fa;
		if (y->is_root){
			is_root=1;
			y->is_root=0;
		}
		else
			z->c[y->getson()]=this;
		int k=getson();
		fa=z;
		y->c[k]=c[k^1],c[k^1]->fa=y;
		c[k^1]=y,y->fa=this;
		siz=y->siz,len1=y->len1,len2=y->len2;
		y->update();
	}
	inline void splay(){
		push();
		while (!is_root){
			if (!fa->is_root){
				if (getson()!=fa->getson())
					rotate();
				else
					fa->rotate();
			}
			rotate();
		}
	}
} d[N],*null=d;
inline Node *access(Node *x){
	Node *y=null;
	for (;x!=null;y=x,x=x->fa){
		x->splay();
		if (x->c[1]!=null){
			x->s.insert(x->c[1]->len1);
			x->mx=max(x->mx,x->c[1]->len1);
		}
		x->c[1]->is_root=1;
		x->c[1]=y;
		y->is_root=0;
		if (x->fa!=null){
			x->fa->s.erase(x->fa->s.find(x->len1));
			if (x->len1==x->fa->mx)
				x->fa->mx=(x->fa->s.empty()?-INF:*x->fa->s.rbegin());
		}
		x->update();
	}
	return y;
}
inline Node *mroot(Node *x){
	Node *t=access(x);
	t->reverse();
	return t;
}
inline void link(Node *u,Node *v){
	mroot(u),u->splay();
	Node *t=mroot(v);
	t->fa=u;
	u->s.insert(t->len1);
	if (t->len1>u->mx){
		u->mx=t->len1;
		u->update();
	}
}
Node *find(Node *t,int k){
	t->pushdown();
	if (k<=t->c[0]->siz)
		return find(t->c[0],k);
	if (k>t->c[0]->siz+1)
		return find(t->c[1],k-t->c[0]->siz-1);
	return t;
}
int main(){
	scanf("%d",&n);
	*null={null,null,null,0,0,0,-1,-1,-INF};
	for (int i=1;i<=n;++i)
		d[i]={null,null,null,1,0,1,0,0,-INF};
	for (int i=1;i<n;++i){
		int u,v;
		scanf("%d%d",&u,&v);
		link(&d[u],&d[v]);
	}
	scanf("%d",&m);
	for (int i=1;i<=m;++i){
		int u,v,l,ans=0;
		scanf("%d%d",&u,&v);
		mroot(&d[u]);
		Node *a=access(&d[v]),*b;
		a=find(a,a->siz>>1),a->splay();
		for (b=a->c[1];b->c[0]!=null;b=b->c[0])
			b->pushdown();
		b->splay();
		b->c[0]=null,b->update();
		ans=max(a->len1,b->len2);
		printf("%d
",ans);
		b->c[0]=a,b->update();
	}
	return 0;
}

---

总结

如果要让(LCT)支持换根操作，并且维护的不满足交换律的东西的时候，要维护它的反向信息。