牛客练习赛61 ——F.苹果树（动态点分）

题意：

   链接：https://ac.nowcoder.com/acm/contest/5026/F

   给定一个由n个点和n-1条边组成的无向连通图（即一棵树），每一条边都有一个权值代表走过它所需要的时间花费，每个点上面有初始有一个苹果，每个苹果有一个成熟度。

       由于苹果不熟或者熟的太透了都不会好吃，所以现在米咔想摘一个成熟度在范围[x,y]的苹果来做苹果派。现在需要你来执行m条操作。

       1 u x：在某一个点u新出现了一个成熟度为x的苹果。（1≤u≤n,1≤x≤10000）

       2 u x y：询问米咔从某一个点u出发，去摘一个成熟度在[x，y]范围内的苹果并回到u的最小时间花费。（1≤u≤n,1≤x≤y≤10000）

n，m（1 ≤ n，m ≤100000）

题解：

     比赛时头铁写了个 O(n*logn*logn*logn) 的树剖+线段树套线段树，结果直接T飞，当场去世...赛后看了聚聚们的代码和题解后，才想起来动态点分治这玩意，差不多快一年没写过了（菜得飞起...

     动态点分裸题。。每个点开个动态线段树，更新时沿点分树的父链跳log次，更新每层父节点的信息。查询时也是沿父链暴跳求最小值。由于是求最小值，查询时有可能会查询同一棵子树内的节点，但这个最小值在从下往上查询时就一定会被记录到，所以就不需要考虑去重问题了。

#include<bits/stdc++.h>
using namespace std;
const long long N=1e5+10,inf=1e18;
int p[N*2],nex[N*2],head[N*2],w[N*2],e;
int A[N],vis[N],siz[N],ma[N],root,h[N],fa[N][20],dis[N][20];
void add(int a,int b,int c){
   p[++e]=b;
   nex[e]=head[a];
   head[a]=e;
   w[e]=c;
}
void dfsroot(int u,int f,int sum){
    int i,v;
    siz[u]=1;
    ma[u]=0;
    for(i=head[u];i;i=nex[i]){
        v=p[i];
        if(vis[v]==0&&v!=f){
            dfsroot(v,u,sum);
            siz[u]+=siz[v];
            ma[u]=max(ma[u],siz[v]);
        }
    }
    ma[u]=max(ma[u],sum-siz[u]);
    if(root==0||ma[root]>ma[u]) root=u;
}
void dfsdis(int u,int f,int ff,int d){
     h[u]++;
     fa[u][h[u]]=ff;
     dis[u][h[u]]=d;
     int i,v;
     for(i=head[u];i;i=nex[i]){
         v=p[i];
         if(v!=f&&vis[v]==0) dfsdis(v,u,ff,d+w[i]);
     }
}
void dfs(int u){
    vis[u]=1;
    int i,v;
    dfsdis(u,u,u,0);
    for(i=head[u];i;i=nex[i]){
        v=p[i];
        if(vis[v]==0){
            root=0;
            dfsroot(v,u,siz[v]);
            dfs(root);
        }
    }
}

struct Tree{
    int L,R;
    long long mi;
}T[30000000];
int cnt,rt[N];
void update(int &x,int l,int r,int k,long long v){
    T[++cnt]=T[x];
    x=cnt;
    T[x].mi=min(T[x].mi,v);
    if(l==r) return;
    int mid=(l+r)/2;
    if(k<=mid) update(T[x].L,l,mid,k,v);
    else update(T[x].R,mid+1,r,k,v);
}
long long qmin(int x,int l,int r,int L,int R){
    if(l==L&&r==R) return T[x].mi;
    int mid=(l+r)/2;
    if(R<=mid) return qmin(T[x].L,l,mid,L,R);
    else if(L>mid) return qmin(T[x].R,mid+1,r,L,R);
    else return min(qmin(T[x].L,l,mid,L,mid),qmin(T[x].R,mid+1,r,mid+1,R));
}
void change(int u,int v){
    for(int i=h[u];i;i--) update(rt[fa[u][i]],1,10000,v,dis[u][i]);
}
long long findans(int u,int l,int r){
    long long ans=inf,i;
    for(i=h[u];i;i--) ans=min(qmin(rt[fa[u][i]],1,10000,l,r)+dis[u][i],ans);
    return ans;
}

int main(){
    int i,n,m,a,b,c,d;
    scanf("%d%d",&n,&m);
    for(i=1;i<=n;i++) scanf("%d",&A[i]);
    for(i=1;i<n;i++){
        scanf("%d%d%d",&a,&b,&c);
        add(a,b,c);
        add(b,a,c);
    }
    T[0].L=T[0].R=0;
    T[0].mi=1e18;
    dfsroot(1,0,n);
    dfs(root);
    for(i=1;i<=n;i++) change(i,A[i]);
    while(m--){
        scanf("%d",&a);
        if(a==1){
            scanf("%d%d",&b,&c);
            change(b,c);
        }
        else{
            scanf("%d%d%d",&b,&c,&d);
            long long ans=findans(b,c,d);
            if(ans>=1e16) ans=-1;
            else ans*=2;
            printf("%lld\n",ans);
        }
    }
}