HDU 3966 Aragorn's Story (树链剖分)

题目地址：HDU 3966
树链剖分模板题。不多说、、
代码以下：

#include <iostream> #include <string.h> #include <math.h> #include <queue> #include <algorithm> #include <stdlib.h> #include <map> #include <set> #include <stdio.h> using namespace std; #define LL long long #define pi acos(⑴.0) #pragma comment(linker,"/STACK:1024000000,1024000000") const int mod=1e9+7; const int INF=0x3f3f3f3f; const double eqs=1e⑶; const int MAXN=50000+10; #define root 1,n,1 #define lson l,mid,rt<<1 #define rson mid+1,r,rt<<1|1 int head[MAXN],cnt,a[MAXN],n; int dep[MAXN],siz[MAXN],son[MAXN],w[MAXN],top[MAXN],tot,fa[MAXN]; int sum[MAXN<<2],lazy[MAXN<<2]; struct node { int u,v,next; }edge[MAXN<<1]; void add(int u,int v) { edge[cnt].v=v; edge[cnt].u=u; edge[cnt].next=head[u]; head[u]=cnt++; } void init() { memset(head,-1,sizeof(head)); cnt=0; memset(dep,0,sizeof(dep)); memset(son,sizeof(son)); tot=0; memset(sum,sizeof(sum)); memset(lazy,sizeof(lazy)); } void dfs1(int u,int p) { siz[u]=1; for(int i=head[u];i!=-1;i=edge[i].next){ int v=edge[i].v; if(v==p) continue ; fa[v]=u; dep[v]=dep[u]+1; dfs1(v,u); siz[u]+=siz[v]; if(siz[son[u]]<siz[v]) son[u]=v; } } void dfs2(int u,int tp) { top[u]=tp;w[u]=++tot; if(son[u]) dfs2(son[u],tp); for(int i=head[u];i!=-1;i=edge[i].next){ int v=edge[i].v; if(v!=fa[u]&&v!=son[u]) dfs2(v,v); } } struct Line_Tree { void PushUp(int rt) { sum[rt]=sum[rt<<1]+sum[rt<<1|1]; } void PushDown(int l,int r,int rt) { int mid=l+r>>1; if(lazy[rt]){ sum[rt<<1]+=(mid-l+1)*lazy[rt]; sum[rt<<1|1]+=(r-mid)*lazy[rt]; lazy[rt<<1]+=lazy[rt]; lazy[rt<<1|1]+=lazy[rt]; lazy[rt]=0; } } void Update(int ll,int rr,int x,int l,int rt) { if(ll<=l&&rr>=r){ lazy[rt]+=x; sum[rt]+=x*(r-l+1); return ; } int mid=l+r>>1; PushDown(l,rt); if(ll<=mid) Update(ll,rr,x,lson); if(rr>mid) Update(ll,rson); PushUp(rt); } int Query(int p,int rt) { if(l==r){ return sum[rt]; } int mid=l+r>>1; PushDown(l,rt); if(p<=mid) return Query(p,lson); return Query(p,rson); } }lt; void Update(int u,int v,int x) { int f1=top[u],f2=top[v]; while(f1!=f2){ if(dep[f1]<dep[f2]){ swap(u,v); swap(f1,f2); } lt.Update(w[f1],w[u],root); u=fa[f1];f1=top[u]; } if(dep[u]<dep[v]){ swap(u,v); } lt.Update(w[v],root); } int main() { int m,q,i,j,u,x; char s[10]; while(scanf("%d%d%d",&n,&m,&q)!=EOF){ for(i=1;i<=n;i++){ scanf("%d",&a[i]); } init(); while(m--){ scanf("%d%d",&u,&v); add(u,v); add(v,u); } dfs1(1,-1); dfs2(1,1); for(i=1;i<=n;i++){ lt.Update(w[i],w[i],a[i],root); } while(q--){ scanf("%s",s); if(s[0]=='I'){ scanf("%d%d%d",&v,&x); Update(u,x); } else if(s[0]=='D'){ scanf("%d%d%d",-x); } else{ scanf("%d",&x); printf("%d ",lt.Query(w[x],root)); } } } return 0; }