SDOI 2016 游戏

树链剖分

线段树维护区间最小值，区间最大值

更新，对于每一个区间，找到当前区间的最小值的最大值，和要更新的值比较，如果比最大值还大，则此数对于以后的询问无任何贡献，直接返回即可，若有贡献，则一直递归到叶子节点，将最值全部更新

询问，直接询问区间最小值即可

/*
 Alice 和 Bob 在玩一个游戏。
游戏在一棵有 n 个点的树上进行。最初，每个点上都只有一个数字，那个数字是 123456789123456789。
有时，Alice 会选择一条从 s 到 t 的路径，在这条路径上的每一个点上都添加一个数字。对于路径上的一个点 r，若 r 与 s 的距离是 dis，那么 Alice 在点 r 上添加的数字是 a×dis+b。
有时，Bob 会选择一条从 s 到 t 的路径。他需要先从这条路径上选择一个点，再从那个点上选择一个数字。
Bob 选择的数字越小越好，但大量的数字让 Bob 眼花缭乱。Bob 需要你帮他找出他能够选择的最小的数字。
*/
#include<bits/stdc++.h>
#define ll long long
#define inf 123456789123456789
#define maxn 100010
using namespace std;
inline int read(){
    int s=,f=;char ch=getchar();
    for(;ch<''||ch>'';ch=getchar())if(ch=='-')f=-;
    for(;ch>=''&&ch<='';ch=getchar())s=s*+ch-'';
    return s*f;
}
int n,m;
int fa[maxn][],dep[maxn],f[maxn];
struct edge{
    int to,next,w;
}G[maxn<<];
int tot,h[maxn];
void add(int x,int y,int z){
    tot++;G[tot].to=y;G[tot].next=h[x];G[tot].w=z;h[x]=tot;
}
int pos[maxn],sz,size[maxn],poz[maxn];
ll t[maxn<<],dis[maxn],mx[maxn<<];
void dfs(int x){
    for(int i=;i<=;++i)
        if(dep[x]<(<<i))break;
        else fa[x][i]=fa[fa[x][i-]][i-];
    for(int i=h[x];i;i=G[i].next){
        int v=G[i].to;
        if(dep[v])continue;
        dep[v]=dep[x]+;
        fa[v][]=x;
        dis[v]=dis[x]+G[i].w;
        dfs(v);
        size[x]+=size[v];
    }size[x]++;
}
void dfs(int x,int ff){
    pos[x]=++sz;f[x]=ff;poz[sz]=x;
    int mx=;
    for(int i=h[x];i;i=G[i].next)
        if(dep[G[i].to]>dep[x]&&size[G[i].to]>size[mx])
            mx=G[i].to;
    if(!mx)return;
    dfs(mx,ff);
    for(int i=h[x];i;i=G[i].next)
        if(dep[G[i].to]>dep[x]&&mx!=G[i].to)
            dfs(G[i].to,G[i].to);
}
int lca(int x,int y){
    if(dep[x]<dep[y])swap(x,y);
    int d=dep[x]-dep[y];
    for(int i=;i<=;++i)
        if(d&(<<i))
            x=fa[x][i];
    if(x==y)return x;
    for(int i=;i>=;--i)
        if(fa[x][i]!=fa[y][i])
            x=fa[x][i],y=fa[y][i];
    if(x==y)return x;
    return fa[x][];
}
void build(int k,int l,int r){
    if(l==r){t[k]=inf;mx[k]=inf;return;}
    int mid=(l+r)>>;
    build(k<<,l,mid);build(k<<|,mid+,r);
    t[k]=min(t[k<<],t[k<<|]);
    mx[k]=max(mx[k<<],mx[k<<|]);
}
void init(){
    n=read();m=read();
    for(int i=;i<n;++i){
        int x=read(),y=read(),z=read();
        add(x,y,z);add(y,x,z);
    }
    dep[]=;dfs();
    dfs(,);
    build(,,n);
}
void work1(int k,int l,int r,int x,int y,ll a,ll b,ll d1,ll d2){
    if(mx[k]<=a*d1+b&&mx[k]<=a*d2+b)return;
    if(l==r){t[k]=a*d1+b,mx[k]=a*d1+b;return;}
    int mid=(l+r)>>;
    if(y<=mid)work1(k<<,l,mid,x,y,a,b,d1,d2);
    else if(x>mid)work1(k<<|,mid+,r,x,y,a,b,d1,d2);
    else work1(k<<,l,mid,x,mid,a,b,d1+dis[poz[y]]-dis[poz[mid]],d2),
         work1(k<<|,mid+,r,mid+,y,a,b,d1,d1+dis[poz[y]]-dis[poz[mid+]]);
    t[k]=min(t[k<<],t[k<<|]);
    mx[k]=max(mx[k<<],mx[k<<|]);
}
void work2(int k,int l,int r,int x,int y,ll a,ll b,ll d1,ll d2){
    if(mx[k]<=a*d1+b&&mx[k]<=a*d2+b)return;
    if(l==r){t[k]=a*d1+b,mx[k]=a*d1+b;return;}
    int mid=(l+r)>>;
    if(y<=mid)work2(k<<,l,mid,x,y,a,b,d1,d2);
    else if(x>mid)work2(k<<|,mid+,r,x,y,a,b,d1,d2);
    else work2(k<<,l,mid,x,mid,a,b,d1,d1+dis[poz[mid]]-dis[poz[x]]),
         work2(k<<|,mid+,r,mid+,y,a,b,d1+dis[poz[mid+]]-dis[poz[x]],d2);
    t[k]=min(t[k<<],t[k<<|]);
    mx[k]=max(mx[k<<],mx[k<<|]);
}
void work(int x,int y,ll a,ll b){
    int k=lca(x,y);
    int u=x;
    while(f[x]!=f[k]){
        work1(,,n,pos[f[x]],pos[x],a,b,dis[u]-dis[x],dis[u]-dis[f[x]]);x=fa[f[x]][];
    }work1(,,n,pos[k],pos[x],a,b,dis[u]-dis[x],dis[u]-dis[k]);
    while(f[y]!=f[k]){
        work2(,,n,pos[f[y]],pos[y],a,b,dis[u]+dis[f[y]]-*dis[k],dis[u]+dis[y]-*dis[k]);y=fa[f[y]][];
    }work2(,,n,pos[k],pos[y],a,b,dis[u]-dis[k],dis[u]+dis[y]-*dis[k]);
}
ll ask(int k,int l,int r,int x,int y){
    if(x<=l&&r<=y)return t[k];
    int mid=(l+r)>>;
    if(y<=mid)return ask(k<<,l,mid,x,y);
    if(x>mid)return ask(k<<|,mid+,r,x,y);
    return min(ask(k<<,l,mid,x,y),ask(k<<|,mid+,r,x,y));
}
void work(int x,int y){
    int k=lca(x,y);
    ll ans=inf;
    while(f[x]!=f[k]){
        ans=min(ans,ask(,,n,pos[f[x]],pos[x]));x=fa[f[x]][];
    }ans=min(ans,ask(,,n,pos[k],pos[x]));
    while(f[y]!=f[k]){
        ans=min(ans,ask(,,n,pos[f[y]],pos[y]));y=fa[f[y]][];
    }ans=min(ans,ask(,,n,pos[k],pos[y]));
    printf("%lld\n",ans);
}
void work(){
    for(int i=;i<=m;++i){
        int opt=read();
        if(opt==){
            int s=read(),t=read(),a=read(),b=read();
            work(s,t,a,b);
        }else{
            int s=read(),t=read();
            work(s,t);
        }
    }
}
int main(){
    init();
    work();
    return ;
}