BZOJ4034T2
4034: [HAOI2015]T2
Time Limit: 10 Sec Memory Limit: 256 MB
Submit: 1649 Solved: 538
Description
有一棵点数为 N 的树，以点 1 为根，且树点有边权。然后有 M 个
操作，分为三种：
操作 1 ：把某个节点 x 的点权增加 a 。
操作 2 ：把某个节点 x 为根的子树中所有点的点权都增加 a 。
操作 3 ：询问某个节点 x 到根的路径中所有点的点权和。
Input
第一行包含两个整数 N, M 。表示点数和操作数。
接下来一行 N 个整数，表示树中节点的初始权值。
接下来 N-1 行每行三个正整数 fr, to ， 表示该树中存在一条边 (fr, to) 。
再接下来 M 行，每行分别表示一次操作。其中第一个数表示该操
作的种类（ 1-3 ） ，之后接这个操作的参数（ x 或者 x a ） 。
Output
对于每个询问操作，输出该询问的答案。答案之间用换行隔开。
Sample Input
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
Sample Output
6
9
13
HINT
对于 100% 的数据， N,M<=100000 ，且所有输入数据的绝对值都不会超过 10^6 。
Source
鸣谢bhiaibogf提供
代码地址：
http://paste.ubuntu.com/15280270/
一开始样例没过，区间修改没有update。。
后来数次RE，要了数据各种测，各种静态查错无果。。
发现是求和的时候又没有update。。pushdown写早了，导致叶子节点也pushdown，于是下标越界。。
突然发现学到了好多。。
感觉人生没爱了。。
不过还有kxj。。
感谢Shallwe&Claris神犇。。
黄学长的代码比较鬼畜。。所以没有%%%%
附上本蒟蒻的代码：
#includebreak;
father[x][i]=father[father[x][i-1]][i-1];
}
for (i=h[x];i;i=edge[i].next)
{
if (vis[edge[i].to])
continue;
deep[edge[i].to]=deep[x]+1;
father[edge[i].to][0]=x;
dfs1(edge[i].to);
size[x]+=size[edge[i].to];
mx[x]=max(mx[x],mx[edge[i].to]);
}
}
void dfs2(int x,int chain)
{
int k=0,i;
sz++;
pos[x]=mx[x]=sz;
belong[x]=chain;
for (i=h[x];i;i=edge[i].next)
if (deep[edge[i].to]>deep[x] && size[edge[i].to]>size[k])
k=edge[i].to;
if (k==0)
return;
dfs2(k,chain);
mx[x]=max(mx[x],mx[k]);
for (i=h[x];i;i=edge[i].next)
if (deep[edge[i].to]>deep[x] && k!=edge[i].to)
{
dfs2(edge[i].to,edge[i].to);
mx[x]=max(mx[x],mx[edge[i].to]);
}
}
void update(int s)
{
node[s].value=node[s*2].value+node[s*2+1].value;
}
void build(int s,int l,int r)
{
int mid;
node[s].l=l;
node[s].r=r;
mid=(l+r)/2;
if (l==r)
return;
build(s*2,l,mid);
build(s*2+1,mid+1,r);
update(s);
}
void paint(int s,long long z,int l,int r)
{
node[s].value+=z*(r-l+1);
delta[s]+=z;
}
void pushdown(int s,int l,int r)
{
int mid;
mid=(l+r)/2;
paint(s*2,delta[s],l,mid);
paint(s*2+1,delta[s],mid+1,r);
delta[s]=0;
}
void change(int s,int x,long long y)
{
int l=node[s].l,r=node[s].r,mid=(l+r)/2;
if (l==r)
{
node[s].value+=y;
return;
}
pushdown(s,l,r);
if (x<=mid)
change(s*2,x,y);
else
change(s*2+1,x,y);
update(s);
}
void insert(int s,int l,int r,int x,int y,long long z)
{
int mid;
if (x<=l && y>=r)
{
paint(s,z,l,r);
return;
}
pushdown(s,l,r);
mid=(l+r)/2;
if (x<=mid)
insert(s*2,l,mid,x,y,z);
if (y>mid)
insert(s*2+1,mid+1,r,x,y,z);
update(s);
}
long long querysum(int s,int x,int y)
{
int l=node[s].l,r=node[s].r,mid=(l+r)/2;
if (l==x && y==r)
return node[s].value;
if (delta[s])
pushdown(s,l,r);
long long t;
if (y<=mid)
t=querysum(s*2,x,y);
else
if (x>mid)
t=querysum(s*2+1,x,y);
else
t=querysum(s*2,x,mid)+querysum(s*2+1,mid+1,y);
update(s);
return t;
}
long long solvesum(int x,int f)
{
long long ans=0;
while (belong[x]!=belong[f])
{
ans+=querysum(1,pos[belong[x]],pos[x]);
x=father[belong[x]][0];
}
ans+=querysum(1,pos[f],pos[x]);
return ans;
}
int main()
{
int i,u,v,p,x,y;
n=read();
m=read();
for (i=1;i<=n;i++)
a[i]=read();
for (i=1;i<=n-1;i++)
{
u=read();
v=read();
add(u,v);
}
dfs1(1);
dfs2(1,1);
build(1,1,n);
for (i=1;i<=n;i++)
change(1,pos[i],a[i]);
for (i=1;i<=m;i++)
{
p=read();
if (p==1)
{
x=read();
y=read();
change(1,pos[x],y);
}
if (p==2)
{
x=read();
y=read();
//printf("%d %d\n",pos[x],mx[x]);
insert(1,1,n,pos[x],mx[x],y);
}
if (p==3)
{
x=read();
printf("%lld\n",solvesum(x,1));
}
}
return 0;
}