SDOI 2016 游戏

时间:2021-08-01 15:18:16

树链剖分

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

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

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

/*
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 ;
}