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神犇。。
黄学长的代码比较鬼畜。。所以没有%%%%
附上本蒟蒻的代码：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int a[100005],n,m,cnt,sz,mx[100005],h[100005],father[100005][17],v[100005],deep[100005],size[100005],pos[100005],belong[100005];
bool vis[100005];
long long delta[400005]={0};
struct data
{
    int to,next;
};
data edge[200005];
struct lol
{
    int l,r;
    long long value;
};
lol node[400005];

int read()
{
    int w=0,c=1;
    char ch=getchar();
    while (ch<'0' || ch>'9')
      {
        if (ch=='-')
          c=-1;
        ch=getchar();
      }
    while (ch>='0' && ch<='9')
      {
        w=w*10+ch-'0';
        ch=getchar();
      }
    return w*c;
}

void add(int u,int v)
{
    cnt++;
    edge[cnt].next=h[u];
    h[u]=cnt;
    edge[cnt].to=v;
    cnt++;
    edge[cnt].next=h[v];
    h[v]=cnt;
    edge[cnt].to=u;
}

void dfs1(int x)
{
    int i;
    size[x]=1;
    vis[x]=true;
    for (i=1;i<=16;i++)
      {
        if (deep[x]<(1<<i))
          break;
        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;
}