【BZOJ1146】网络管理（整体二分）
【BZOJ1146】网络管理（整体二分）
题面
题解
要看树套树的戳这里
毕竟是：智商不够数据结构来补
所以，
我们来当一回智商够的选手
听说主席树的题目大部分都可以整体二分
这题丢进去整体二分就行了
只需要查询树上贡献的前缀和
但是，对于这种带修改的
不能够直接修改
每一个修改操作必须拆成两个：
一个删掉原来的修改，一个加入现在的修改
否则就会重复计算贡献从而\(gg\)
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define ll long long
#define RG register
#define MAX 81000
#define lson (now<<1)
#define rson (now<<1|1)
inline int read()
{
RG int x=0,t=1;RG char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=-1,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return x*t;
}
struct Line{int v,next;}e[MAX<<1];
int h[MAX],cnt=1;
inline void Add(int u,int v){e[cnt]=(Line){v,h[u]};h[u]=cnt++;}
int dfn[MAX],low[MAX],fa[MAX],size[MAX],hson[MAX],dep[MAX],top[MAX];
int tim,n,Q;
int LCA(int u,int v)
{
while(top[u]!=top[v])
dep[top[u]]<dep[top[v]]?v=fa[top[v]]:u=fa[top[u]];
return dep[u]<dep[v]?u:v;
}
void dfs1(int u,int ff)
{
size[u]=1;fa[u]=ff;dep[u]=dep[ff]+1;
for(int i=h[u];i;i=e[i].next)
{
if(e[i].v==ff)continue;
dfs1(e[i].v,u);
size[u]+=size[e[i].v];
if(size[e[i].v]>size[hson[u]])hson[u]=e[i].v;
}
}
void dfs2(int u,int tp)
{
top[u]=tp;dfn[u]=++tim;
if(hson[u])dfs2(hson[u],tp);
for(int i=h[u];i;i=e[i].next)
{
if(e[i].v==fa[u]||e[i].v==hson[u])continue;
dfs2(e[i].v,e[i].v);
}
low[u]=tim;
}
int sum,ans[MAX];
struct Op{int k,a,b,id,opt;}O[MAX<<1],p1[MAX<<1],p2[MAX<<1];
int c[MAX],TT[MAX];
int lowbit(int x){return x&(-x);}
void add(int x,int w){while(x<=n)c[x]+=w,x+=lowbit(x);}
int getsum(int x){int ret=0;while(x)ret+=c[x],x-=lowbit(x);return ret;}
int jump(int u,int v)
{
int ret=0;
while(top[u]!=top[v])
{
if(dep[top[u]]<dep[top[v]])swap(u,v);
ret+=getsum(dfn[u])-getsum(dfn[top[u]]-1);
u=fa[top[u]];
}
if(dep[u]<dep[v])swap(u,v);
ret+=getsum(dfn[u])-getsum(dfn[v]-1);
return ret;
}
void Work(int L,int R,int l,int r)
{
if(L>R)return;
if(l==r)
{
for(int i=L;i<=R;++i)
if(!O[i].opt&&ans[O[i].id]!=-1)ans[O[i].id]=l;
return;
}
int mid=(l+r)>>1;
int t1=0,t2=0;
for(int i=L;i<=R;++i)
{
if(O[i].k==0)
{
if(O[i].b<=mid)p1[++t1]=O[i];
else p2[++t2]=O[i],add(dfn[O[i].a],O[i].id);
}
else
{
int ss=jump(O[i].a,O[i].b);
if(ss>=O[i].k)p2[++t2]=O[i];
else O[i].k-=ss,p1[++t1]=O[i];
}
}
for(int i=L;i<=R;++i)if(O[i].opt&&O[i].b>mid)add(dfn[O[i].a],-O[i].id);
for(int i=1;i<=t1;++i)O[L+i-1]=p1[i];
for(int i=1;i<=t2;++i)O[L+t1+i-1]=p2[i];
Work(L,L+t1-1,l,mid);Work(L+t1,R,mid+1,r);
}
int main()
{
n=read();Q=read();
int Cnt=0;
for(int i=1;i<=n;++i)
{
TT[i]=read();
O[++Cnt]=(Op){0,i,TT[i],1,1};
}
for(int i=1;i<n;++i)
{
int u=read(),v=read();
Add(u,v);Add(v,u);
}
dfs1(1,0);dfs2(1,1);
for(int i=1;i<=Q;++i)
{
int k=read(),a=read(),b=read();
if(k)
{
O[++Cnt]=(Op){k,a,b,++sum,0};
int len=dep[O[Cnt].a]+dep[O[Cnt].b]-2*dep[LCA(O[Cnt].a,O[Cnt].b)]+1;
if(O[Cnt].k>len)ans[O[Cnt].id]=-1;
}
else
{
O[++Cnt]=(Op){0,a,TT[a],-1,1};
O[++Cnt]=(Op){0,a,TT[a]=b,1,1};
}
}
Work(1,Cnt,1,1e8);
for(int i=1;i<=sum;++i)
ans[i]!=-1?printf("%d\n",ans[i]):puts("invalid request!");
return 0;
}
