树的统计 ZJOI2008

树的统计 zjoi2008 —树链剖分

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

输入：

输入文件的第一行为一个整数n，表示节点的个数。
接下来n – 1行，每行2个整数a和b，表示节点a和节点b之间有一条边相连。
接下来一行n个整数，第i个整数wi表示节点i的权值。
接下来1行，为一个整数q，表示操作的总数。
接下来q行，每行一个操作，以“CHANGE u t”或者“QMAX u v”或者“QSUM u v”的形式给出。

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

样例输入：
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

样例输出：
4
1
2
2
10
6
5
6
5
16

数据范围：
对于100％的数据，保证1<=n<=30000，0<=q<=200000；中途操作中保证每个节点的权值w在-30000到30000之间。

---

解题思路：
我也是够了，一道树链剖分的模板题用了我整整一个晚上，智商不够啊！而且还是被一个SB错误卡了下去。。。言归正传，正宗的树链剖分模板题，线段树上维护即可，判断轻重子节点拉成一条链，用num数组很巧妙地在树链和线段树之间转换，非常玄学的操作，废话不多说，上代码。

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
#define INF 0x3f3f3f3f
#define LD (o<<1)
#define RD ((o<<1)|1)
#define mem(a) memset(a,0,sizeof(a))
#define N 30010
struct Tree{int l,r,maxv,sum;}T[4*N];//线段树
struct Edge{int v,next;}E[N];//边
struct Node{
    int dep,siz,fa,w,hson,num,top;
    Node(){dep=siz=fa=w=hson=num=top=0;}
}P[N];//点
/*
dep深度
size节点数
fa父亲
w节点权值
hson重孩子
num在线段树编号
top表示每条链第一个位置
*/
int first[N],pred[N];
/*
first 邻接表头指针
pred线段树的编号是原来什么编号
*/
int n,m,last,tot;//tot记录线段树总编号
char c[10];
void init(){
    last=tot=0;
    mem(first);mem(pred);
}
void add(int u,int v){
    E[++last].v=v;
    E[last].next=first[u];
    first[u]=last;
}
void dfs1(int s){//深搜出每个节点的重子节点，顺便建树
    P[s].siz=1;//当前为根节点的子树大小
    P[s].hson=0;//默认重孩子为0
    for(int i=first[s];i;i=E[i].next){//邻接表
        int &son=E[i].v;
        if(son==P[s].fa)continue;
        P[son].dep=P[s].dep+1;
        P[son].fa=s;
        dfs1(son);
        P[s].siz+=P[son].siz;
        if(P[P[s].hson].siz<P[son].siz)//打擂台判断重子节点
            P[s].hson=son;
    }
}
void dfs2(int s,int tp){//深搜出重链
    P[s].top=tp;//每一条链的起点
    P[s].num=++tot;//重链节点编号
    pred[tot]=s;
    if(P[s].hson)dfs2(P[s].hson,tp);//有重孩子继续深搜
    for(int i=first[s];i;i=E[i].next){
        int &son=E[i].v;
        if(son!=P[s].fa&&son!=P[s].hson)
            dfs2(son,son);//重新搜索一链
    }
}
void pushup(int o){
    T[o].sum=T[LD].sum+T[RD].sum;
    T[o].maxv=max(T[LD].maxv,T[RD].maxv);
}
void build(int o,int l,int r){
    T[o].l=l;T[o].r=r;
    if(l==r){T[o].sum=T[o].maxv=P[pred[l]].w;return;}
    int mid=(l+r)>>1;
    build(LD,l,mid);
    build(RD,mid+1,r);
    pushup(o);
}
void update(int o,int p,int val){
    if(T[o].l==p&&T[o].r==p){T[o].sum=T[o].maxv=val;return;}
    int mid=(T[o].l+T[o].r)>>1;
    if(p<=mid)update(LD,p,val);
    if(p>=mid+1)update(RD,p,val);
    pushup(o);
}
int query_sum(int o,int l,int r){
    if(l<=T[o].l&&T[o].r<=r)return T[o].sum;
    int mid=(T[o].l+T[o].r)>>1;
    if(r<=mid)return query_sum(LD,l,r);
    if(l>mid)return query_sum(RD,l,r);
    if(l<=mid&&r>mid)return query_sum(LD,l,mid)+query_sum(RD,mid+1,r);
}
int query_maxv(int o,int l,int r){
    if(l<=T[o].l&&T[o].r<=r)return T[o].maxv;
    int mid=(T[o].l+T[o].r)>>1;
    if(r<=mid)return query_maxv(LD,l,r);
    if(l>mid)return query_maxv(RD,l,r);
    if(l<=mid&&r>mid)return max(query_maxv(LD,l,mid),query_maxv(RD,mid+1,r));
}
//求链 目的是想办法跳到一条重链上来
//由于一颗树只有一个完整重莲，如果2个点在一颗子树上，那么先统计轻边，自然会转到重边上来。
int get_sum(int x,int y){
    int ans=0;
    while(P[x].top!=P[y].top){//不是一条重链
        if(P[P[x].top].dep<P[P[y].top].dep)swap(x,y);
        //链的首点深度大的在前
        ans+=query_sum(1,P[P[x].top].num,P[x].num);
        //统计x这条链
        x=P[P[x].top].fa;//x变为这条链顶点的父亲
    }
    if(P[x].num>P[y].num)swap(x,y);
    ans+=query_sum(1,P[x].num,P[y].num);
    return ans;
}
int get_maxv(int x,int y){
    int ans=-INF;
    while(P[x].top!=P[y].top){//不是一条重链
        if(P[P[x].top].dep<P[P[y].top].dep)swap(x,y);
        //链的首点深度大的在前
        ans=max(ans,query_maxv(1,P[P[x].top].num,P[x].num));
        //统计x这条链
        x=P[P[x].top].fa;//x变为这条链顶点的父亲
    }
    if(P[x].num>P[y].num)swap(x,y);
    ans=max(ans,query_maxv(1,P[x].num,P[y].num));
    return ans;
}
int main(){
    init();
    scanf("%d",&n);
    for(int i=1;i<n;i++){
        int a,b;scanf("%d%d",&a,&b);
        add(a,b);
        add(b,a);//双向建边
    }
    for(int i=1;i<=n;i++)scanf("%d",&P[i].w);
    dfs1(1);
    dfs2(1,1);
    build(1,1,tot);//构建线段树
    scanf("%d",&m);
    for(int i=1;i<=m;i++){
        int x,y;scanf("%s%d%d",c,&x,&y);
        if(c[0]=='C'){P[x].w=y;update(1,P[x].num,y);}
        else{
            if(c[1]=='M')printf("%d\n",get_maxv(x,y));
            if(c[1]=='S')printf("%d\n",get_sum(x,y));
        }
    }
    return 0;
}