【bzoj 1036】 树的统计Count 【ZJOI2008】

Description

　　一棵树上有n个节点，编号分别为1到n，每个节点都有一个权值w。我们将以下面的形式来要求你对这棵树完成
一些操作： I. CHANGE u t : 把结点u的权值改为t II. QMAX u v: 询问从点u到点v的路径上的节点的最大权值 I
II. QSUM u v: 询问从点u到点v的路径上的节点的权值和 注意：从点u到点v的路径上的节点包括u和v本身

Input

　　输入的第一行为一个整数n，表示节点的个数。接下来n – 1行，每行2个整数a和b，表示节点a和节点b之间有
一条边相连。接下来n行，每行一个整数，第i行的整数wi表示节点i的权值。接下来1行，为一个整数q，表示操作
的总数。接下来q行，每行一个操作，以“CHANGE u t”或者“QMAX u v”或者“QSUM u v”的形式给出。 
对于100％的数据，保证1<=n<=30000，0<=q<=200000；中途操作中保证每个节点的权值w在-30000到30000之间。

Output

　　对于每个“QMAX”或者“QSUM”的操作，每行输出一个整数表示要求输出的结果。

Sample Input

4
1 2
2 3
4 1
4 2 1 3
12
QMAX 3 4
QMAX 3 3
QMAX 3 2
QMAX 2 3
QSUM 3 4
QSUM 2 1
CHANGE 1 5
QMAX 3 4
CHANGE 3 6
QMAX 3 4
QMAX 2 4
QSUM 3 4

Sample Output

4
1
2
2
10
6
5
6
5
16

这道题是树链剖分的模板题，下面是程序：

#include<stdio.h>
#include<iostream>
using namespace std;
const int N=30005;
struct edge{
	int v,next;
}e[N];
struct treenode{
	int lc,rc,mx,sum;
};
int max(int a,int b){
	return a>b?a:b;
}
struct Segment_tree{
	int k;
	treenode t[N<<2];
	Segment_tree(){
		k=1;
	}
	void build(int rt,int l,int r){
		if(l==r){
			t[rt].lc=t[rt].rc=0;
			return;
		}
		int m=l+r>>1;
		build(t[rt].lc=++k,l,m);
		build(t[rt].rc=++k,m+1,r);
		t[rt].mx=max(t[t[rt].lc].mx,t[t[rt].rc].mx);
		t[rt].sum=t[t[rt].lc].sum+t[t[rt].rc].sum;
	}
	void change(int rt,int l,int r,int &x,int &w){
		if(l==r){
			t[rt].sum=t[rt].mx=w;
			return;
		}
		int m=l+r>>1;
		if(x<=m){
			change(t[rt].lc,l,m,x,w);
		}
		else{
			change(t[rt].rc,m+1,r,x,w);
		}
		t[rt].mx=max(t[t[rt].lc].mx,t[t[rt].rc].mx);
		t[rt].sum=t[t[rt].lc].sum+t[t[rt].rc].sum;
	}
	int asksum(int rt,int l,int r,int a,int b){
		if(l==a&&r==b){
			return t[rt].sum;
		}
		int m=l+r>>1;
		if(b<=m){
			return asksum(t[rt].lc,l,m,a,b);
		}
		if(a>m){
			return asksum(t[rt].rc,m+1,r,a,b);
		}
		return asksum(t[rt].lc,l,m,a,m)+asksum(t[rt].rc,m+1,r,m+1,b);
	}
	int askmx(int rt,int l,int r,int a,int b){
		if(l==a&&r==b){
			return t[rt].mx;
		}
		int m=l+r>>1;
		if(b<=m){
			return askmx(t[rt].lc,l,m,a,b);
		}
		if(a>m){
			return askmx(t[rt].rc,m+1,r,a,b);
		}
		return max(askmx(t[rt].lc,l,m,a,m),askmx(t[rt].rc,m+1,r,m+1,b));
	}
}t;
int n,sum[N],son[N],top[N],dep[N],fa[N],dfn[N],rank[N],head[N],k,dfs_clock;
void add(int u,int v){
	++k;
	e[k].v=v;
	e[k].next=head[u];
	head[u]=k;
}
void dfs1(int u,int f,int d){
	dep[u]=d;
	sum[u]=1;
	fa[u]=f;
	son[u]=0;
	int i,v;
	for(i=head[u];i;i=e[i].next){
		v=e[i].v;
		if(v!=f){
			dfs1(v,u,d+1);
			sum[u]+=sum[v];
			if(sum[v]>sum[son[u]]){
				son[u]=v;
			}
		}
	}
}
void dfs2(int u,int t){
	top[u]=t;
	dfn[u]=++dfs_clock;
	rank[dfs_clock]=u;
	if(!son[u]){
		return;
	}
	dfs2(son[u],t);
	int i,v;
	for(i=head[u];i;i=e[i].next){
		v=e[i].v;
		if(v!=fa[u]&&v!=son[u]){
			dfs2(v,v);
		}
	}
}
int asksum(int x,int y){
	int ans=0;
	while(top[x]!=top[y]){
		if(dep[top[x]]<dep[top[y]]){
			swap(x,y);
		}
		ans+=t.asksum(1,1,n,dfn[top[x]],dfn[x]);
		x=fa[top[x]];
	}
	if(dfn[x]>dfn[y]){
		swap(x,y);
	}
	ans+=t.asksum(1,1,n,dfn[x],dfn[y]);
	return ans;
}
int askmx(int x,int y){
	int ans=-0x7fffffff;
	while(top[x]!=top[y]){
		if(dep[top[x]]<dep[top[y]]){
			swap(x,y);
		}
		ans=max(ans,t.askmx(1,1,n,dfn[top[x]],dfn[x]));
		x=fa[top[x]];
	}
	if(dfn[x]>dfn[y]){
		swap(x,y);
	}
	ans=max(ans,t.askmx(1,1,n,dfn[x],dfn[y]));
	return ans;
}
char ch(){
	char c=getchar();
	while(c!='U'&&c!='X'&&c!='E'){
		c=getchar();
	}
	return c;
}
void read(int &s){
	s=0;
	int f=1;
	char c=getchar();
	while((c<'0'||c>'9')&&c!='-'){
		c=getchar();
	}
	if(c=='-'){
		f=-1;
		c=getchar();
	}
	while(c>='0'&&c<='9'){
		s=s*10+c-'0';
		c=getchar();
	}
	s*=f;
}
int main(){
	int i,u,v,m;
	char c;
	scanf("%d",&n);
	for(i=1;i<n;i++){
		scanf("%d%d",&u,&v);
		add(u,v);
		add(v,u);
	}
	dfs1(1,-1,1);
	dfs2(1,1);
	t.build(1,1,n);
	for(i=1;i<=n;i++){
		scanf("%d",&m);
		t.change(1,1,n,dfn[i],m);
	}
	scanf("%d",&m);
	while(m--){
		c=ch();
		read(u);
		read(v);
		if(c=='X'){
			printf("%d\n",askmx(u,v));
			continue;
		}
		if(c=='U'){
			printf("%d\n",asksum(u,v));
			continue;
		}
		t.change(1,1,n,dfn[u],v);
	}
	return 0;
}