点了动态点分治的科技树,这道题是树形态不变的动态点分治,形态变化的话...待会补
考虑点分治过程中的这样一种结构:按递归层次把当前层的重心与上层重心互相连接,这就是点分治树,容易看出它的树高只有$O(\log_2n)$
对这棵树进行讨论,先考虑没有修改怎么做,对于点分树上的一个点$x$,我们遍历它的每个儿子$u$的子树并找到$u$的子树内的所有点到$x$的最大距离,记为$far_u$,那么$x$的所有儿子$u$的$far_u$中最大和次大加起来就是经过$x$的答案
因为有修改所以我们需要把$u$子树内的所有点到$x$的距离用一个堆存在节点$u$(不妨称此堆为堆$1$),同时我们需要把$x$的所有儿子$u$的堆$1$堆顶用一个堆存在$x$(不妨称为堆$2$),再用堆$3$存所有(堆$2$的最大值和次大值之和),此时堆$3$的堆顶就是答案
到这里,修改就变得很简单了,改堆$1$的同时更新堆$2$堆$3$即可
实现的时候如果某个节点是可用的,它的堆$2$内要有一个$0$表示路径从$x$出发,修改时也要记得加$0$或删$0$
#include<stdio.h>
#include<queue>
using namespace std;
const int inf=2147483647;
struct heap{
priority_queue<int>h,d;
void push(int v){h.push(v);}
void erase(int v){d.push(v);}
void adj(){
while(!d.empty()&&h.top()==d.top()){
h.pop();
d.pop();
}
}
int top(){
adj();
return h.top();
}
void pop(){
adj();
h.pop();
}
int sec(){
int x,y;
x=top();
pop();
y=top();
push(x);
return y;
}
int size(){return h.size()-d.size();}
}h1[100010],h2[100010],al;
int h[100010],to[200010],nex[200010],M;
void add(int a,int b){
M++;
to[M]=b;
nex[M]=h[a];
h[a]=M;
}
int fa_[100010][17],dep[100010];
void dfs(int x){
for(int i=h[x];i;i=nex[i]){
if(to[i]!=fa_[x][0]){
fa_[to[i]][0]=x;
dep[to[i]]=dep[x]+1;
dfs(to[i]);
}
}
}
int lca(int x,int y){
int i;
if(dep[x]<dep[y])swap(x,y);
for(i=16;i>=0;i--){
if(dep[fa_[x][i]]>=dep[y])x=fa_[x][i];
}
if(x==y)return x;
for(i=16;i>=0;i--){
if(fa_[x][i]!=fa_[y][i]){
x=fa_[x][i];
y=fa_[y][i];
}
}
return fa_[x][0];
}
int dis(int x,int y){return dep[x]+dep[y]-2*dep[lca(x,y)];}
int siz[100010],n;
bool v[100010];
void dfs1(int fa,int x){
n++;
siz[x]=1;
for(int i=h[x];i;i=nex[i]){
if(!v[to[i]]&&to[i]!=fa){
dfs1(x,to[i]);
siz[x]+=siz[to[i]];
}
}
}
int mn,cn;
void dfs2(int fa,int x){
int i,k;
k=0;
for(i=h[x];i;i=nex[i]){
if(!v[to[i]]&&to[i]!=fa){
dfs2(x,to[i]);
k=max(k,siz[to[i]]);
}
}
k=max(k,n-siz[x]);
if(k<mn){
mn=k;
cn=x;
}
}
void dfs3(heap&p,int s,int fa,int x){
p.push(dis(x,s));
for(int i=h[x];i;i=nex[i]){
if(!v[to[i]]&&to[i]!=fa)dfs3(p,s,x,to[i]);
}
}
int fa[100010];
int solve(int f,int x){
n=0;
dfs1(0,x);
mn=inf;
dfs2(0,x);
x=cn;
if(f)dfs3(h1[x],f,0,x);
fa[x]=f;
v[x]=1;
h2[x].push(0);
for(int i=h[x];i;i=nex[i]){
if(!v[to[i]])h2[x].push(h1[solve(x,to[i])].top());
}
if(h2[x].size()>1)al.push(h2[x].top()+h2[x].sec());
return x;
}
void modify(int x,int v,bool f){
if(h2[fa[x]].size()>1)al.erase(h2[fa[x]].top()+h2[fa[x]].sec());
if(h2[fa[x]].size()>0&&h1[x].size()>0)h2[fa[x]].erase(h1[x].top());
f?h1[x].push(v):h1[x].erase(v);
if(h1[x].size()>0)h2[fa[x]].push(h1[x].top());
if(h2[fa[x]].size()>1)al.push(h2[fa[x]].top()+h2[fa[x]].sec());
}
void change(int x){
v[x]^=1;
if(h2[x].size()>1)al.erase(h2[x].top()+h2[x].sec());
v[x]?h2[x].push(0):h2[x].erase(0);
if(h2[x].size()>1)al.push(h2[x].top()+h2[x].sec());
for(int i=x;fa[i];i=fa[i])modify(i,dis(fa[i],x),v[x]);
}
int main(){
int n,m,i,j,x,y,sum;
char s[5];
scanf("%d",&n);
for(i=1;i<n;i++){
scanf("%d%d",&x,&y);
add(x,y);
add(y,x);
}
dep[1]=1;
dfs(1);
for(j=1;j<17;j++){
for(i=1;i<=n;i++)fa_[i][j]=fa_[fa_[i][j-1]][j-1];
}
solve(0,1);
scanf("%d",&m);
sum=n;
while(m--){
scanf("%s",s);
if(s[0]=='G')
printf("%d\n",sum<2?sum-1:al.top());
else{
scanf("%d",&x);
v[x]?(sum--):(sum++);
change(x);
}
}
}