2019CCPC女生赛 D.Tree题解【树链剖分+线段树】

Sample Input
4 4
2 3 4 5
1 2
2 3
2 4
0 3 4
0 1 3
1 2 3
1 1 4

Sample Output
2
4

我们模拟赛时还不会树链剖分，用dfs稳的超时。于是赛后学习了一下树剖，并且摸清了线段树维护除了区间加、区间乘，输出区间和以外的线段树维护方法。

正解：树链剖分

通过给树的结点重新编号，使得树中任意路径转化成logn段连续的区间
转换为区间后用线段树进行维护
重儿子：结点最多的子树，若相同则任选一个
轻儿子：剩下的所有结点
重边：重儿子和父结点之间的边
重链：重边构成的重链
重儿子在重链里面，轻儿子在以这个轻儿子为开头往下的重链里面
定理：树中任意一条路径均可拆分成Ologn个连续区间（重链）。
按照树的dfs序来将树转化为序列，要求优先遍历当前结点的重儿子。这样能使得每一条重链上的点的编号连续。
求两个点之间的距离，类似找最近公共祖先。就是找到自己所在的重链的开头往上寻找，一直沿着重链走，直到到达同一条重链，将两点之间的部分抠出来，得到logn条重链。得到logn个区间。
访问每个区间logn,总的时间复杂度为(logn)2
线段树维护开方后区间和

#include <iostream>
#include <cstring>
#include <algorithm>

using namespace std;

typedef long long LL;
const int N = 100010, M = N * 2;

int n, m;
int w[N], h[N], e[M], ne[M], idx;
int id[N], nw[N], cnt;
int dep[N], sz[N], top[N], fa[N], son[N];
struct Tree
{
    int l, r;
    LL add, sum;
}tr[N * 4];

void add(int a, int b)
{
    e[idx] = b, ne[idx] = h[a], h[a] = idx ++ ;
}

void dfs1(int u, int father, int depth)
{
    dep[u] = depth, fa[u] = father, sz[u] = 1;
    for (int i = h[u]; ~i; i = ne[i])
    {
        int j = e[i];
        if (j == father) continue;
        dfs1(j, u, depth + 1);
        sz[u] += sz[j];
        if (sz[son[u]] < sz[j]) son[u] = j;
    }
}

void dfs2(int u, int t)
{
    id[u] = ++ cnt, nw[cnt] = w[u], top[u] = t;
    if (!son[u]) return;
    dfs2(son[u], t);
    for (int i = h[u]; ~i; i = ne[i])
    {
        int j = e[i];
        if (j == fa[u] || j == son[u]) continue;
        dfs2(j, j);
    }
}

void pushup(int u)
{
    tr[u].sum = tr[u << 1].sum + tr[u << 1 | 1].sum;
}

void pushdown(int u)
{
    auto &root = tr[u], &left = tr[u << 1], &right = tr[u << 1 | 1];
    if (root.add)
    {
        left.add += root.add, left.sum += root.add * (left.r - left.l + 1);
        right.add += root.add, right.sum += root.add * (right.r - right.l + 1);
        root.add = 0;
    }
}

void build(int u, int l, int r)
{
    tr[u] = {l, r, 0, nw[r]};
    if (l == r) return;
    int mid = l + r >> 1;
    build(u << 1, l, mid), build(u << 1 | 1, mid + 1, r);
    pushup(u);
}

void update(int u, int l, int r)
{
    if (tr[u].l == tr[u].r)//当区间长度为1时只对tr[u].sum修改
    {
        tr[u].sum = sqrt(tr[u].sum);
        return;
    }
    pushdown(u);
    int mid = tr[u].l + tr[u].r >> 1;
    if (l <= mid) update(u << 1, l, r);
    if (r > mid) update(u << 1 | 1, l, r);
    pushup(u);
}

LL query(int u, int l, int r)
{
    if (l <= tr[u].l && r >= tr[u].r) return tr[u].sum;
    pushdown(u);
    int mid = tr[u].l + tr[u].r >> 1;
    LL res = 0;
    if (l <= mid) res += query(u << 1, l, r);
    if (r > mid) res += query(u << 1 | 1, l, r);
    return res;
}

void update_path(int u, int v)
{
    while (top[u] != top[v])
    {
        if (dep[top[u]] < dep[top[v]]) swap(u, v);
        update(1, id[top[u]], id[u]);
        u = fa[top[u]];
    }
    if (dep[u] < dep[v]) swap(u, v);
    update(1, id[v], id[u]);
}

LL query_path(int u, int v)
{
    LL res = 0;
    while (top[u] != top[v])
    {
        if (dep[top[u]] < dep[top[v]]) swap(u, v);
        res += query(1, id[top[u]], id[u]);
        u = fa[top[u]];
    }
    if (dep[u] < dep[v]) swap(u, v);
    res += query(1, id[v], id[u]);
    return res;
}

int main()
{
    scanf("%d%d", &n,&m);
    for (int i = 1; i <= n; i ++ ) scanf("%d", &w[i]);
    memset(h, -1, sizeof h);
    for (int i = 0; i < n - 1; i ++ )
    {
        int a, b;
        scanf("%d%d", &a, &b);
        add(a, b), add(b, a);
    }
    dfs1(1, -1, 1);
    dfs2(1, 1);
    build(1, 1, n);

    while (m -- )
    {
        int t, u, v, k;
        scanf("%d%d", &t, &u);
        if (t == 0)
        {
            scanf("%d", &v);
            update_path(u, v);
        }
        else if (t == 1)
        {
            scanf("%d", &v);
            printf("%lld\n", query_path(u, v));
        }
    }

    return 0;
}