hdu_5029_relief grain(树链剖分)

题目连接：http://acm.hdu.edu.cn/showproblem.php?pid=5029

题意：给你一个树，然后给你两点，将这两点之间的点涂上颜色，问涂色最多的那个颜色是什么，如果数量相同，就输出编号小的颜色

题解：树链剖分，这个题很巧妙，线段树维护的是一个颜色线段树，然后如果要涂[l,r]的颜色a，就在l处标记a +1，在r+1处标记a  -1。然后从左往右依次更新线段树，并query，就能得出答案，我的跑了1400+ms,鸟神的非递归线段树只跑我一半的时间

详解请看：鸟神的这题题解：relief grain

#include<cstdio>
#include<cstring>
#define F(i,a,b) for(int i=a;i<=b;++i)
#define root 0,100001,1
#define ls l,m,rt<<1
#define rs m+1,r,rt<<1|1
#pragma comment(linker, "/STACK:102400000,102400000")
const int N=100010;
int n,m,x,y,z,idx,ans[N],sum[N*4],nxt[10*N],g[N],v[N*10],w[N*10],ed,pre[N],siz[N],dep[N],hson[N],top[N],tid[N],ftid[N];

inline void adg(int x,int y,int z){v[++ed]=y,w[ed]=z,nxt[ed]=g[x],g[x]=ed;}
//树链部分
void dfs1(int u,int fa){
    pre[u]=fa,siz[u]=1,dep[u]=dep[fa]+1,hson[u]=0;
    for(int i=g[u];~i;i=nxt[i]){
        int vv=v[i];
        if(vv!=fa){
            dfs1(vv,u);
            if(siz[vv]>siz[hson[u]])hson[u]=vv;
            siz[u]+=siz[vv];
        }
    }
}

void dfs2(int u,int tp){
    tid[u]=idx,ftid[idx++]=u,top[u]=tp;
    if(hson[u])dfs2(hson[u],tp);
    for(int i=g[u];~i;i=nxt[i]){
        int vv=v[i];
        if(vv!=pre[u]&&vv!=hson[u])dfs2(vv,vv);
    }
}

void Update(int u,int v,int c){
    int fu=top[u],fv=top[v];
    while(fu!=fv){
        if(dep[fu]>=dep[fv])
            adg(tid[fu],c,1),adg(tid[u]+1,c,-1),u=pre[fu],fu=top[u];
        else adg(tid[fv],c,1),adg(tid[v]+1,c,-1),v=pre[fv],fv=top[v];
    }
    if(dep[u]>dep[v])u=u^v,v=u^v,u=u^v;
    adg(tid[u],c,1),adg(tid[v]+1,c,-1);
}
//线段树部分
void update(int x,int v,int l,int r,int rt){
    if(l==r){sum[rt]+=v;return;}
    int m=(l+r)>>1;
    if(x<=m)update(x,v,ls);else update(x,v,rs);
    sum[rt]=sum[rt<<1]>sum[rt<<1|1]?sum[rt<<1]:sum[rt<<1|1];
}

int query(int l,int r,int rt){
    if(l==r)return l;
    int m=(l+r)>>1;
    if(sum[rt]==sum[rt<<1])return query(ls);
    else return query(rs);
}

int main(){
    while(~scanf("%d%d",&n,&m),n+m){
        memset(g,-1,sizeof(g)),ed=0;
        F(i,1,n-1)scanf("%d%d",&x,&y),adg(x,y,0),adg(y,x,0);
        dfs1(1,0),idx=1,dfs2(1,1);
        F(i,0,4*N-1)sum[i]=0;
        memset(g,-1,sizeof(g)),ed=0;
        F(i,1,m)scanf("%d%d%d",&x,&y,&z),Update(x,y,z);
        F(u,1,n){
            for(int i=g[u];~i;i=nxt[i])update(v[i],w[i],root);
            ans[ftid[u]]=query(root);
        }
        F(i,1,n)printf("%d
",ans[i]);
    }
    return 0;
}