【BZOJ】[HAOI2015]树上操作-DFS序

传送门:点击打开链接

 

题意:

有一棵点数为 N 的树,以点 1 为根,且树点有边权。有 M 个操作,分为三种:

操作 1 :把某个节点 x 的点权增加 a 。
操作 2 :把某个节点 x 为根的子树中所有点的点权都增加 a 。

操作 3 :询问某个节点 x 到根的路径中所有点的点权和。

数据范围:

 N,M<=100000 ,且所有输入数据的绝对值都不会超过 10^6 。

DFS序的板子题,网上博客很多,都比较好理解。

代码实现的时候用线段树来二分排dfs序。要注意的是,dfs序in和out要同时更新,注意正负。

用的是朴素做法,据说可以树链剖分,可我不会啊。才学就打一个朴素版吧。

//3.11:我来补代码了:树链剖分:点击打开链接

代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e5+2;
int n,m,root=1,cnt,tot,ini[N],num;
int in[N],out[N],df[600010];
int to[N<<1],next[N<<1],head[N<<1];

struct node{
	int l,r;
	int fu,zh;//zh->zheng_fu->fu
	ll sum,add;
}t[2400010];

inline void link(int u,int v){
	++tot;
	next[tot]=head[u];
	head[u]=tot;
	to[tot]=v;
}

inline void pushdown(int x){
  if(t[x].add==0) return;
    t[x<<1].sum+=1ll*(t[x<<1].zh-t[x<<1].fu)*t[x].add;
    t[x<<1|1].sum+=1ll*(t[x<<1|1].zh-t[x<<1|1].fu)*t[x].add;
   	t[x<<1].add+=t[x].add;
   	t[x<<1|1].add+=t[x].add;
   	t[x].add=0;
}

inline void maintain(int x){
	t[x].sum=t[x<<1].sum+t[x<<1|1].sum;
}

inline void dfs(int fa,int k)
{
	df[++cnt]=k;in[k]=cnt;
    for(int i=head[k];i;i=next[i]){
		if(to[i]!=fa)
		  dfs(k,to[i]);
    }
    df[++cnt]=-k;out[k]=cnt;
}

inline void build(int k,int l,int r){
	t[k].l=l,t[k].r=r;t[k].sum=t[k].add=0;
	if(l==r){
		if(df[l]>0)
		  t[k].sum=(ll)ini[df[l]],t[k].zh++;
		else t[k].sum=-(ll)ini[-df[l]],t[k].fu++;
	}else{
		int mid=(l+r)>>1;
		build(k<<1,l,mid);
		build(k<<1|1,mid+1,r);
		maintain(k);
		t[k].zh=t[k<<1].zh+t[k<<1|1].zh;
		t[k].fu=t[k<<1].fu+t[k<<1|1].fu;
	}
}

inline void ad(int k,int L,int R,ll addv)
{
	if(t[k].l>=L && t[k].r<=R){
		t[k].sum+=1ll*(t[k].zh-t[k].fu)*addv;
	    t[k].add+=addv;
	}else{
		pushdown(k);
		int mid=(t[k].l+t[k].r)>>1;
	    if(R<=mid) ad(k<<1,L,R,addv);
	    else if(L>mid) ad(k<<1|1,L,R,addv);
	    else {
		  ad(k<<1,L,mid,addv);
		  ad(k<<1|1,mid+1,R,addv);
		}
	    maintain(k);
	}
}

inline ll query(int k,int L,int R)
{
	if(t[k].l>=L && t[k].r<=R)
		return t[k].sum;
	int mid=(t[k].l+t[k].r)>>1;pushdown(k);
	if(R<=mid) return query(k<<1,L,R);
	else if(L>mid)
		return query(k<<1|1,L,R);
	else return 
		 query(k<<1,L,mid) + query(k<<1|1,mid+1,R);
}

int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
	  scanf("%d",&ini[i]);
	for(int x,y,i=1;i<n;i++){
		scanf("%d%d",&x,&y);
		link(x,y);
		link(y,x);
	}
	dfs(0,root);
	build(root,1,cnt);
	while(m--){
		int op,x,addx;
		scanf("%d",&op);
		if(op==1){
			scanf("%d%d",&x,&addx);
			ad(root,in[x],in[x],(ll)addx); 
		    ad(root,out[x],out[x],(ll)addx);
		}else if(op==2){
			scanf("%d%d",&x,&addx);
			int j=out[x];
		    ad(root,in[x],out[x],(ll)addx);
		}else{
			scanf("%d",&x);
			printf("%lld\n",query(root,1,in[x]));
		}
	}
	return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值