树链剖分基本步骤
总体步骤:init(); dfs1(1); dfs2(1,1); build(1,1,n); cl(); dfs1 处理出size[sj](子节点数),son[sj](重儿子),dep[sj](深度),vi[sj](点权,把边权放在子节点上)。 dfs2 处理出top[sj](所在链的链首,轻儿子从子节点重新开始,重儿子从父节点继承),id[sj](点出现的位置,用cnt标记),pos[sj](cnt位置的边连着哪一个点)。 id[x]=++cnt; pos[cnt]=x; if(son[x]) dfs2(son[x],y); for(int i=h[x];i!=-1;i=b[i].ne) if(b[i].v!=fa[x]&&b[i].v!=son[x]) dfs2(b[i].v,b[i].v); build if(z==y) mx[x]=vi[pos[y]]; 这个位置链的最大值等于这个位置的链所连点的权值。一般线段树方法建树。和线段树有关的数组要开四倍四倍四倍。。。 cl 分change和query Change 改变的边编号要按输入数据检索,从线段树上找到这条边上作为子节点的点改变它的点权。 int now=(b1<<1)-1; int fr=b[now].u,to=b[now].v; if(fa[fr]==to) change(id[fr],b2,1,1,n); else change(id[to],b2,1,1,n); Query query函数即为一般的线段树查询,把不在一条链上的两个点向一条链上靠拢,时刻注意要求哪点在下层。如果本来就是点权,则id[x]不应该+1且id[x]==id[y]的情况是可行的。 int fx=top[x],fy=top[y],res=1<<31; while(fx^fy) { if(dep[fx]<dep[fy]) { jh(x,y); jh(fx,fy); } res=qd(res,query(id[fx],id[x],1,1,n)); x=fa[fx]; fx=top[x]; } if(dep[x]>dep[y]) jh(x,y); if(id[x]<id[y]) res=qd(res,query(id[x]+1,id[y],1,1,n)); return res; 零散的小函数 qd取较大值,jh取地址交换,add邻接表加边(要有u、v、w、ne)
因为我太弱了不写就会很快忘掉QAQ
南风知我意,吹梦到西洲。