HZOJ Dash Speed

测试点1～2：暴力。

测试点3～4：可以将边按r从大到小排序不断加入，然后用并茶几维护深度。好像也可以用猫树做。

好吧其他的部分分并没有看懂。

正解：

线段树分治，求出每个速度的答案。

对于速度区间$[L,R]$，将完全包含这个区间的边加入，对于其余的边，按照和mid的关系分到左右儿子，这里的一条边有可能同时分到两个儿子所以直接塞vector就行了。

那么到达叶子区间时，满足条件的树的结构已经出来了，答案就是这些联通块的直径。

那么大体思路已经清晰了，如何动态维护树的结构呢？lct！！！并查集即可，维护联通块直径的两个端点，合并时分6中情况讨论即可，注意由于递归完左儿子后还要处理右儿子，所以左儿子中连的边都要清掉，于是并查集不能路径压缩，按秩合并即可（然而我故意把size大的合并到size小的还是A掉了）。

#include<iostream>
#include<cstring>
#include<cstdio>
#include<vector>
#define re register
#define co const
#define rec re co
#define LL long long
using namespace std;
struct edge
{
    int u,v,nxt,l,r;
    #define u(x) ed[x].u
    #define v(x) ed[x].v
    #define l(x) ed[x].l
    #define r(x) ed[x].r
    #define n(x) ed[x].nxt
    friend bool operator < (edge a,edge b)
    {return a.r>b.r;}
}ed[500010];
int first[500010],num_e=1;
#define f(x) first[x]
int n,m;
int f[200010][21],dep[200010],size[200010];
int ttop[200010],son[200010],dfn[200010],id[200010],cnt;
void dfs1(int x)
{
    size[x]=1;
    for(int i=f(x);i;i=n(i))
    if(v(i)!=f[x][0])
    {
        f[v(i)][0]=x;dep[v(i)]=dep[x]+1;
        dfs1(v(i));size[x]+=size[v(i)];
        if(size[v(i)]>size[son[x]])son[x]=v(i);
    }
}
void dfs2(int x,int t)
{
    ttop[x]=t;dfn[x]=++cnt;id[cnt]=x;
    if(son[x])dfs2(son[x],t);
    for(int i=f(x);i;i=n(i))
    if(v(i)!=f[x][0]&&v(i)!=son[x])
    dfs2(v(i),v(i));
}
inline int LCA(int x,int y)
{
    while(ttop[x]!=ttop[y])
    {
        if(dep[ttop[x]]>dep[ttop[y]])swap(x,y);
        y=f[ttop[y]][0];
    }
    if(dep[x]>dep[y])swap(x,y);
    return x;
}
#define st sta[top]
#define st2 sta2[top2]
#define st3 sta3[top3]
struct node{int id,pu,pv;}sta[500010];int top;
struct node2{int id,of;}sta2[500010];int top2;
struct node3{int id,ad;}sta3[500010];int top3;
vector<int> in[800010];
int fa[200010],u[200010],v[200010],al[200010],siz[200010];
int get(re int x){return fa[x]==x?x:get(fa[x]);}
LL lenth(re int u,re int v){return dep[u]+dep[v]-dep[LCA(u,v)]*2;}
void solve(re int l,re int r,rec int x,re int res)
{
    int mid=(l+r)>>1,tem=top,tem2=top2,tres=res,tem3=top3;
    for(int i=0;i<in[x].size();i++)
    if(l(in[x][i])<=l&&r(in[x][i])>=r)
    {
        
        int f1=get(u(in[x][i])),f2=get(v(in[x][i]));if(f1==f2)continue;
        int gl;
        if(siz[f1]>siz[f2])
        {
            fa[f1]=f2;gl=f2;siz[f2]+=siz[f1];
            sta[++top]=(node){f2,u[f2],v[f2]};
            sta2[++top2]=(node2){f1,f1};
            sta3[++top3]=(node3){f2,siz[f1]};
        }
        else
        {
            fa[f2]=f1;gl=f1;siz[f1]+=siz[f2];
            sta[++top]=(node){f1,u[f1],v[f1]};
            sta2[++top2]=(node2){f2,f2};
            sta3[++top3]=(node3){f1,siz[f2]};
        }
        int t1=u[f1],t2=u[f2],t3=v[f1],t4=v[f2];
        int max1=lenth(t1,t2),
            max2=lenth(t1,t3),
            max3=lenth(t1,t4),
            max4=lenth(t2,t3),
            max5=lenth(t2,t4),
            max6=lenth(t3,t4);
        int maxn=max(max1,max2);maxn=max(maxn,max3);
        maxn=max(maxn,max4);maxn=max(maxn,max5);maxn=max(maxn,max6);
             if(maxn==max1)u[gl]=t1,v[gl]=t2;
        else if(maxn==max2)u[gl]=t1,v[gl]=t3;
         else if(maxn==max3)u[gl]=t1,v[gl]=t4;
        else if(maxn==max4)u[gl]=t2,v[gl]=t3;
        else if(maxn==max5)u[gl]=t2,v[gl]=t4;
        else if(maxn==max6)u[gl]=t3,v[gl]=t4;
        res=max(res,maxn);
    }
    else 
    {
        if(l(in[x][i])<=mid)in[x*2].push_back(in[x][i]);
        if(r(in[x][i])> mid)in[x*2+1].push_back(in[x][i]);
    }
    if(l==r)
    {
        al[l]=res;
        while(top>tem)u[st.id]=st.pu,v[st.id]=st.pv,top--;
        while(top2>tem2)fa[st2.id]=st2.of,top2--;
        while(top3>tem3)siz[st3.id]-=st3.ad,top3--;
        res=tres;
        return;
    }
    solve(l,mid,x*2,res);solve(mid+1,r,x*2+1,res);
    while(top>tem)u[st.id]=st.pu,v[st.id]=st.pv,top--;
    while(top2>tem2)fa[st2.id]=st2.of,top2--;
    while(top3>tem3)siz[st3.id]-=st3.ad,top3--;
    res=tres;
}
inline int read();
inline void add(rec int u,rec int v,rec int l,rec int r);
signed main()
{    
//    freopen("speed1.in","r",stdin);
//    freopen("out.out","w",stdout);
    
    n=read(),m=read();
    int tu,tv,tl,tr;
    for(re int i=1;i<n;i++)
    {
        tu=read(),tv=read(),tl=read(),tr=read();
        add(tu,tv,tl,tr);add(tv,tu,tl,tr);
    }
    dfs1(1),dfs2(1,1);
    for(re int i=1;i<=20;i++)    
        for(re int j=1;j<=n;j++)
            f[j][i]=f[f[j][i-1]][i-1];
    for(re int i=1;i<=n;i++)fa[i]=u[i]=v[i]=i,siz[i]=1;
    for(re int i=2;i<=num_e;i+=2)in[1].push_back(i);
    solve(1,n,1,0);
    int q;
    for(re int i=1;i<=m;i++)
    {
        q=read();
        printf("%d\n",al[q]);
    }
}
inline int read()
{
    int s=0,f=1;char a=getchar();
    while(a<'0'||a>'9'){if(a=='-')f=-1;a=getchar();}
    while(a>='0'&&a<='9'){s=s*10+a-'0';a=getchar();}
    return s*f;
}
inline void add(rec int u,rec int v,rec int l,rec int r)
{
    ++num_e;
    u(num_e)=u;
    v(num_e)=v;
    l(num_e)=l;
    r(num_e)=r;
    n(num_e)=f(u);
    f(u)=num_e;
}