luogu3250 网络 (整体二分+树上差分+树状数组)

首先整体二分，问题变成是否存在经过一个点的满足条件的路径

那么我对于每个路径(a,b,lca)，在树状数组的dfn[a]++,dfn[b]++,dfn[lca]--,dfn[fa[lca]--]

然后直接查那个点的子树和就行了

 #include<bits/stdc++.h>
 #define CLR(a,x) memset(a,x,sizeof(a))
 #define MP make_pair
 using namespace std;
 typedef long long ll;
 typedef unsigned long long ull;
 typedef pair<int,int> pa;
 const int maxn=2e5+,maxm=4e5+;
 
 inline ll rd(){
     ll x=;char c=getchar();int neg=;
     while(c<''||c>''){if(c=='-') neg=-;c=getchar();}
     while(c>=''&&c<='') x=x*+c-'',c=getchar();
     return x*neg;
 }
 
 int N,M,eg[maxn*][],egh[maxn],ect;
 int dfn[maxn][],tot,fa[maxn][],dep[maxn];
 int tr[maxn],ans[maxm],que[maxm];
 
 inline int lowbit(int x){return x&(-x);}
 inline void add(int x,int y){
     for(;x&&x<=N;x+=lowbit(x)) tr[x]+=y;
 }
 inline int query(int x){
     int re=;for(;x;x-=lowbit(x)) re+=tr[x];return re;
 }
 
 struct Node{
     int d,a,b,v,lca;
     Node(int x=,int y=,int z=,int l=,int k=){
         d=x,a=y,b=z,v=l,lca=k;
     }
     inline void cover(int x){
         add(dfn[a][],d*x);
         add(dfn[b][],d*x);
         add(dfn[lca][],-d*x);
         add(dfn[fa[lca][]][],-d*x);
     }
 }op[maxm],tmp[maxm];
 
 inline void adeg(int a,int b){
     eg[++ect][]=b,eg[ect][]=egh[a],egh[a]=ect;
 }
 
 inline void dfs(int x){
     for(int i=;fa[x][i]&&fa[fa[x][i]][i];i++)
         fa[x][i+]=fa[fa[x][i]][i];
     dfn[x][]=++tot;
     for(int i=egh[x];i;i=eg[i][]){
         int b=eg[i][];if(b==fa[x][]) continue;
         dep[b]=dep[x]+,fa[b][]=x;dfs(b);
     }dfn[x][]=tot;
 }
 
 inline int getlca(int x,int y){
     if(dep[x]<dep[y]) swap(x,y);
     for(int i=log2(dep[x]-dep[y]);i>=&&dep[x]!=dep[y];i--){
         if(dep[fa[x][i]]>=dep[y])
             x=fa[x][i];
     }
     if(x==y) return x;
     for(int i=log2(dep[x]);i>=;i--){
         if(fa[x][i]!=fa[y][i])
             x=fa[x][i],y=fa[y][i];
     }return fa[x][];
 }
 
 inline void solve(int l,int r,int vl,int vr){
     if(l>r||vl>vr) return;
     int a=l-,b=r+,mid=vl+vr>>;
     int cnt=;
 
     for(int i=l;i<=r;i++){
         if(op[i].d){
             if(op[i].v>=mid) tmp[--b]=op[i],op[i].cover(),cnt+=op[i].d;
             else tmp[++a]=op[i];
         }else{
             int re=query(dfn[op[i].a][])-query(dfn[op[i].a][]-);
             if(cnt-re>=) tmp[--b]=op[i],ans[op[i].b]=mid;
             else tmp[++a]=op[i];
         }
     }
 
     for(int i=l;i<=r;i++){
         if(op[i].d&&op[i].v>=mid) op[i].cover(-);
     }
 
     for(int i=l;i<=a;i++) op[i]=tmp[i];
     for(int i=r;i>=b;i--) op[i]=tmp[r-i+b];
     solve(l,a,vl,mid-);solve(b,r,mid+,vr);
 }
 
 int main(){
     // freopen("network4.in","r",stdin);
     int i,j,k;
     N=rd(),M=rd();
     for(i=;i<N;i++){
         int a=rd(),b=rd();
         adeg(a,b);adeg(b,a);
     }
     dep[]=;dfs();
     for(i=,j=;i<=M;i++){
         int o=rd();
         if(o==){
             int a=rd(),b=rd(),c=rd();
             op[i]=Node(,a,b,c,getlca(a,b));
         }else if(o==){
             int t=rd();
             op[i]=Node(-,op[t].a,op[t].b,op[t].v,op[t].lca);
         }else{
             op[i]=Node(,rd(),i,,);
             que[++j]=i;
         }
     }
     CLR(ans,-);
     solve(,M,,1e9);
     for(i=;i<=j;i++){
         printf("%d\n",ans[que[i]]);
     }
     return ;
 }