[HAOI2015]树上操作（树链剖分，线段树）

题目描述

有一棵点数为 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 ） 。

输出格式：

对于每个询问操作，输出该询问的答案。答案之间用换行隔开。

思路：

树剖板子题，不解释了……

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<bitset>
#include<queue>
#include<algorithm>
#define int long long
#define rii register int i
#define rij register int j
using namespace std;
int n,m;
struct ljb{
    int to,nxt;
}x[];
struct node{
    int val,lazy;
}y[];
int f[],sd[],head[],val[],size[];
int weison[],cnt,bnt,nid[],nval[],ntop[];
int res;
void add(int from,int to)
{
    cnt++;
    x[cnt].to=to;
    x[cnt].nxt=head[from];
    head[from]=cnt;
}
void dfs1(int wz,int fa,int s)
{
    sd[wz]=s;
    f[wz]=fa;
    size[wz]=;
    int maxn=;
    for(rii=head[wz];i!=;i=x[i].nxt)
    {
        int to=x[i].to;
        if(to==fa)
        {
            continue;
        }
        dfs1(to,wz,s+);
        size[wz]+=size[to];
        if(size[to]>maxn)
        {
            maxn=size[to];
            weison[wz]=to;
        }
    }
}
void dfs2(int wz,int top)
{
    bnt++;
    nid[wz]=bnt;
    nval[bnt]=val[wz];
    ntop[wz]=top;
    if(weison[wz]==)
    {
        return;
    }
    dfs2(weison[wz],top);
    for(rii=head[wz];i!=;i=x[i].nxt)
    {
        int to=x[i].to;
        if(weison[wz]==to||f[wz]==to)
        {
            continue;
        }
        dfs2(to,to);
    }
}
void build(int bh,int l,int r)
{
    if(l==r)
    {
        y[bh].val=nval[l];
        return;
    }
    int mid=(l+r)/;
    build(bh*,l,mid);
    build(bh*+,mid+,r);
    y[bh].val+=y[bh*].val+y[bh*+].val;
}
void pushdown(int bh,int cd)
{
    y[bh*].lazy+=y[bh].lazy;
    y[bh*+].lazy+=y[bh].lazy;
    y[bh*].val+=y[bh].lazy*(cd-(cd/));
    y[bh*+].val+=y[bh].lazy*(cd/);
    y[bh].lazy=;
}
void updata(int bh,int nl,int nr,int l,int r,int val)
{
    int len=(nr-nl+);
    if(l<=nl&&nr<=r)
    {
        y[bh].lazy+=val;
        y[bh].val+=val*len;
        return;
    }
    if(y[bh].lazy!=)
    {
        pushdown(bh,len);
    }
    int mid=(nl+nr)/;
    if(l<=mid)
    {
        updata(bh*,nl,mid,l,r,val);
    }
    if(r>mid)
    {
        updata(bh*+,mid+,nr,l,r,val);
    }
    y[bh].val=(y[bh*].val+y[bh*+].val);
}
void query(int bh,int nl,int nr,int l,int r)
{
    if(l<=nl&&r>=nr)
    {
        res+=y[bh].val;
        return;
    }
    int mid=(nl+nr)/;
    if(y[bh].lazy!=)
    {
        pushdown(bh,nr-nl+);
    }
    if(l<=mid)
    {
        query(bh*,nl,mid,l,r);
    }
    if(r>mid)
    {
        query(bh*+,mid+,nr,l,r);
    }
}
int querylj(int from,int to)
{
    int ans=;
    while(ntop[from]!=ntop[to])
    {
        if(sd[ntop[from]]<sd[ntop[to]])
        {
            swap(from,to);
        }
        res=;
        query(,,n,nid[ntop[from]],nid[from]);
        ans+=res;
        from=f[ntop[from]];
    }
    if(sd[from]>sd[to])
    {
        swap(from,to);
    }
    res=;
    query(,,n,nid[from],nid[to]);
    ans+=res;
    return ans;
}
void addson(int wz,int val)
{
    updata(,,n,nid[wz],nid[wz]+size[wz]-,val);
}
signed main()
{
    scanf("%lld%lld",&n,&m);
    for(rii=;i<=n;i++)
    {
        scanf("%lld",&val[i]);
    }
    for(rii=;i<n;i++)
    {
        int from,to;
        scanf("%lld%lld",&from,&to);
        add(from,to);
        add(to,from);
    }
    dfs1(,,);
    dfs2(,);
    build(,,n);
    for(rii=;i<=m;i++)
    {
        int opt,from,to,val,kkk;
        scanf("%lld",&opt);
        if(opt==)
        {
            scanf("%lld%lld",&from,&val);
            updata(,,n,nid[from],nid[from],val);
        }
        if(opt==)
        {
            scanf("%lld%lld",&from,&val);
            addson(from,val);
        }
        if(opt==)
        {
            scanf("%lld",&from);
            kkk=querylj(from,);
            printf("%lld\n",kkk);
        }
    }
}