首页 > 其他 > 详细

P3384 轻重链剖分(树剖模板)

时间:2020-08-02 21:44:06      阅读:76      评论:0      收藏:0      [点我收藏+]

如题,已知一棵包含 NN 个结点的树(连通且无环),每个节点上包含一个数值,需要支持以下操作:

操作 11: 格式: 1\ x\ y\ z1 x y z 表示将树从 xx 到 yy 结点最短路径上所有节点的值都加上 zz。

操作 22: 格式: 2\ x\ y2 x y 表示求树从 xx 到 yy 结点最短路径上所有节点的值之和。

操作 33: 格式: 3\ x\ z3 x z 表示将以 xx 为根节点的子树内所有节点值都加上 zz。

操作 44: 格式: 4\ x4 x 表示求以 xx 为根节点的子树内所有节点值之和

#include<bits/stdc++.h>
using namespace std;
const int maxn=2e5+100;
typedef long long ll;
int n,m,r,mod;
vector<int> g[maxn];

int son[maxn];//重儿子编号 
int id[maxn];//新编号
int fa[maxn];//父亲节点
int cnt;
int dep[maxn];
int size[maxn];//子树大小
int top[maxn];//当前链顶节点
int w[maxn];//初始点权数组
int wt[maxn]; 


//线段树部分

struct node {
    int l,r;
    ll sum;
    ll lazy;
}segTree[maxn*4];
void build (int i,int l,int r) {
    segTree[i].l=l;
    segTree[i].r=r;
    if (l==r) {
        segTree[i].sum=wt[l];
        return;
    }
    int mid=(l+r)>>1;
    build(i<<1,l,mid);
    build(i<<1|1,mid+1,r);
    segTree[i].sum=segTree[i<<1].sum+segTree[i<<1|1].sum;
    segTree[i].sum%=mod;
}
void spread (int i) {
    if (segTree[i].lazy) {
        segTree[i<<1].sum+=segTree[i].lazy*(segTree[i<<1].r-segTree[i<<1].l+1);
        segTree[i<<1].sum%=mod;
        segTree[i<<1|1].sum+=segTree[i].lazy*(segTree[i<<1|1].r-segTree[i<<1|1].l+1);
        segTree[i<<1].sum%=mod;
        segTree[i<<1].lazy+=segTree[i].lazy;
        segTree[i<<1|1].lazy+=segTree[i].lazy;
        segTree[i].lazy=0;
    }
} 
void update (int i,int l,int r,int val) {
    //区间更新
    if (l<=segTree[i].l&&segTree[i].r<=r) {
        segTree[i].sum+=(ll)val*(segTree[i].r-segTree[i].l+1);
        segTree[i].sum%=mod;
        segTree[i].lazy+=val;
        return;
    } 
    spread(i);
    int mid=(segTree[i].l+segTree[i].r)>>1;
    if (l<=mid) 
        update(i<<1,l,r,val);
    if (r>mid)
        update(i<<1|1,l,r,val);
    segTree[i].sum=segTree[i<<1].sum+segTree[i<<1|1].sum;
    segTree[i].sum%=mod;
}
ll query (int i,int l,int r) {
    if (l<=segTree[i].l&&r>=segTree[i].r)
        return segTree[i].sum%mod;
    spread(i);
    int mid=(segTree[i].l+segTree[i].r)>>1;
    ll ans=0;
    if (l<=mid) 
        ans+=query(i<<1,l,r);
    if (r>mid)
        ans+=query(i<<1|1,l,r);
    return ans%mod; 
}

//树剖部分
int qRange (int x,int y) {
    //查询路径权值和 
    int ans=0;
    while (top[x]!=top[y]) {
        if (dep[top[x]]<dep[top[y]]) swap(x,y);
        ans+=query(1,id[top[x]],id[x]);
        //加上x到x所在链的顶端这一段区间的点权和
        ans%=mod;
        x=fa[top[x]];//把x跳到x所在链顶端的那个点的上面一个点 
    }
    if (dep[x]>dep[y]) swap(x,y);
    ans+=query(1,id[x],id[y]);
    return ans%mod;
} 
void upRange (int x,int y,int k) {
    //修改整条路径 
    k%=mod;
    while (top[x]!=top[y]) {
        if (dep[top[x]]<dep[top[y]]) swap(x,y);
        update(1,id[top[x]],id[x],k);
        x=fa[top[x]];
    }
    if (dep[x]>dep[y]) swap(x,y);
    update(1,id[x],id[y],k);
}
int qSon (int x) {
    //查询子树 
    return query(1,id[x],id[x]+size[x]-1);
}
int upSon (int x,int k) {
    //修改子树 
    update(1,id[x],id[x]+size[x]-1,k);
} 
void dfs1 (int x,int f,int deep) {
    dep[x]=deep;
    fa[x]=f;
    size[x]=1;
    int maxson=-1;//记录重儿子的儿子数
    for (int y:g[x]) {
        if (y==f) continue;
        dfs1(y,x,deep+1);
        size[x]+=size[y];
        if (size[y]>maxson) son[x]=y,maxson=size[y];
    } 
}
void dfs2 (int x,int topf) {
    //topf表示当前链的最顶端的节点
    id[x]=++cnt;
     wt[cnt]=w[x];
     top[x]=topf;
     if (!son[x]) return;
     dfs2(son[x],topf);//先处理重儿子
     for (int y:g[x]) {
         if (y==fa[x]||y==son[x]) continue;
         dfs2(y,y);//对于每个轻儿子都有一条属于它自己的链 
     } 
}
int main () {
    scanf("%d%d%d%d",&n,&m,&r,&mod);
    for (int i=1;i<=n;i++) scanf("%d",w+i);
    for (int i=1;i<n;i++) {
        int x,y;
        scanf("%d%d",&x,&y);
        g[x].push_back(y);
        g[y].push_back(x);
    }
    dfs1(r,0,1);
    dfs2(r,r);
    build(1,1,n);
    while (m--) {
        int op;
        scanf("%d",&op);
        if (op==1) {
            int x,y,z;
            scanf("%d%d%d",&x,&y,&z);
            upRange(x,y,z);
        }
        else if (op==2) {
            int x,y;
            scanf("%d%d",&x,&y);
            printf("%d\n",qRange(x,y));
        }
        else if (op==3) {
            int x,y;
            scanf("%d%d",&x,&y);
            upSon(x,y);
        }
        else {
            int x;
            scanf("%d",&x);
            printf("%d\n",qSon(x));
        }
    }
}

 

 

P3384 轻重链剖分(树剖模板)

原文:https://www.cnblogs.com/zhanglichen/p/13420203.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!