[bzoj3631][JLOI2014]松鼠的新家 树链剖分

Description

松鼠的新家是一棵树，前几天刚刚装修了新家，新家有n个房间，并且有n-1根树枝连接，每个房间都可以相互到达，且俩个房间之间的路线都是唯一的。天哪，他居然真的住在“树”上。松鼠想邀请小熊维尼前来参观，并且还指定一份参观指南，他希望维尼能够按照他的指南顺序，先去a1，再去a2，……，最后到an，去参观新家。

可是这样会导致维尼重复走很多房间，懒惰的维尼不听地推辞。可是松鼠告诉他，每走到一个房间，他就可以从房间拿一块糖果吃。维尼是个馋家伙，立马就答应了。

现在松鼠希望知道为了保证维尼有糖果吃，他需要在每一个房间各放至少多少个糖果。因为松鼠参观指南上的最后一个房间an是餐厅，餐厅里他准备了丰盛的大餐，所以当维尼在参观的最后到达餐厅时就不需要再拿糖果吃了。

Input

第一行一个整数n，表示房间个数

第二行n个整数，依次描述a1-an

接下来n-1行，每行两个整数x，y，表示标号x和y的两个房间之间有树枝相连。

Output

一共n行，第i行输出标号为i的房间至少需要放多少个糖果，才能让维尼有糖果吃。

Sample Input

5
1 4 5 3 2
1 2
2 4
2 3
4 5

Sample Output

1
2
1
2
1

HINT

2<= n <=300000

Source

 

思维难度:NOIP+

代码难度:省选

算法:树链剖分

#include<cstdio>
#include<iostream>
using namespace std;
const int Maxn=300005;
int fa[Maxn],deep[Maxn],sz[Maxn],son[Maxn],top[Maxn],idx[Maxn],h[Maxn],cnt,tcnt,n,root,a[Maxn],cut[Maxn];
struct node{
	int v,next;
}e[Maxn*2];
struct tree{
	int x,c,l,r;
}t[Maxn*4];
void add(int u,int v){
	cnt++;
	e[cnt].v=v;
	e[cnt].next=h[u];
	h[u]=cnt;
}
inline int read(){
	int x=0,f=1;char c=getchar();
	while(c<'0'||c>'9'){
		if(c=='-')f=-1;
		c=getchar();
	}
	while(c<='9'&&c>='0'){
		x=x*10+c-'0';
		c=getchar();
	}
	return x*f;
}
inline int dfs1(int u){
	sz[u]=1;
	int heavy=0,maxn=0;
	for(int i=h[u];i;i=e[i].next){
		int v=e[i].v;
		if(deep[v])continue;
		deep[v]=deep[u]+1;fa[v]=u;
		int tmp=dfs1(v);
		if(tmp>maxn){
			maxn=tmp;
			heavy=v;
		}
		sz[u]+=tmp;
	}
	son[u]=heavy;
	return sz[u];
}
inline void dfs2(int u){
	idx[u]=++tcnt;
	if(!son[u])return;
	top[son[u]]=top[u];dfs2(son[u]);
	for(int i=h[u];i;i=e[i].next){
		int v=e[i].v;
		if(idx[v])continue;
		top[v]=v;
		dfs2(v);
	}
}
inline int lson(int rt){
	return rt<<1;
}
inline int rson(int rt){
	return rt<<1|1;
}
inline int pushup(int rt){
	t[rt].x=t[lson(rt)].x+t[rson(rt)].x;
}
inline void build(int l,int r,int rt){
	t[rt].l=l;t[rt].r=r;
	if(l==r)return;
	int mid=l+r>>1;
	build(l,mid,lson(rt));
	build(mid+1,r,rson(rt));
	pushup(rt);
}
inline void pushdown(int rt){
	if(t[rt].l==t[rt].r)return;
	if(t[rt].c==0)return;
	t[lson(rt)].c+=t[rt].c;
	t[rson(rt)].c+=t[rt].c;
	t[lson(rt)].x+=t[rt].c*(t[lson(rt)].r-t[lson(rt)].l+1);
	t[rson(rt)].x+=t[rt].c*(t[rson(rt)].r-t[rson(rt)].l+1);
	t[rt].c=0;
}
inline void update(int l,int r,int rt){
	if(t[rt].l>=l&&t[rt].r<=r){
		t[rt].x+=(t[rt].r-t[rt].l+1);
		t[rt].c+=1;
		return;
	}
	pushdown(rt);
	int mid=t[rt].l+t[rt].r>>1;
	if(l<=mid){
		update(l,r,lson(rt));
	}
	if(r>mid){
		update(l,r,rson(rt));
	}
	pushup(rt);
}
inline int query(int l,int rt){
	if(t[rt].l==l&&t[rt].r==l){
		return t[rt].x;
	}
	pushdown(rt);
	int mid=t[rt].l+t[rt].r>>1;
	if(l<=mid){
		return query(l,lson(rt));
	}
	else{
		return query(l,rson(rt));
	}
}
inline void Tupdate(int x,int y){
	while(top[x]!=top[y]){
		if(deep[top[x]]<deep[top[y]])swap(x,y);
		update(idx[top[x]],idx[x],1);
		x=fa[top[x]];
	}
	if(idx[x]>idx[y])swap(x,y);update(idx[x],idx[y],1);
}
int main(){
	int u,v;
	n=read();
	for(int i=1;i<=n;++i)a[i]=read();root=a[1];
	for(int i=1;i<=n-1;++i){
		u=read();v=read();
		add(u,v);add(v,u);
	}
	deep[root]=1;top[root]=root;fa[root]=0;
	dfs1(root);dfs2(root);build(1,n,1);
	for(int i=2;i<=n;++i){
		Tupdate(a[i],a[i-1]);
		cut[a[i]]++;
	}
	for(int i=1;i<=n;++i){
		printf("%d\n",query(idx[i],1)-cut[i]);
	}
	return 0;
}