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序的作用就是将树分成一整块一整块的区间,这样就可以线段树维护