[HDU4123]Bob’s Race

题目大意：
给定一棵$n$个点并且有边权的树，每个点的权值为该点能走的最远长度，并输入$m$个询问，每次询问最多有多少个编号连续的点，他们的最大最小点权差小于等于$Q$。

思路：
两趟DP（DFS）求出每个点能走的最远长度，然后用ST算法预处理出每一段最大最小值。
对于每组询问，用尺取法求出最大值。
注意log2不能用cmath里面的函数（尤其是不能在for语句上），不然会TLE，最好是自己实现。

#include<cstdio>
#include<cctype>
#include<vector>
inline int getint() {
    char ch;
    while(!isdigit(ch=getchar()));
    int x=ch^'0';
    while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
    return x;
}
inline int log2(const float x){
    return ((unsigned&)x>>23&255)-127;
}
const int N=50001,logN=16;
struct Edge {
    int to,w;
};
std::vector<Edge> e[N];
inline void add_edge(const int u,const int v,const int w) {
    e[u].push_back((Edge){v,w});
}
int w[N][3],son[N];
void dfs1(const int x,const int par) {
    w[x][0]=w[x][1]=0;
    for(unsigned i=0;i<e[x].size();i++) {
        int &y=e[x][i].to;
        if(y==par) continue;
        dfs1(y,x);
        if(w[y][0]+e[x][i].w>w[x][0]) {
            w[x][1]=w[x][0];
            w[x][0]=w[y][0]+e[x][i].w;
            son[x]=y;
        }
        else if(w[y][0]+e[x][i].w>w[x][1]) {
            w[x][1]=w[y][0]+e[x][i].w;
        }
    }
}
void dfs2(const int x,const int par) {
    for(unsigned i=0;i<e[x].size();i++) {
        int &y=e[x][i].to;
        if(y==par) continue;
        w[y][2]=std::max(w[x][2],son[x]!=y?w[x][0]:w[x][1])+e[x][i].w;
        dfs2(y,x);
    }
}
int max[N][logN],min[N][logN];
inline int getsum(const int l,const int r) {
    int k=log2(r-l+1);
    return std::max(max[l][k],max[r-(1<<k)+1][k])-std::min(min[l][k],min[r-(1<<k)+1][k]);
}
int main() {
    for(;;) {
        int n=getint(),m=getint();
        if(!n&&!m) return 0;
        for(int i=1;i<n;i++) {
            int u=getint(),v=getint(),w=getint();
            add_edge(u,v,w);
            add_edge(v,u,w);
        }
        dfs1(1,0);
        dfs2(1,0);
        for(int i=1;i<=n;i++) {
            max[i][0]=min[i][0]=std::max(w[i][0],w[i][2]);
        }
        for(int j=1;1<<j<=n;j++) {
            for(int i=1;i<=n-(1<<j)+1;i++) {
                max[i][j]=std::max(max[i][j-1],max[i+(1<<(j-1))][j-1]);
                min[i][j]=std::min(min[i][j-1],min[i+(1<<(j-1))][j-1]);
            }
        }
        while(m--) {
            int q=getint(),ans=0;
            for(int l=1,r=1;l<=r;l++) {
                for(;r<=n&&getsum(l,r)<=q;r++);
                ans=std::max(ans,r-l);
            }
            printf("%d
",ans);
        }
        for(int i=1;i<=n;i++) e[i].clear();
    }
}