洛谷P3178 [HAOI2015]樹上操作
阿新 • • 發佈:2017-12-23
增加 long ons 絕對值 fin pen () ont getch
題目描述
有一棵點數為 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 。
樹鏈剖分的裸題
每次暴力更改就好
註意這題需要開long long
#include<iostream> #include<cstdio> #include<cstring> #define ls k<<1 #define rs k<<1|1 #define LL long long using namespace std; const LL MAXN=1e6+10; inline char nc() { static char buf[MAXN],*p1=buf,*p2=buf; return p1==p2&&(p1=(p2=buf)+fread(buf,1,MAXN,stdin),p1==p2)?EOF:++*p1; } inline LL read() {char c=getchar();LL x=0,f=1; while(c<‘0‘||c>‘9‘){if(c==‘-‘)f=-1;c=getchar();} while(c>=‘0‘&&c<=‘9‘){x=x*10+c-‘0‘,c=getchar();} return x*f; } LL root=1; struct node { LL u,v,w,nxt; }edge[MAXN]; LL head[MAXN]; LL num=1; inline void AddEdge(LL x,LL y) { edge[num].u=x; edge[num].v=y; edge[num].nxt=head[x]; head[x]=num++; } struct Tree { LL l,r,f,w,siz; }T[MAXN]; LL a[MAXN],b[MAXN],tot[MAXN],idx[MAXN],deep[MAXN],son[MAXN],top[MAXN],fa[MAXN],cnt=0; void update(LL k) { T[k].w=T[ls].w+T[rs].w; } void PushDown(LL k) { if(!T[k].f) return ; T[ls].w+=T[k].f*T[ls].siz; T[rs].w+=T[k].f*T[rs].siz; T[ls].f+=T[k].f; T[rs].f+=T[k].f; T[k].f=0; } LL dfs1(LL now,LL f,LL dep) { deep[now]=dep; tot[now]=1; fa[now]=f; LL maxson=-1; for(LL i=head[now];i!=-1;i=edge[i].nxt) { if(edge[i].v==f) continue; tot[now]+=dfs1(edge[i].v,now,dep+1); if(tot[edge[i].v]>maxson) maxson=tot[edge[i].v],son[now]=edge[i].v; } return tot[now]; } void dfs2(LL now,LL topf) { idx[now]=++cnt; a[cnt]=b[now]; top[now]=topf; if(!son[now]) return ; dfs2(son[now],topf); for(LL i=head[now];i!=-1;i=edge[i].nxt) if(!idx[edge[i].v]) dfs2(edge[i].v,edge[i].v); } void Build(LL k,LL ll,LL rr) { T[k].l=ll;T[k].r=rr;T[k].siz=rr-ll+1; if(ll==rr) { T[k].w=a[ll]; return ; } LL mid=(ll+rr)>>1; Build(ls,ll,mid); Build(rs,mid+1,rr); update(k); } void PointAdd(LL k,LL pos,LL val) { if(T[k].l==T[k].r) { T[k].w+=val; return ; } PushDown(k); LL mid=(T[k].l+T[k].r)>>1; if(pos<=mid) PointAdd(ls,pos,val); if(pos>mid) PointAdd(rs,pos,val); update(k); } void IntervalAdd(LL k,LL ll,LL rr,LL val) { if(ll<=T[k].l&&T[k].r<=rr) { T[k].w+=T[k].siz*val; T[k].f+=val; return ; } PushDown(k); LL mid=(T[k].l+T[k].r)>>1; if(ll<=mid) IntervalAdd(ls,ll,rr,val); if(rr>mid) IntervalAdd(rs,ll,rr,val); update(k); } LL IntervalAsk(LL k,LL ll,LL rr) { LL ans=0; if(ll<=T[k].l&&T[k].r<=rr) { ans+=T[k].w; return ans; } PushDown(k); LL mid=(T[k].l+T[k].r)>>1; if(ll<=mid) ans+=IntervalAsk(ls,ll,rr); if(rr>mid) ans+=IntervalAsk(rs,ll,rr); return ans; } LL TreeSum(LL x,LL y) { LL ans=0; while(top[x]!=top[y])//不在同一條鏈內 { if(deep[top[x]]<deep[top[y]]) swap(x,y); ans+=IntervalAsk(1,idx[top[x]],idx[x]); x=fa[top[x]]; } if(deep[x]>deep[y]) swap(x,y); ans+=IntervalAsk(1,idx[x],idx[y]); return ans; } int main() { #ifdef WIN32 freopen("a.in","r",stdin); #else #endif memset(head,-1,sizeof(head)); LL N=read(),M=read(); for(LL i=1;i<=N;i++) b[i]=read(); for(LL i=1;i<=N-1;i++) { LL x=read(),y=read(); AddEdge(x,y);AddEdge(y,x); } dfs1(root,0,1); dfs2(root,root); Build(1,1,N); while(M--) { LL opt=read(),x,val; if(opt==1) { x=read(),val=read(); PointAdd(1,idx[x],val); } else if(opt==2) { x=read(),val=read(); IntervalAdd(1,idx[x],idx[x]+tot[x]-1,val); } else { x=read(); printf("%lld\n",TreeSum(root,x)); } } return 0; }
洛谷P3178 [HAOI2015]樹上操作