[SHOI 2012] 魔法树

[题目链接]

https://www.lydsy.com/JudgeOnline/problem.php?id=2836

[算法]

树链剖分

时间复杂度 : O(NlogN ^ 2)

[代码]

#include<bits/stdc++.h>
using namespace std;
#define MAXN 100010
typedef long long LL;

struct edge
{
    int to , nxt;
} e[MAXN << ];

int n , tot , timer;
int head[MAXN] , size[MAXN] , son[MAXN]  , fa[MAXN] , dfn[MAXN] , depth[MAXN] , top[MAXN];

template <typename T> inline void chkmax(T &x , T y) { x = max(x , y); }
template <typename T> inline void chkmin(T &x , T y) { x = min(x , y); }
template <typename T> inline void read(T &x)
{
    T f = ; x = ;
    char c = getchar();
    for (; !isdigit(c); c = getchar()) if (c == '-') f = -f;
    for (; isdigit(c); c = getchar()) x = (x << ) + (x << ) + c - '';
    x *= f;
}
struct Segment_Tree
{
    struct Node
    {
        int l , r;
        LL sum , tag;
    } Tree[MAXN << ];
    inline void build(int index , int l , int r)
    {
        Tree[index].l = l;
        Tree[index].r = r;
        Tree[index].tag = Tree[index].sum = ;
        if (l == r) return;
        int mid = (l + r) >> ;
        build(index <<  , l , mid);
        build(index <<  |  , mid +  , r);
    }
    inline void pushdown(int index)
    {
        int l = Tree[index].l , r = Tree[index].r;
        int mid = (l + r) >> ;
        Tree[index << ].sum += (mid - l + ) * Tree[index].tag;
        Tree[index <<  | ].sum += (r - mid) * Tree[index].tag;
        Tree[index << ].tag += Tree[index].tag;
        Tree[index <<  | ].tag += Tree[index].tag;
        Tree[index].tag = ;
    }
    inline void update(int index)
    {
        Tree[index].sum = Tree[index << ].sum + Tree[index <<  | ].sum;
    }
    inline void modify(int index , int l , int r , LL value)
    {
        if (Tree[index].l == l && Tree[index].r == r)
        {
            Tree[index].sum += value * (r - l + );
            Tree[index].tag += value;
            return;
        }
        pushdown(index);
        int mid = (Tree[index].l + Tree[index].r) >> ;
        if (mid >= r) modify(index <<  , l , r , value);
        else if (mid +  <= l) modify(index <<  |  , l , r , value);
        else
        {
            modify(index <<  , l , mid , value);
            modify(index <<  |  , mid +  , r , value);
        }
        update(index);
    }
    inline LL query(int index , int l , int r)
    {
        if (Tree[index].l == l && Tree[index].r == r) return Tree[index].sum;
        pushdown(index);
        int mid = (Tree[index].l + Tree[index].r) >> ;
        if (mid >= r) return query(index <<  , l , r);
        else if (mid +  <= l) return query(index <<  |  , l , r);
        else return query(index <<  , l , mid) + query(index <<  |  , mid +  , r);
    }
} SGT;
inline void addedge(int u , int v)
{
    ++tot;
    e[tot] = (edge){v , head[u]};
    head[u] = tot;
}
inline void dfs1(int u)
{
    son[u] = -;
    size[u] = ;
    for (int i = head[u]; i; i = e[i].nxt)
    {
        int v = e[i].to;
        if (v == fa[u]) continue;
        fa[v] = u;
        depth[v] = depth[u] + ;
        dfs1(v);
        size[u] += size[v];
        if (son[u] == - || size[v] > size[son[u]]) son[u] = v;
    }
}
inline void dfs2(int u , int tp)
{
    dfn[u] = ++timer;
    top[u] = tp;
    if (son[u] != -) dfs2(son[u] , tp);
    for (int i = head[u]; i; i = e[i].nxt)
    {
        int v = e[i].to;
        if (v == fa[u] || v == son[u]) continue;
        dfs2(v , v);
    }
}
inline void modify(int u , int v , LL d)
{
    int tu = top[u] , tv = top[v];
    while (tu != tv)
    {
        if (depth[tu] > depth[tv])
        {
            swap(u , v);
            swap(tu , tv);
        }
        SGT.modify( , dfn[tv] , dfn[v] , d);
        v = fa[tv]; tv = top[v];
    }
    if (depth[u] > depth[v]) swap(u , v);
    SGT.modify(, dfn[u] , dfn[v] , d);
}

int main()
{

    scanf("%d" , &n);
    for (int i = ; i < n; i++)
    {
        int u , v;
        scanf("%d%d" , &u , &v);
        addedge(u , v);
        fa[v] = u;
    }
    dfs1();
    dfs2( , );
    SGT.build( ,  , n);
    int q;
    scanf("%d" , &q);
    while (q--)
    {
        char op[];
        scanf("%s" , &op);
        if (op[] == 'A')
        {
            int u , v;
            LL d;
            scanf("%d%d%lld" , &u , &v , &d);
            modify(u , v , d);
        } else
        {
            int u;
            scanf("%d" , &u);
            printf("%lld\n" , SGT.query( , dfn[u] , dfn[u] + size[u] - ));
        }
    }

    return ;
}