1. 程式人生 > 實用技巧 >樹鏈剖分模板

樹鏈剖分模板

題目連結:https://www.luogu.com.cn/problem/P3178

題目描述

有一棵點數為 N 的樹,以點 1 為根,且樹點有邊權。然後有 M 個操作,分為三種:

  • 操作 1 :把某個節點 x 的點權增加 a 。
  • 操作 2 :把某個節點 x 為根的子樹中所有點的點權都增加 a 。
  • 操作 3 :詢問某個節點 x 到根的路徑中所有點的點權和。

輸入格式

第一行包含兩個整數 N, M 。表示點數和運算元。
接下來一行 N 個整數,表示樹中節點的初始權值。
接下來 N-1 行每行兩個正整數 from, to , 表示該樹中存在一條邊 (from, to) 。
再接下來 M 行,每行分別表示一次操作。其中第一個數表示該操作的種類( 1-3 ) ,之後接這個操作的引數( x 或者 x a ) 。

輸出格式

對於每個詢問操作,輸出該詢問的答案。答案之間用換行隔開。

輸入輸出樣例

輸入 #1
5 5
1 2 3 4 5
1 2
1 4
2 3
2 5
3 3
1 2 1
3 5
2 1 2
3 3
輸出 #1
6
9
13

說明/提示

對於 100% 的資料, N,M<=100000 ,且所有輸入資料的絕對值都不會超過 10^6 。

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e5+9;
const double ep=1e-6;
const int mod=998244353;
const int MAX=1e9; #define mk make_pair #define PII pair<int,int> #define PLL pair<ll,ll> #define pb push_back typedef long long ll; ll n,m,b[maxn],size[maxn],d[maxn]; ll f[maxn][20]; vector<int>vec[maxn]; ll be[maxn]; ll pp[maxn],cnt; struct node{ ll l,r,sum,lazy; }a[maxn<<2
]; void dfs(int x,int p) { size[x]=1; for(int j:vec[x]) { if(j==p)continue; f[j][0]=x; d[j]=d[x]+1; dfs(j,x); size[x]+=size[j]; } } void dfs1(int x,int chain,int p) { be[x]=chain; pp[x]=++cnt; int pos=0; for(int j:vec[x]) { if(j==p)continue; if(size[j]>size[pos])pos=j; } if(!pos)return; dfs1(pos,chain,x); for(int j:vec[x]) { if(j==p||j==pos)continue; dfs1(j,j,x); } } int lca(int a,int b) { if(d[a]<d[b])swap(a,b); for(int i=17;i>=0;i--)if((d[a]-d[b])>>i)a=f[a][i]; if(a==b)return a; for(int i=17;i>=0;i--)if(f[a][i]!=f[b][i])a=f[a][i],b=f[b][i]; return f[a][0]; } void build(int k,int l,int r) { a[k].l=l,a[k].r=r; if(l==r)return; int mid=(l+r)>>1; build(k<<1,l,mid); build(k<<1|1,mid+1,r); } void pushup(int k) { a[k].sum=a[k<<1].sum+a[k<<1|1].sum; } void pushdown(int k) { if(a[k].l==a[k].r) { a[k].lazy=0;return; } a[k<<1].sum+=(a[k<<1].r-a[k<<1].l+1)*a[k].lazy; a[k<<1|1].sum+=(a[k<<1|1].r-a[k<<1|1].l+1)*a[k].lazy; a[k<<1].lazy+=a[k].lazy; a[k<<1|1].lazy+=a[k].lazy; a[k].lazy=0; } void change(int k,int l,int r,int x) { if(a[k].l>=l&&a[k].r<=r) { a[k].sum+=(a[k].r-a[k].l+1)*x; a[k].lazy+=x; return; } if(a[k].lazy)pushdown(k); int mid=(a[k].l+a[k].r)>>1; if(l<=mid)change(k<<1,l,r,x); if(r>mid)change(k<<1|1,l,r,x); pushup(k); } ll query_SUM(int k,int l,int r) { if(a[k].lazy)pushdown(k); if(a[k].l>=l&&a[k].r<=r) { return a[k].sum; } int mid=(a[k].l+a[k].r)>>1; if(r<=mid)return query_SUM(k<<1,l,r); else if(l>mid)return query_SUM(k<<1|1,l,r); else return query_SUM(k<<1,l,mid)+query_SUM(k<<1|1,mid+1,r); } ll query(int x) { ll sum1=0; while(be[x]!=1) { sum1+=query_SUM(1,pp[be[x]],pp[x]); x=f[be[x]][0]; } sum1+=query_SUM(1,1,pp[x]); return sum1; } int main() { scanf("%lld%lld",&n,&m); for(int i=1;i<=n;i++)scanf("%lld",&b[i]); for(int i=1;i<n;i++) { int u,v; scanf("%d%d",&u,&v); vec[u].pb(v); vec[v].pb(u); } dfs(1,0); for(int i=1;i<=17;i++)for(int j=1;j<=n;j++)f[j][i]=f[f[j][i-1]][i-1]; dfs1(1,1,-1); build(1,1,n); for(int i=1;i<=n;i++)change(1,pp[i],pp[i],b[i]); while(m--) { int op,x,y; scanf("%d%d",&op,&x); if(op==1) { scanf("%d",&y); change(1,pp[x],pp[x],y); } if(op==2) { scanf("%d",&y); change(1,pp[x],pp[x]+size[x]-1,y); } if(op==3) { printf("%lld\n",query(x)); } } }