题面
Description
Jiajia和Wind是一对恩爱的夫妻,并且他们有很多孩子。
某天,Jiajia、Wind和孩子们决定在家里玩捉迷藏游戏。
他们的家很大且构造很奇特,由N个屋子和N-1条双向走廊组成,
这N-1条走廊的分布使得任意两个屋子都互相可达。
游戏是这样进行的,孩子们负责躲藏,Jiajia负责找,而Wind负责操纵这N个屋子的灯。
在起初的时候,所有的灯都没有被打开。
每一次,孩子们只会躲藏在没有开灯的房间中,
但是为了增加刺激性,孩子们会要求打开某个房间的电灯或者关闭某个房间的电灯。
为了评估某一次游戏的复杂性,
Jiajia希望知道可能的最远的两个孩子的距离(即最远的两个关灯房间的距离)。
我们将以如下形式定义每一种操作:
Input
输入文件hide.in第一行包含一个整数N,表示房间的个数,房间将被编号为1,2,3…N的整数。
接下来N-1行每行两个整数a, b,表示房间a与房间b之间有一条走廊相连。
接下来一行包含一个整数Q,表示操作次数。.
接着Q行,每行一个操作,如上文所示。
Output
对于每一个操作Game,输出一个非负整数到hide.out,表示最远的两个关灯房间的距离。
若只有一个房间是关着灯的,输出0;若所有房间的灯都开着,输出-1。
Sample Input
8
1 2
2 3
3 4
3 5
3 6
6 7
6 8
7
G
C 1
G
C 2
G
C 1
G
Sample Output
4
3
3
4
Hint
对于20%的数据,N≤50, Q≤100;
对于60%的数据,N≤3000, Q≤10000;
对于100%的数据,N≤100000, Q≤500000。
题目分析
动态询问多组点对信息,所以使用动态点分治。
根据点分治的基本思路,
对于一个重心,我们通常记录经过重心的路径来统计答案。
这里我们沿用这种思路,设亮着的房间为白点,黑着的房间为黑点。
对于某个重心\(w\),记录每个子树\(v\)中的黑点到\(w\)的最长链(每个子树只能贡献一个,避免重复)。
计算该重心对答案的贡献时直接取堆\(c[w]\)中的最大值+次大值即可。
由于有删除操作,我们要能够迅速更新最长链的方法。
于是,我们对每个节点\(v\),再开一个堆\(up\),记录以\(v\)为根的子树到\(v\)的\(parent\)的最长链。
注意:
上面的\(parent\)指的是点分树上的父亲。
然后再维护一个\(ans\)堆记录答案。
还有一个注意事项,由于要支持堆中的删除操作,需要额外建一个堆维护删除信息。
代码实现
#include
#include
#include
#include
#include
#include
#include
#include
#define MAXN 0x7fffffff
typedef long long LL;
const int N=100005;
using namespace std;
inline int Getint(){register int x=0,f=1;register char ch=getchar();while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}while(isdigit(ch)){x=x*10+ch-'0';ch=getchar();}return x*f;}
int h[N],cnt;
struct Edge{int to,next;}g[N<<1];
void AddEdge(int x,int y){g[++cnt].to=y,g[cnt].next=h[x],h[x]=cnt;}
int fa[N],dep[N];
int size[N],son[N],tp[N];
void Dfs1(int x){
size[x]=1;
for(int i=h[x];i;i=g[i].next){
int to=g[i].to;
if(to==fa[x])continue;
fa[to]=x,dep[to]=dep[x]+1;
Dfs1(to),size[x]+=size[to];
if(size[to]>size[son[x]])son[x]=to;
}
}
void Dfs2(int x,int top){
tp[x]=top;
if(son[x])Dfs2(son[x],top);
for(int i=h[x];i;i=g[i].next){
int to=g[i].to;
if(to==fa[x]||to==son[x])continue;
Dfs2(to,to);
}
}
int LCA(int x,int y){
while(tp[x]^tp[y]){
if(dep[tp[x]]q,del;
void push(int x){q.push(x);}
void erase(int x){del.push(x);}
int top(){
while(del.size()&&q.top()==del.top())q.pop(),del.pop();
return q.top();
}
void Pop(){
while(del.size()&&q.top()==del.top())q.pop(),del.pop();
q.pop();
}
int sec(){
int tmp=top();Pop();
int x=top();push(tmp);
return x;
}
int size(){return q.size()-del.size();}
}c[N],up[N],ans;
int n,num,sum,prt[N];
int rt,sz[N],mx[N];
bool vis[N],light[N];
void to_ans(int v,bool f){
if(c[v].size()>1){
int x=c[v].top()+c[v].sec();
f?ans.push(x):ans.erase(x);
}
}
void Explore(int x,int fa,int top){
up[top].push(dis(x,prt[top]));
for(int i=h[x];i;i=g[i].next){
int to=g[i].to;
if(vis[to]||to==fa)continue;
Explore(to,x,top);
}
}
void Get_size(int x,int fa) {
sz[x]=1;
for(int i=h[x];i;i=g[i].next) {
int to=g[i].to;
if(to==fa||vis[to]) continue;
Get_size(to,x);
sz[x]+=sz[to];
}
}
void Get_root(int x,int fa){
mx[x]=0;
for(int i=h[x];i;i=g[i].next){
int to=g[i].to;
if(to==fa||vis[to])continue;
Get_root(to,x);
mx[x]=max(mx[x],sz[to]);
}
mx[x]=max(mx[x],sum-sz[x]);
if(mx[rt]>mx[x])rt=x;
}
void Solve(int x){
vis[x]=1,c[x].push(0);
for(int i=h[x];i;i=g[i].next){
int to=g[i].to;
if(vis[to]||to==prt[x])continue;
Get_size(to,0),sum=sz[to];
rt=0,Get_root(to,0);
prt[rt]=x;
Explore(rt,0,rt);
c[x].push(up[rt].top());
Solve(rt);
}
to_ans(x,1);
}
void Change(int x,bool f){
f?(num--):(num++);
to_ans(x,0);
f?c[x].erase(0):c[x].push(0);
to_ans(x,1);
for(int i=x;prt[i];i=prt[i]){
to_ans(prt[i],0);
if(up[i].size())c[prt[i]].erase(up[i].top());
f?up[i].erase(dis(x,prt[i])):up[i].push(dis(x,prt[i]));
if(up[i].size())c[prt[i]].push(up[i].top());
to_ans(prt[i],1);
}
}
int main(){
num=n=Getint();
for(int i=1;i