[HAOI2015]树上操作

[HAOI2015]树上操作

https://ac.nowcoder.com/acm/problem/19995

思路

比蓝题的树剖模板题少了一个函数...

不会树剖的可以看看洛谷日报以及oiwike.

代码

我复习了一下...
code大概有点长...

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e5+5; 
int w[N];
vector<int>g[N];
int f[N],sz[N],son[N],dep[N];
void dfs(int u,int fa)
{
    sz[u]=1;dep[u]=dep[fa]+1;f[u]=fa;
    for(int v:g[u])
    {
        if(v==fa)    continue;
        dfs(v,u);
        sz[u]+=sz[v];
        if(sz[v]>sz[son[u]])    son[u]=v;
    }
}

int dfn[N],idx[N],top[N],id;
void DFS(int u,int tp)
{
    dfn[u]=++id;
    idx[id]=u;
    top[u]=tp;
    if(!son[u])    return;
    DFS(son[u],tp);
    for(int v:g[u])
    {
        if(!dfn[v])
            DFS(v,v);
    }
}

struct SegTree{
    ll l,r,sum,lazy,len;
}tr[N<<2];

void pushup(int u)
{
    tr[u].sum=tr[u<<1].sum+tr[u<<1|1].sum;
}

void pushdown(int u)
{
    tr[u<<1].sum+=tr[u<<1].len*tr[u].lazy;
    tr[u<<1|1].sum+=tr[u<<1|1].len*tr[u].lazy;
    tr[u<<1].lazy+=tr[u].lazy;
    tr[u<<1|1].lazy+=tr[u].lazy;
    tr[u].lazy=0;
}

void build(int u,int l,int r)
{
    tr[u].l=l,tr[u].r=r,tr[u].len=(r-l+1);
    if(l==r)
    {
        tr[u].sum=w[idx[l]];
        return;
    }
    int mid=(l+r)>>1;
    build(u<<1,l,mid);
    build(u<<1|1,mid+1,r);
    pushup(u);
}

void modify(int u,int L,int R,ll val)
{
    if(L>tr[u].r||R<tr[u].l)    return;
    if(L<=tr[u].l&&tr[u].r<=R)
    {
        tr[u].sum+=tr[u].len*val;
        tr[u].lazy+=val;
        return;
    }
    pushdown(u);
    modify(u<<1,L,R,val);
    modify(u<<1|1,L,R,val);
    pushup(u);
}

ll query(int u,int L,int R)
{
    if(L>tr[u].r||R<tr[u].l)    return 0;
    if(L<=tr[u].l&&tr[u].r<=R)
    {
        return tr[u].sum;
    }
    pushdown(u);
    return query(u<<1,L,R)+query(u<<1|1,L,R);
}

ll ask(int u,int v)//询问u到v的路径权值和.
{
    ll res=0;
    while(top[u]!=top[v])
    {
        if(dep[top[u]]<dep[top[v]])    swap(u,v);
        res+=query(1,dfn[top[u]],dfn[u]);
        u=f[top[u]];
    }
    if(dep[u]<dep[v])    swap(u,v);
    res+=query(1,dfn[v],dfn[u]);
    return res;
} 

int main()
{
    int n,m;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
        scanf("%d",&w[i]);
    for(int i=1;i<n;i++)
    {
        int u,v;
        scanf("%d%d",&u,&v);
        g[u].push_back(v);
        g[v].push_back(u); 
    }
    dfs(1,1);
    DFS(1,1);
    build(1,1,n);
    while(m--)
    {
        int op,x,a;
        scanf("%d",&op);
        if(op==1)
        {
            scanf("%d%d",&x,&a);//x节点权值+a 
            modify(1,dfn[x],dfn[x],a);
        }
        else if(op==2)
        {
            scanf("%d%d",&x,&a);//x的子树权值都+a.
            modify(1,dfn[x],dfn[x]+sz[x]-1,a); 
        }
        else
        {
            scanf("%d",&x);//询问x到根的权值和.
            printf("%lld\n",ask(1,x));
        }
    }
    return 0;
}
lpt的小屋 文章被收录于专栏

我想要一份甜甜的爱情

全部评论
太棒了,我很喜欢的码风
点赞 回复 分享
发布于 2022-10-04 22:39 澳大利亚

相关推荐

11-04 14:10
东南大学 Java
_可乐多加冰_:去市公司包卖卡的
点赞 评论 收藏
分享
10-13 17:47
门头沟学院 Java
wulala.god:图一那个善我面过,老板网上找的题库面的
点赞 评论 收藏
分享
1 收藏 评论
分享
牛客网
牛客企业服务