poj2763 树链剖分（线段树）

注意这里都是把边放到线段树中，所以lca的时候，要注意如果top[x]==top[y] && x==y 的时候已经完成了。

仔细想想边和点的不同之处！！！

#include<map>
#include<queue>
#include<stack>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define INF 1000000007
#define mod 100000000
#define ll long long
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
using namespace std;
const int MAXN = ;
struct node
{
    int to;
    int val;
    int next;
}edge[MAXN*];
int a[MAXN][];
int fa[MAXN],pre[MAXN],ind,w[MAXN],top[MAXN],siz[MAXN],son[MAXN],deq[MAXN],cnt;//w[]表示它和他父亲的连线在线段树种的位置
int vis[MAXN],n,m,tree[MAXN<<];
void add(int x,int y,int z)
{
    edge[ind].to = y;
    edge[ind].val = z;
    edge[ind].next = pre[x];
    pre[x] = ind++;
}
void dfs1(int rt,int pa,int d)
{
    vis[rt] = ;
    siz[rt] = ;
    son[rt] = -;
    deq[rt] = d;
    fa[rt] = pa;
    for(int i = pre[rt]; i != -; i = edge[i].next){
        int t = edge[i].to;
        if(!vis[t]){
            dfs1(t,rt,d+);
            siz[rt] += siz[t];
            if(siz[t] > siz[son[rt]]){
                son[rt] = t;
            }
        }
    }
}
void dfs2(int rt,int tp)
{
    vis[rt] = ;
    w[rt] = ++ cnt;
    top[rt] = tp;
    if(son[rt] != -){
        dfs2(son[rt],tp);
    }
    for(int i = pre[rt]; i != -; i = edge[i].next){
        int t = edge[i].to;
        if(!vis[t] && t != son[rt]){
            dfs2(t,t);
        }
    }
}
void pushup(int rt)
{
    tree[rt] = tree[rt<<] + tree[rt<<|];
}
void updata(int p,int val,int l,int r,int rt)
{
    if(l == r){
        tree[rt] = val;
        return ;
    }
    int m = (l + r) >> ;
    if(m >= p){
        updata(p,val,lson);
    }
    else {
        updata(p,val,rson);
    }
    pushup(rt);
}
int query(int L,int R,int l,int r,int rt)
{
    if(L <= l && r <= R){
        return tree[rt];
    }
    int m = (l + r) >> ;
    int ans = ;
    if(m >= L){
        ans += query(L,R,lson);
    }
    if(m < R){
        ans += query(L,R,rson);
    }
    return ans;
}
int lca(int x,int y)
{
    int ans = ;
    while(top[x] != top[y]){
        if(deq[top[x]] < deq[top[y]]){
            swap(x,y);
        }
        ans += query(w[top[x]],w[x],,cnt,);
        x = fa[top[x]];
    }
    if(x == y){//同一点没有边，仔细想想这就是为什么边和点的不同地方！
        return ans;
    }
    if(deq[x] < deq[y]){
        swap(x,y);
    }

    ans += query(w[son[y]],w[x],,cnt,);//既然top[x] == top[y]说明x 和 y在同一条重边上 所以son[y]也在该路径上
                                         //并且是树根的孩子，又w[]表示该点和其父亲相连的点的边在线段树中的编号，
                                         //所以需要这步，这样能得出答案。
    return ans;
}
int main()
{
    int s,p;
    while(~scanf("%d%d%d",&n,&p,&s)){
        int x,y,z;
        ind = ;
        memset(pre,-,sizeof(pre));
        for(int i = ; i < n; i++){
            scanf("%d%d%d",&a[i][],&a[i][],&a[i][]);
            add(a[i][],a[i][],a[i][]);
            add(a[i][],a[i][],a[i][]);
        }
        memset(vis,,sizeof(vis));
        dfs1(,,);//求deq siz son fa

        cnt = ;
        memset(vis,,sizeof(vis));
        dfs2(,);//求top w
        memset(tree,,sizeof(tree));
        for(int i = ; i < n; i++){
            if(deq[a[i][]] < deq[a[i][]]){//让a[i][0]在下面
                swap(a[i][],a[i][]);
            }
            updata(w[a[i][]],a[i][],,cnt,);//将a[i][0]和他父亲相连的边插入到线段树中
        }
        int q;
        while(p--){
            scanf("%d",&q);
            if(q == ){
                scanf("%d",&x);
                printf("%d\n",lca(s,x));
                s = x;
            }
            else {
                scanf("%d%d",&x,&y);
                updata(w[a[x][]],y,,cnt,);
            }
        }
    }
    return ;
}