树链剖分学习

啊，好无聊所以又爬了一下午的树。

这次是树链剖分，讲道理，从例题出发，SPOJ QTREE。

题意就是动态修改变权值，求树上两点之间的所有边权的最大值。

我的体会：

其实就是在树上dfs一遍，记录下每个节点的dfs序，只不过这里的dfs是优先去选最长的链的（即那个儿子节点的子树节点数目越多就去访问谁）。代码如下：

void buildTree(int v,int tp)
{
    T[v] = ++ord, top[v] = tp;
    // cout << v << " " << T[v] << " " << top[v] << endl;
    if(son[v]) {
        buildTree(son[v],top[v]);
    }
    for(int i = head[v];i != -1;i = e[i].next) {
        int u = e[i].to;
        if(u != fa[v] && u != son[v]) {
            buildTree(u,u);
        }
    }
}

每条点维护有top节点，即该点所在的链的头结点，然后这些链通过dfs序就可以映射到线段树上面从而维护一些信息。然后既然知道了每个点在线段树的位置，以及每个点所在的链的头，那么其余的操作就是类似于线段树啦，区间查询更新之类的。。。

求两点之间的最大值：

int findMax(int u,int v)
{
    //T数组存的是节点在线段树的位置，top数组存的是链的头结点编号
    //dep数组存的是节点的深度
    int f1 = top[u], f2 = top[v], tmp = 0;
    while (f1 != f2)
    {
        //**这里是比较top的深度
        if (dep[f1] < dep[f2]) {
            swap(f1, f2); 
            swap(u, v);
        }
        tmp = max(tmp, query(1, 1, ord, T[f1], T[u]));
        u = fa[f1]; f1 = top[u];
    }
    if (u == v) return tmp;
    if (dep[u] > dep[v]) swap(u, v);
    return max(tmp, query(1, 1, ord, T[son[u]], T[v]));
}

还是很好懂的哈，有问题的可以戳 这里