[BZOJ 1036][ZJOI2008]樹的統計Count（樹鏈剖分）

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

Solution

---

發現還是模板_(:з」∠)_

#include<iostream>
#include<cstring>
#include<cstdio>
#define MAXN 30005
#define inf 0x3f3f3f3f
#define Max(a,b) (a>b?a:b)
using namespace std;
int n,q,cnt=0,head[MAXN],sz=0,ans;
int siz[MAXN],father[MAXN],dep[MAXN],top[MAXN],pos[MAXN];
struct Node1{
    int next,to;
}Edges[MAXN*2];
void addedge(int u,int v)
{
    Edges[++cnt].next=head[u];
    head[u]=cnt;
    Edges[cnt].to=v;
}
void dfs1(int u)
{
    siz[u]=1;
    for(int i=head[u];~i;i=Edges[i].next)
    {
        int v=Edges[i].to;
        if(v==father[u])continue;
        father[v]=u;dep[v]=dep[u]+1;
        dfs1(v);
        siz[u]+=siz[v];
    }
}
void dfs2(int u,int t)
{
    sz++;pos[u]=sz;
    top[u]=t;
    int k,maxsiz=0;
    for(int i=head[u];~i;i=Edges[i].next)
    {
        int v=Edges[i].to;
        if(v!=father[u]&&siz[v]>maxsiz)
        maxsiz=siz[v],k=v;
    }
    if(!maxsiz)return;
    dfs2(k,t);
    for(int i=head[u];~i;i=Edges[i].next)
    {
        int v=Edges[i].to;
        if(v!=father[u]&&v!=k)
        dfs2(v,v);
    }
}
struct Node2{
    int l,r,sum,maxn;
}t[MAXN*4];
void build(int idx,int l,int r)
{
    t[idx].l=l,t[idx].r=r;
    if(l==r)return;
    int mid=(l+r)/2;
    build(idx*2,l,mid);
    build(idx*2+1,mid+1,r);
}
void change(int idx,int a,int b)
{
    if(t[idx].l==t[idx].r)
    {
        t[idx].maxn=t[idx].sum=b;
        return;
    }
    int mid=(t[idx].l+t[idx].r)/2;
    if(a<=mid)change(idx*2,a,b);
    else change(idx*2+1,a,b);
    t[idx].maxn=Max(t[idx*2].maxn,t[idx*2+1].maxn);
    t[idx].sum=t[idx*2].sum+t[idx*2+1].sum;
}
void qmax(int idx,int a,int b)
{
    if(a<=t[idx].l&&b>=t[idx].r)
    {
        ans=Max(ans,t[idx].maxn);
        return;
    }
    int mid=(t[idx].l+t[idx].r)/2;
    if(b<=mid)qmax(idx*2,a,b);
    else if(a>mid)qmax(idx*2+1,a,b);
    else
    {
        qmax(idx*2,a,b);
        qmax(idx*2+1,a,b);
    }
}
void qsum(int idx,int a,int b)
{
    if(a<=t[idx].l&&b>=t[idx].r)
    {
        ans+=t[idx].sum;
        return;
    }
    int mid=(t[idx].l+t[idx].r)/2;
    if(b<=mid)qsum(idx*2,a,b);
    else if(a>mid)qsum(idx*2+1,a,b);
    else
    {
        qsum(idx*2,a,b);
        qsum(idx*2+1,a,b);
    }
}
int main()
{
    memset(head,-1,sizeof(head));
    scanf("%d",&n);
    int a,b;
    for(int i=1;i<n;i++)
    {
        scanf("%d%d",&a,&b);
        addedge(a,b);
        addedge(b,a);
    }
    dfs1(1);dfs2(1,1);
    build(1,1,sz);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&a);
        change(1,pos[i],a);
    }
    scanf("%d",&q);
    char opt[15];
    for(int i=1;i<=q;i++)
    {
        scanf("%s%d%d",&opt,&a,&b);
        if(!memcmp(opt,"CHANGE",strlen(opt)))
        {
            change(1,pos[a],b);
        }
        else if(!memcmp(opt,"QMAX",strlen(opt)))
        {
            ans=-inf;
            while(top[a]!=top[b])
            {
                if(dep[top[a]]<dep[top[b]])swap(a,b);
                qmax(1,pos[top[a]],pos[a]);
                a=father[top[a]];
            }
            if(dep[a]>dep[b])swap(a,b);
            qmax(1,pos[a],pos[b]);
            printf("%d\n",ans);
        }
        else if(!memcmp(opt,"QSUM",strlen(opt)))
        {
            ans=0;
            while(top[a]!=top[b])
            {
                if(dep[top[a]]<dep[top[b]])swap(a,b);
                qsum(1,pos[top[a]],pos[a]);
                a=father[top[a]];
            }
            if(dep[a]>dep[b])swap(a,b);
            qsum(1,pos[a],pos[b]);
            printf("%d\n",ans);
        }
    }
    return 0;
}