HDU 3966 Aragorn's Story (树链剖分)

题目地址：HDU 3966
树链剖分模板题。不多说、、
代码如下：

#include <iostream>
#include <string.h>
#include <math.h>
#include <queue>
#include <algorithm>
#include <stdlib.h>
#include <map>
#include <set>
#include <stdio.h>
using namespace std;
#define LL long long
#define pi acos(-1.0)
#pragma comment(linker, "/STACK:1024000000,1024000000")
const int mod=1e9+7;
const int INF=0x3f3f3f3f;
const double eqs=1e-3;
const int MAXN=50000+10;
#define root 1, n, 1
#define lson l, mid, rt<<1
#define rson mid+1, r, rt<<1|1
int head[MAXN], cnt, a[MAXN], n;
int dep[MAXN], siz[MAXN], son[MAXN], w[MAXN], top[MAXN], tot, fa[MAXN];
int sum[MAXN<<2], lazy[MAXN<<2];
struct node
{
        int u, v, next;
}edge[MAXN<<1];
void add(int u, int v)
{
        edge[cnt].v=v;
        edge[cnt].u=u;
        edge[cnt].next=head[u];
        head[u]=cnt++;
}
void init()
{
        memset(head,-1,sizeof(head));
        cnt=0;
        memset(dep,0,sizeof(dep));
        memset(son,0,sizeof(son));
        tot=0;
        memset(sum,0,sizeof(sum));
        memset(lazy,0,sizeof(lazy));
}
void dfs1(int u, int p)
{
        siz[u]=1;
        for(int i=head[u];i!=-1;i=edge[i].next){
                int v=edge[i].v;
                if(v==p) continue ;
                fa[v]=u;
                dep[v]=dep[u]+1;
                dfs1(v,u);
                siz[u]+=siz[v];
                if(siz[son[u]]<siz[v]) son[u]=v;
        }
}
void dfs2(int u, int tp)
{
        top[u]=tp;w[u]=++tot;
        if(son[u]) dfs2(son[u],tp);
        for(int i=head[u];i!=-1;i=edge[i].next){
                int v=edge[i].v;
                if(v!=fa[u]&&v!=son[u])
                        dfs2(v,v);
        }
}
struct Line_Tree
{
        void PushUp(int rt)
        {
                sum[rt]=sum[rt<<1]+sum[rt<<1|1];
        }
        void PushDown(int l, int r, int rt)
        {
                int mid=l+r>>1;
                if(lazy[rt]){
                        sum[rt<<1]+=(mid-l+1)*lazy[rt];
                        sum[rt<<1|1]+=(r-mid)*lazy[rt];
                        lazy[rt<<1]+=lazy[rt];
                        lazy[rt<<1|1]+=lazy[rt];
                        lazy[rt]=0;
                }
        }
        void Update(int ll, int rr, int x, int l, int r, int rt)
        {
                if(ll<=l&&rr>=r){
                        lazy[rt]+=x;
                        sum[rt]+=x*(r-l+1);
                        return ;
                }
                int mid=l+r>>1;
                PushDown(l,r,rt);
                if(ll<=mid) Update(ll,rr,x,lson);
                if(rr>mid) Update(ll,rr,x,rson);
                PushUp(rt);
        }
        int Query(int p, int l, int r, int rt)
        {
                if(l==r){
                        return sum[rt];
                }
                int mid=l+r>>1;
                PushDown(l,r,rt);
                if(p<=mid) return Query(p,lson);
                return Query(p,rson);
        }
}lt;
void Update(int u, int v, int x)
{
        int f1=top[u], f2=top[v];
        while(f1!=f2){
                if(dep[f1]<dep[f2]){
                        swap(u,v);
                        swap(f1,f2);
                }
                lt.Update(w[f1],w[u],x,root);
                u=fa[f1];f1=top[u];
        }
        if(dep[u]<dep[v]){
                swap(u,v);
        }
        lt.Update(w[v],w[u],x,root);
}
int main()
{
        int m, q, i, j, u, v, x;
        char s[10];
        while(scanf("%d%d%d",&n,&m,&q)!=EOF){
                for(i=1;i<=n;i++){
                        scanf("%d",&a[i]);
                }
                init();
                while(m--){
                        scanf("%d%d",&u,&v);
                        add(u,v);
                        add(v,u);
                }
                dfs1(1,-1);
                dfs2(1,1);
                for(i=1;i<=n;i++){
                        lt.Update(w[i],w[i],a[i],root);
                }
                while(q--){
                        scanf("%s",s);
                        if(s[0]=='I'){
                                scanf("%d%d%d",&u,&v,&x);
                                Update(u,v,x);
                        }
                        else if(s[0]=='D'){
                                scanf("%d%d%d",&u,&v,&x);
                                Update(u,v,-x);
                        }
                        else{
                                scanf("%d",&x);
                                printf("%d\n",lt.Query(w[x],root));
                        }
                }
        }
        return 0;
}