dfs维护树上信息

dfs序可以维护树上很多信息

包涵以下几种(较为基本):

1.单点修改,子树查询

2.单点修改,树链查询

3.树链修改,单点查询

4.树链修改,子树和查询

5.子树和修改,树链查询

树链是指树上任意两点构成的路径


1.单点修改,子树查询

//在dfs序上 单点修改,区间查询
//即dfs为下标,这样保证一颗子树对应一个区间,这里用树状数组维护 
#include<bits/stdc++.h>
#define N 100005
#define M 200005
using namespace std;
int first[N],next[M],to[M],tot;
int a[N],st[N],ed[N],c[N<<2],sign;
int n,m,dep[N];
int read(){
	int cnt=0,f=1;char ch=0;
	while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
	while(isdigit(ch))cnt=cnt*10+(ch-'0'),ch=getchar();
	return cnt*f;
}
void add(int x,int y){
	next[++tot]=first[x],first[x]=tot,to[tot]=y;
}
void dfs(int u,int f){
	st[u]=++sign;//子树起点
	for(int i=first[u];i;i=next[i]){
		int t=to[i];if(t==f) continue;
		dep[t]=dep[u]+1;
		dfs(t,u);
	}ed[u]=sign;//子树终点
}
void update(int x,int val){for(;x<=n;x+=x&-x) c[x]+=val;}
int quary(int x){int ans=0;for(;x;x-=x&-x) ans+=c[x];return ans;}
int main(){
	n=read(),m=read();;
	for(int i=1;i<=n-1;i++){
		int x=read(),y=read();
		add(x,y),add(y,x);
	}
	dep[1]=1,dfs(1,0);
	while(m--){
		int op=read(),x=read();
		if(op==1){int val=read();update(st[x],val);}
		else cout<<quary(ed[x])-quary(st[x]-1)<<endl;
	}return 0;
}

 2.单点修改,树链查询

/*
树链和(x,y) 
为x到根的节点权值和+y到根的节点权值和-lca(x,y)到根的权值和-fa(lca(x,y))到根的权值和
于是变为单点查询 

对于修改a 当b为a的子孙时,a对b有贡献
于是每修改一个a 我们把a的子树都加上一个val

-> 区间修改子树,单点查询 
*/ 
#include<bits/stdc++.h>
#include<cstdio>
#define N 100005
#define M 200005
using namespace std;
int first[N],next[M],to[M],tot;
int a[N],st[N],ed[N],sign;
int n,m,dep[N],fa[N][20];
struct Node{
	int l,r,val,tag;
}t[N<<2];
int read(){
	int cnt=0,f=1;char ch=0;
	while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
	while(isdigit(ch))cnt=cnt*10+(ch-'0'),ch=getchar();
	return cnt*f;
}
void add(int x,int y){
	next[++tot]=first[x],first[x]=tot,to[tot]=y;
}
void dfs(int u,int f){
	for(int i=1;i<=18;i++)
		fa[u][i]=fa[fa[u][i-1]][i-1]; 
	st[u]=++sign;
	for(int i=first[u];i;i=next[i]){
		int t=to[i];if(t==f) continue;
		dep[t]=dep[u]+1,fa[t][0]=u;
		dfs(t,u);
	}ed[u]=sign;
}
int lca(int x,int y){
	if(dep[x]<dep[y])swap(x,y);
	for(int i=18;i>=0;i--)
		if(dep[fa[x][i]]>=dep[y]) x=fa[x][i];
	if(x==y) return x;
	for(int i=18;i>=0;i--)
		if(fa[x][i]!=fa[y][i]) 
			x=fa[x][i],y=fa[y][i];
	return fa[x][0];
}
void build(int o,int l,int r){
	t[o].l=l,t[o].r=r;
	if(l==r) return;
	int mid=(l+r)>>1;
	build(o<<1,l,mid),build(o<<1|1,mid+1,r);
} 
void spread(int o){
	if(t[o].tag){
		t[o<<1].val+=t[o].tag*(t[o<<1].r-t[o<<1].l+1);
		t[o<<1|1].val+=t[o].tag*(t[o<<1|1].r-t[o<<1|1].l+1);
		t[o<<1].tag+=t[o].tag;
		t[o<<1|1].tag+=t[o].tag;
		t[o].tag=0;
	}
}
void update(int o,int L,int R,int val){
	if(L<=t[o].l&&t[o].r<=R){
		t[o].val+=val*(t[o].r-t[o].l+1);
		t[o].tag+=val;return;
	}
	spread(o);
	int mid=(t[o].l+t[o].r)>>1;
	if(L<=mid) update(o<<1,L,R,val);
	if(R>mid) update(o<<1|1,L,R,val);
	t[o].val=t[o<<1].val+t[o<<1|1].val;
}
int quary(int o,int pos){
	if(t[o].l==t[o].r) return t[o].val;
	spread(o);
	int mid=(t[o].l+t[o].r)>>1;
	if(pos<=mid) return quary(o<<1,pos);
	else return quary(o<<1|1,pos);
}
int main(){
	n=read(),m=read();
	build(1,1,n);
	for(int i=1;i<=n-1;i++){
		int x=read(),y=read();
		add(x,y),add(y,x);
	}
	dep[1]=1,dfs(1,0);
	while(m--){
		int op=read(),x=read(),y=read();
		if(op==1) update(1,st[x],ed[x],y);
		else{
			int x1=quary(1,st[x]);//x到根的点权和 
			int y1=quary(1,st[y]);//y到根的点权和
			int l=lca(x,y); 
			int l1=quary(1,st[l]);//减去lca到根的点权和
			int l2=quary(1,st[fa[l][0]]);//减去lca's father的
			int ans=x1+y1-l1-l2; 
			cout<<ans<<endl;
		}
	}
}

3.树链修改,单点查询

/*
对于修改a 当b为a的父亲时,a对b有贡献
于是对于查询单点
可以转化为子树节点的贡献和
-> 单点修改,区间查询 
*/ 
#include<bits/stdc++.h>
#include<cstdio>
#define N 100005
#define M 200005
using namespace std;
int first[N],next[M],to[M],tot;
int a[N],st[N],ed[N],sign;
int n,m,dep[N],fa[N][20];
struct Node{
	int l,r,val;
}t[N<<2];
int read(){
	int cnt=0,f=1;char ch=0;
	while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
	while(isdigit(ch))cnt=cnt*10+(ch-'0'),ch=getchar();
	return cnt*f;
}
void add(int x,int y){
	next[++tot]=first[x],first[x]=tot,to[tot]=y;
}
void dfs(int u,int f){
	for(int i=1;i<=18;i++)
		fa[u][i]=fa[fa[u][i-1]][i-1]; 
	st[u]=++sign;
	for(int i=first[u];i;i=next[i]){
		int t=to[i];if(t==f) continue;
		dep[t]=dep[u]+1,fa[t][0]=u;
		dfs(t,u);
	}ed[u]=sign;
}
int lca(int x,int y){
	if(dep[x]<dep[y])swap(x,y);
	for(int i=18;i>=0;i--)
		if(dep[fa[x][i]]>=dep[y]) x=fa[x][i];
	if(x==y) return x;
	for(int i=18;i>=0;i--)
		if(fa[x][i]!=fa[y][i]) 
			x=fa[x][i],y=fa[y][i];
	return fa[x][0];
}
void build(int o,int l,int r){
	t[o].l=l,t[o].r=r;
	if(l==r) return;
	int mid=(l+r)>>1;
	build(o<<1,l,mid),build(o<<1|1,mid+1,r);
} 
void update(int o,int pos,int val){
	if(t[o].l==t[o].r){
		t[o].val+=val;return;
	}
	int mid=(t[o].l+t[o].r)>>1;
	if(pos<=mid) update(o<<1,pos,val);
	else update(o<<1|1,pos,val);
	t[o].val=t[o<<1].val+t[o<<1|1].val;
}
int quary(int o,int L,int R){
	if(L<=t[o].l&&t[o].r<=R) return t[o].val;
	int mid=(t[o].l+t[o].r)>>1,ans=0;
	if(L<=mid) ans+=quary(o<<1,L,R);
	if(R>mid) ans+=quary(o<<1|1,L,R);
	return ans;
}
int main(){
	n=read(),m=read();
	build(1,1,n);
	for(int i=1;i<=n-1;i++){
		int x=read(),y=read();
		add(x,y),add(y,x);
	}
	dep[1]=1,dfs(1,0);
	while(m--){
		int op=read(),x=read();
		if(op==1){
			int y=read(),val=read();
			update(1,st[x],val);
			update(1,st[y],val);
			int l=lca(x,y);
			update(1,st[l],-val);
			if(l!=1) update(1,st[fa[l][0]],-val);  
		}
		else cout<<quary(1,st[x],ed[x])<<endl;
	}
}

4.树链修改,子树和查询

/*
同样分成4个 
x对y的贡献 (x在y的子树中)
为 w[x] * (dep[x] - dep[y] + 1) = w[x] * dep[x] - w[x] * (dep[y] - 1) 
所以维护两个线段树 一个维护w[x]的和 一个维护w[x]*dep[x] 

树链修改->单点修改
子树查询-> 区间查w再乘上dep[y],区间查dep[x]*w[x] 
*/
#include<bits/stdc++.h>
#include<cstdio>
#define N 100005
#define M 200005
using namespace std;
int first[N],next[M],to[M],tot;
int a[N],st[N],ed[N],sign;
int n,m,dep[N],fa[N][20];
struct Node{
	int l,r,val;
}t1[N<<2];
struct Node1{
	int l,r,val;
}t2[N<<2];
int read(){
	int cnt=0,f=1;char ch=0;
	while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
	while(isdigit(ch))cnt=cnt*10+(ch-'0'),ch=getchar();
	return cnt*f;
}
void add(int x,int y){
	next[++tot]=first[x],first[x]=tot,to[tot]=y;
}
void dfs(int u,int f){
	for(int i=1;i<=18;i++)
		fa[u][i]=fa[fa[u][i-1]][i-1]; 
	st[u]=++sign;
	for(int i=first[u];i;i=next[i]){
		int t=to[i];if(t==f) continue;
		dep[t]=dep[u]+1,fa[t][0]=u;
		dfs(t,u);
	}ed[u]=sign;
}
int lca(int x,int y){
	if(dep[x]<dep[y])swap(x,y);
	for(int i=18;i>=0;i--)
		if(dep[fa[x][i]]>=dep[y]) x=fa[x][i];
	if(x==y) return x;
	for(int i=18;i>=0;i--)
		if(fa[x][i]!=fa[y][i]) 
			x=fa[x][i],y=fa[y][i];
	return fa[x][0];
}
void build(int o,int l,int r){
	t1[o].l=l,t1[o].r=r; t2[o].l=l,t2[o].r=r;
	if(l==r) return;
	int mid=(l+r)>>1;
	build(o<<1,l,mid),build(o<<1|1,mid+1,r);
} 
void update(int o,int pos,int dep,int val){
	if(t1[o].l==t1[o].r){
		t1[o].val+=val;
		t2[o].val+=val*dep;
		return;
	}
	int mid=(t1[o].l+t1[o].r)>>1;
	if(pos<=mid) update(o<<1,pos,dep,val);
	else update(o<<1|1,pos,dep,val);
	t1[o].val=t1[o<<1].val+t1[o<<1|1].val;
	t2[o].val=t2[o<<1].val+t2[o<<1|1].val; 
}
int quary1(int o,int L,int R){
	if(L<=t1[o].l&&t1[o].r<=R) return t1[o].val;
	int mid=(t1[o].l+t1[o].r)>>1,ans=0;
	if(L<=mid) ans+=quary1(o<<1,L,R);
	if(R>mid) ans+=quary1(o<<1|1,L,R);
	return ans;
}
int quary2(int o,int L,int R){
	if(L<=t2[o].l&&t2[o].r<=R) return t2[o].val;
	int mid=(t2[o].l+t2[o].r)>>1,ans=0;
	if(L<=mid) ans+=quary2(o<<1,L,R);
	if(R>mid) ans+=quary2(o<<1|1,L,R);
	return ans;
}
int main(){
	n=read(),m=read();
	build(1,1,n);
	for(int i=1;i<=n-1;i++){
		int x=read(),y=read();
		add(x,y),add(y,x);
	}
	dep[1]=1,dfs(1,0);
	while(m--){
		int op=read(),x=read();
		if(op==1){
			int y=read(),val=read();
			update(1,st[x],dep[x],val);
			update(1,st[y],dep[y],val);
			int l=lca(x,y);
			update(1,st[l],dep[l],-val);
			if(l!=1) update(1,st[fa[l][0]],dep[fa[l][0]],-val);
		}
		else cout<<quary2(1,st[x],ed[x])-(dep[x]-1)*quary1(1,st[x],ed[x])<<endl;
	}return 0;
}

5.子树和修改,树链查询

/*
修改x的权值,当且仅当y在x的子树中时对有贡献
贡献为 w[x] * (dep[y] - dep[x] + 1) = w[x] * (dep[y] + 1) - w[x] * dep[x]
*/ 

#include<bits/stdc++.h>
#include<cstdio>
#define N 100005
#define M 200005
using namespace std;
int first[N],next[M],to[M],tot;
int a[N],st[N],ed[N],sign;
int n,m,dep[N],fa[N][20];
struct Node{
	int l,r,val,tag;
}t1[N<<2];
struct Node1{
	int l,r,val,tag;
}t2[N<<2];
int read(){
	int cnt=0,f=1;char ch=0;
	while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
	while(isdigit(ch))cnt=cnt*10+(ch-'0'),ch=getchar();
	return cnt*f;
}
void add(int x,int y){
	next[++tot]=first[x],first[x]=tot,to[tot]=y;
}
void dfs(int u,int f){
	for(int i=1;i<=18;i++)
		fa[u][i]=fa[fa[u][i-1]][i-1]; 
	st[u]=++sign;
	for(int i=first[u];i;i=next[i]){
		int t=to[i];if(t==f) continue;
		dep[t]=dep[u]+1,fa[t][0]=u;
		dfs(t,u);
	}ed[u]=sign;
}
int lca(int x,int y){
	if(dep[x]<dep[y])swap(x,y);
	for(int i=18;i>=0;i--)
		if(dep[fa[x][i]]>=dep[y]) x=fa[x][i];
	if(x==y) return x;
	for(int i=18;i>=0;i--)
		if(fa[x][i]!=fa[y][i]) 
			x=fa[x][i],y=fa[y][i];
	return fa[x][0];
}
void build(int o,int l,int r){
	t1[o].l=l,t1[o].r=r; t2[o].l=l,t2[o].r=r;
	if(l==r) return;
	int mid=(l+r)>>1;
	build(o<<1,l,mid),build(o<<1|1,mid+1,r);
} 
void spread(int o){
	if(t1[o].tag){
		t1[o<<1].val+=t1[o].tag*(t1[o<<1].r-t1[o<<1].l+1);
		t1[o<<1|1].val+=t1[o].tag*(t1[o<<1|1].r-t1[o<<1|1].l+1);
		t1[o<<1].tag+=t1[o].tag;
		t1[o<<1|1].tag+=t1[o].tag;
		t1[o].tag=0;
	}
	if(t2[o].tag){
		t2[o<<1].val+=t2[o].tag*(t2[o<<1].r-t2[o<<1].l+1);
		t2[o<<1|1].val+=t2[o].tag*(t2[o<<1|1].r-t2[o<<1|1].l+1);
		t2[o<<1].tag+=t2[o].tag;
		t2[o<<1|1].tag+=t2[o].tag;
		t2[o].tag=0;
	}
}
void update(int o,int L,int R,int dep,int val){
	if(L<=t1[o].l&&t1[o].r<=R){
		t1[o].val+=val*(t1[o].r-t1[o].l+1);
		t1[o].tag+=val;
		t2[o].val+=val*dep*(t2[o].r-t2[o].l+1);
		t2[o].tag+=val*dep;
		return;
	}
	spread(o);
	int mid=(t1[o].l+t1[o].r)>>1;
	if(L<=mid) update(o<<1,L,R,dep,val);
	if(R>mid) update(o<<1|1,L,R,dep,val);
	t1[o].val=t1[o<<1].val+t1[o<<1|1].val;
	t2[o].val=t2[o<<1].val+t2[o<<1|1].val; 
}
int quary1(int o,int L,int R){
	if(L<=t1[o].l&&t1[o].r<=R) return t1[o].val;
	spread(o);
	int mid=(t1[o].l+t1[o].r)>>1,ans=0;
	if(L<=mid) ans+=quary1(o<<1,L,R);
	if(R>mid) ans+=quary1(o<<1|1,L,R);
	return ans;
}
int quary2(int o,int L,int R){
	if(L<=t2[o].l&&t2[o].r<=R) return t2[o].val;
	spread(o);
	int mid=(t2[o].l+t2[o].r)>>1,ans=0;
	if(L<=mid) ans+=quary2(o<<1,L,R);
	if(R>mid) ans+=quary2(o<<1|1,L,R);
	return ans;
}
int main(){
	n=read(),m=read();
	build(1,1,n);
	for(int i=1;i<=n-1;i++){
		int x=read(),y=read();
		add(x,y),add(y,x);
	}
	dep[1]=1,dfs(1,0);
	while(m--){
		int op=read(),x=read(),y=read();
		if(op==1){
			update(1,st[x],ed[x],dep[x],y);
		}
		else{
			int ans=0;
			int x1=(dep[x]+1)*quary1(1,st[x],st[x])-quary2(1,st[x],st[x]);
			int y1=(dep[y]+1)*quary1(1,st[y],st[y])-quary2(1,st[y],st[y]);
			int l=lca(x,y);
			int l1=(dep[l]+1)*quary1(1,st[l],st[l])-quary2(1,st[l],st[l]);
			ans=x1+y1-l1;
			if(l!=1) 
                        ans-=(dep[fa[l][0]]+1)*quary1(1,st[fa[l][0]],st[fa[l][0]])-quary2(1,st[fa[l][0]],st[fa[l][0]]);
			cout<<ans<<endl;
		}
	}return 0;
}

总结一下

类型修改的方法答案
1.单点修改,子树查询单点修改区间查询
2.单点修改,树链查询修改整棵子树的值,因为它对整棵子树有贡献4个单点查询
3.树链修改,单点查询树链修改变成4个单点修改

子树(区间)查询,

因为子树中的树链对这个点有影响

4.树链修改,子树和查询同样拆成4个,分别插入两个线段树,一个为w,一个为dep[当前插入的点]*w, 单点修改子树(区间)查询
5.子树和修改,树链查询

子树直接区间修改,一个修改w

另一个修改dep[当前根节点]*w的值

4个单点查询

注意: 所有线段树的下标均为dfs序,因为这样一颗子树对应一整个区间

dfs序的作用就是将树分成一整块一整块的区间,这样就可以线段树维护

 

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

FSYo

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值