题目描述

如题，已知一棵包含N个结点的树（连通且无环），每个节点上包含一个数值，需要支持以下操作：

操作1： 格式： 1 x y z 表示将树从x到y结点最短路径上所有节点的值都加上z

操作2： 格式： 2 x y 表示求树从x到y结点最短路径上所有节点的值之和

操作3： 格式： 3 x z 表示将以x为根节点的子树内所有节点值都加上z

操作4： 格式： 4 x 表示求以x为根节点的子树内所有节点值之和

输入输出格式

输入格式：

第一行包含4个正整数N、M、R、P，分别表示树的结点个数、操作个数、根节点序号和取模数（即所有的输出结果均对此取模）。

接下来一行包含N个非负整数，分别依次表示各个节点上初始的数值。

接下来N-1行每行包含两个整数x、y，表示点x和点y之间连有一条边（保证无环且连通）

接下来M行每行包含若干个正整数，每行表示一个操作，格式如下：

操作1： 1 x y z

操作2： 2 x y

操作3： 3 x z

操作4： 4 x

输出格式：

输出包含若干行，分别依次表示每个操作2或操作4所得的结果（对P取模）

输入输出样例

输入样例#1： 复制

5 5 2 24

7 3 7 8 0

1 2

1 5

3 1

4 1

3 4 2

3 2 2

4 5

1 5 1 3

2 1 3
输出样例#1： 复制

2

21

说明

时空限制：1s，128M

数据规模：

对于30%的数据： N≤10,M≤10 N \leq 10, M \leq 10 N≤10,M≤10

对于70%的数据： N≤103,M≤103 N \leq {10}^3, M \leq {10}^3 N≤103,M≤103

对于100%的数据： N≤105,M≤105 N \leq {10}^5, M \leq {10}^5 N≤105,M≤105

（ 其实，纯随机生成的树LCA+暴力是能过的，可是，你觉得可能是纯随机的么233 ）

样例说明：

树的结构如下：

各个操作如下：

故输出应依次为2、21(重要的事情说三遍：记得取模)

解题思路：

　　裸的树链剖分的模板。。。

推荐写法：

#include<bits/stdc++.h>

using namespace std;

#define re register int

#define ll long long

#define INF 0x3f3f3f3f

#define maxn 100009

#define maxm

inline ll read()

{

    ll x=,f=;char ch=getchar();

    while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}

    while(ch>=''&&ch<=''){x=(x<<)+(x<<)+(ll)(ch-'');ch=getchar();}

    return x*f;

}

int id[maxn],rk[maxn],top[maxn],son[maxn],size[maxn],depth[maxn],fa[maxn];

int head[maxn],val[maxn],sum[maxn<<],add[maxn<<];

int n,m,k,ans,tot,cnt,root,base;

struct edge

{

    int to,nxt;

}p[maxn<<];

#define ls(x) x<<1

#define rs(x) x<<1|1

void add_edge(re x,re y)

{

    p[++cnt]={y,head[x]},head[x]=cnt;

}

void dfs1(int u,int father)

{

    depth[u]=depth[father]+,fa[u]=father,size[u]=;

    for(int i=head[u];i;i=p[i].nxt)

    {

        int v=p[i].to;

        if(v==father)    continue;

        dfs1(v,u);

        size[u]+=size[v];

        if(size[v]>size[son[u]])

            son[u]=v;

    }

}

void dfs2(int u,int father)

{

    top[u]=father,id[u]=++tot,rk[tot]=u;

    if(!son[u])    return ;

    dfs2(son[u],father);

    for(int i=head[u];i;i=p[i].nxt)

    {

        int v=p[i].to;

        if(v!=son[u]&&v!=fa[u])

            dfs2(v,v);

    }

}

int mod(int x)

{

    return x>=base?x%base:x;

}

void push_up(int x)

{

    sum[x]=sum[ls(x)]+sum[rs(x)];

    sum[x]=mod(sum[x]);

}

void built(int x,int l,int r)

{

    if(l==r)

    {

        sum[x]=val[rk[l]];

        return ;

    }

    int mid=(l+r)>>;

    built(ls(x),l,mid);

    built(rs(x),mid+,r);

    push_up(x);

}

void pass(int x,int l,int r,int k)

{

    add[x]+=k,add[x]=mod(add[x]);

    sum[x]+=(r-l+)*k,sum[x]=mod(sum[x]);

}

void push_down(int x,int l,int r)

{

    if(!add[x]) return ;

    int mid=(l+r)>>;

    pass(ls(x),l,mid,add[x]);

    pass(rs(x),mid+,r,add[x]);

    add[x]=;

}

int Ask(int x,int l,int r,int nl,int nr)

{

    int res=;

    if(nl<=l&&r<=nr)

        return sum[x];

    push_down(x,l,r);

    int mid=(l+r)>>;

    if(nl<=mid)

        res+=Ask(ls(x),l,mid,nl,nr),res=mod(res);

    if(mid<nr)

        res+=Ask(rs(x),mid+,r,nl,nr),res=mod(res);

    push_up(x);

    return res;

}

void Add(int x,int l,int r,int nl,int nr,int k)

{

    if(nl<=l&&r<=nr)

    {

        add[x]+=k,add[x]=mod(add[x]);

        sum[x]+=(r-l+)*k,sum[x]=mod(sum[x]);

        return ;

    }

    push_down(x,l,r);

    int mid=(l+r)>>;

    if(nl<=mid)

        Add(ls(x),l,mid,nl,nr,k);

    if(mid<nr)

        Add(rs(x),mid+,r,nl,nr,k);

    push_up(x);

} 

void Work1(int x,int y,int z)

{

    int fx=top[x],fy=top[y];

    while(fx!=fy)

    {

        if(depth[fx]>depth[fy])

            swap(x,y),swap(fx,fy);

        Add(,,tot,id[fy],id[y],z);

        y=fa[fy],fy=top[y];

    }

    if(id[x]>id[y])

        swap(x,y);

    Add(,,tot,id[x],id[y],z);

}

int Work2(int x,int y)

{

    int res=,fx=top[x],fy=top[y];

    while(fx!=fy)

    {

        if(depth[fx]>depth[fy])    swap(x,y),swap(fx,fy);

        res+=Ask(,,tot,id[fy],id[y]),res=mod(res);

        y=fa[fy],fy=top[y];

    }

    if(id[x]>id[y])    swap(x,y);

    res+=Ask(,,tot,id[x],id[y]),res=mod(res);

}

int main()

{

//    freopen("1.in","r",stdin);

//    freopen("1.out","w",stdout);

    n=read(),m=read(),root=read(),base=read();

    for(re i=;i<=n;i++)

        val[i]=read();

    for(re i=;i<n;i++)

    {

        re x=read(),y=read();

        add_edge(x,y),add_edge(y,x);

    }

    dfs1(root,root),dfs2(root,root);

    built(,,tot);

    for(re i=;i<=m;i++)

    {

        re opt=read(),x,y,z;

        if(opt==)

        {

            x=read(),y=read(),z=read();

            Work1(x,y,z);

        }

        if(opt==)

        {

            x=read(),y=read();

            printf("%d\n",Work2(x,y));

        }

        if(opt==)

        {

            x=read(),z=read();

            Add(,,tot,id[x],id[x]+size[x]-,z);

        }

        if(opt==)

        {

            x=read();

            printf("%d\n",Ask(,,tot,id[x],id[x]+size[x]-));

        }

    }

    fclose(stdin);

    fclose(stdout);

    return ;

}
#include<bits/stdc++.h>

using namespace std;

#define re register int

#define ll long long

#define INF 0x3f3f3f3f

#define maxn 100009

#define maxm

#define ls(x) x<<1

#define rs(x) x<<1|1

inline ll read()

{

    ll x=,f=;char ch=getchar();

    while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}

    while(ch>=''&&ch<=''){x=(x<<)+(x<<)+(ll)(ch-'');ch=getchar();}

    return x*f;

}

ll base;

int n,m,k,ans,tot,cnt,root;

struct edge

{

    int to,nxt;

}p[maxn<<];

int head[maxn],id[maxn],fa[maxn],depth[maxn],size[maxn],son[maxn],rk[maxn],top[maxn];

ll sum[maxn<<],add[maxn<<],val[maxn];

void add_edge(int x,int y)

{

    p[++cnt].to=y,p[cnt].nxt=head[x],head[x]=cnt;

}

void dfs1(int u,int father)

{

    depth[u]=depth[father]+,size[u]=,fa[u]=father;

    for(int i=head[u];i;i=p[i].nxt)

    {

        int v=p[i].to;

        if(v==father)    continue;

        dfs1(v,u);

        size[u]+=size[v];

        if(size[v]>size[son[u]])

            son[u]=v;

    }

}

void dfs2(int u,int father)

{

    top[u]=father;

    id[u]=++tot;

    rk[tot]=u;

    if(!son[u])

        return ;

    dfs2(son[u],father);

    for(int i=head[u];i;i=p[i].nxt)

    {

        int v=p[i].to;

        if(v!=son[u]&&v!=fa[u])

            dfs2(v,v);

    }

}

void push_up(int x)

{

    sum[x]=(sum[ls(x)]+sum[rs(x)])%base;

}

void built(int x,int l,int r)

{

    if(l==r)

    {

        sum[x]=val[rk[l]];

        return ;

    }

    int mid=(l+r)>>;

    built(ls(x),l,mid);

    built(rs(x),mid+,r);

    push_up(x);

}

void pass(int x,int l,int r,ll k)

{

    add[x]+=k,add[x]%=base;

    sum[x]+=(r-l+)*k,sum[x]%=base;

}

void push_down(int x,int l,int r)

{

    if(!add[x])    return ;

    int mid=(l+r)>>;

    pass(ls(x),l,mid,add[x]);

    pass(rs(x),mid+,r,add[x]);

    add[x]=;

}

ll Query(int x,int l,int r,int nl,int nr)

{

    ll res=;

    if(nl<=l&&r<=nr)

        return sum[x];

    push_down(x,l,r);

    int mid=(l+r)>>;

    if(nl<=mid)

        res+=Query(ls(x),l,mid,nl,nr),res%=base;

    if(mid<nr)

        res+=Query(rs(x),mid+,r,nl,nr),res%=base;

    push_up(x);

    return res%base;

}

void Add(int x,int l,int r,int nl,int nr,ll k)

{

    if(nl<=l&&r<=nr)

    {

        add[x]+=k,add[x]%=base;

        sum[x]+=(r-l+)*k,sum[x]%=base;

        return ;

    }

    push_down(x,l,r);

    int mid=(l+r)>>;

    if(nl<=mid)

        Add(ls(x),l,mid,nl,nr,k);

    if(mid<nr)

        Add(rs(x),mid+,r,nl,nr,k);

    push_up(x);

}

void Work1(int x,int y,ll z)

{

    int fx=top[x],fy=top[y];

    while(fx!=fy)

    {

        if(depth[fx]>=depth[fy])

        {

            Add(,,tot,id[fx],id[x],z);

            x=fa[fx],fx=top[x];

        }

        else

        {

            Add(,,tot,id[fy],id[y],z);

            y=fa[fy],fy=top[y];

        }

    }

    if(id[x]<=id[y])

        Add(,,tot,id[x],id[y],z);

    else

        Add(,,tot,id[y],id[x],z);

}

ll Work2(int x,int y)

{

    ll ans=;

    int fx=top[x],fy=top[y];

    while(fx!=fy)

    {

        if(depth[fx]>=depth[fy])

        {

            ans+=Query(,,tot,id[fx],id[x]);

            x=fa[fx];fx=top[x];

        }

        else

        {

            ans+=Query(,,tot,id[fy],id[y]);

            y=fa[fy];fy=top[y];

        }

    }

    if(id[x]<=id[y])

        ans+=Query(,,tot,id[x],id[y]),ans%=base;

    else

        ans+=Query(,,tot,id[y],id[x]),ans%=base;

    return ans%base;

}

void Work3(int x,int y)

{

    Add(,,tot,id[x],id[x]+size[x]-,y);

}

ll Work4(int x)

{

    ll res=Query(,,tot,id[x],id[x]+size[x]-);

    return res%base;

}

int main()

{

     freopen("1.in","r",stdin);

     freopen("1.out","w",stdout);

    n=read(),m=read(),root=read(),base=read();

    for(int i=;i<=n;i++)

        val[i]=read()%base;

    for(int i=;i<n;i++)

    {

        int  x=read(),y=read();

        add_edge(x,y),add_edge(y,x);

    }

    dfs1(root,root),dfs2(root,root);

    built(,,tot);

    for(int i=;i<=m;i++)

    {

        int opt=read(),x,y;

        ll z;

        if(opt==)

        {

            x=read(),y=read(),z=read();

            Work1(x,y,z);

        }

        if(opt==)

        {

            x=read(),y=read();

            printf("%lld\n",Work2(x,y));

        }

        if(opt==)

        {

            x=read(),z=read();

            Work3(x,z);

        }

        if(opt==)

        {

            x=read();

            printf("%lld\n",Work4(x));

        }

    }

    fclose(stdin);

    fclose(stdout);

    return ;

}

Luogu 3384 【模板】树链剖分的更多相关文章

  1. [luogu P3384] [模板]树链剖分

    [luogu P3384] [模板]树链剖分 题目描述 如题,已知一棵包含N个结点的树(连通且无环),每个节点上包含一个数值,需要支持以下操作: 操作1: 格式: 1 x y z 表示将树从x到y结点 ...

  2. luoguP3384 [模板]树链剖分

    luogu P3384 [模板]树链剖分 题目 #include<iostream> #include<cstdlib> #include<cstdio> #inc ...

  3. 【Luogu P3384】树链剖分模板

    树链剖分的基本思想是把一棵树剖分成若干条链,再利用线段树等数据结构维护相关数据,可以非常暴力优雅地解决很多问题. 树链剖分中的几个基本概念: 重儿子:对于当前节点的所有儿子中,子树大小最大的一个儿子就 ...

  4. [洛谷P3384] [模板] 树链剖分

    题目传送门 显然是一道模板题. 然而索引出现了错误,狂wa不止. 感谢神犇Dr_J指正.%%%orz. 建线段树的时候,第44行. 把sum[p]=bv[pos[l]]%mod;打成了sum[p]=b ...

  5. 模板 树链剖分BFS版本

    //点和线段树都从1开始 //边使用vector vector<int> G[maxn]; ],num[maxn],iii[maxn],b[maxn],a[maxn],top[maxn], ...

  6. 洛谷3384&bzoj1036树链剖分

    值得注意的是: 一个点的子树是存在一起的...也就是说我们修改子树的时候只用... /********************************************************* ...

  7. P3384 [模板] 树链剖分

    #include <bits/stdc++.h> using namespace std; typedef long long ll; int n, m, rt, mod, cnt, to ...

  8. 『题解』洛谷P3384 【模板】树链剖分

    Problem Portal Portal1: Luogu Description 如题,已知一棵包含\(N\)个结点的树(连通且无环),每个节点上包含一个数值,需要支持以下操作: 操作\(1\): ...

  9. 树链剖分详解（洛谷模板 P3384）

    洛谷·[模板]树链剖分 写在前面 首先,在学树链剖分之前最好先把 LCA.树形DP.DFS序 这三个知识点学了 emm还有必备的 链式前向星.线段树 也要先学了. 如果这三个知识点没掌握好的话,树链剖 ...

随机推荐

  1. python学习之 getattr vs __get__ vs __getattr __ vs __getattribute__ vs __getitem__

    1. getattr.setattr.hasattr getattr比较常用,与setattr和hasattr一起出现,他们也是最容易理解的,下面是他的用法: class Profile(): nam ...

  2. 解决select下拉框禁用（设置disabled属性）,后台获取值为空

    如果下拉框设置disabled属性后,提交表单到后台,后台获取的下拉框的值为空,以下有三种解决获取不到下拉框选项值的方法. 有下拉框html如:<select name="select ...

  3. 利用request和re抓取猫眼电影排行

    import requests import re import time def get_one_page(url): headers = { 'User-Agent': 'Mozilla/5.0 ...

  4. 【译】索引进阶（十一）：SQL SERVER中的索引碎片【上篇】

    原文链接:传送门. 第十章节我们分析了索引的内部结构.有了这些关于索引结构的知识,我们便可以分析索引碎片了:其产生的原因,如何防止,以及何时可以不去关注它们. 一些背景知识 / 复习 以下知识对于理解 ...

  5. Java位运算原理及使用讲解

    前言日常开发中位运算不是很常用,但是巧妙的使用位运算可以大量减少运行开销,优化算法.举个例子,翻转操作比较常见,比如初始值为1,操作一次变为0,再操作一次变为1.可能的做法是使用三木运算符,判断原始值 ...

  6. git知识总结-3.gitignore文件说明

    1.前言 一般我们总会有些文件无需纳入 Git 的管理,也不希望它们总出现在未跟踪文件列表. 通常都是些自动生成的文件,比如日志文件,或者编译过程中创建的临时文件等. 在这种情况下,我们可以创建一个名 ...

  7. day17——其他内置函数

    zip函数: print(list(zip(('a','b','c'),(1,2,3)))) p={'name':'alex','age':18,'gender':'none'} print(list ...

  8. 高可用Redis(九)：Redis Sentinel

    1.主从复制高可用的问题 主从复制高可用的作用 1.为master提供备份,当master宕机时,slave有完整的备份数据 2.对master实现分流,实现读写分离 但是主从架构有一个问题 1.如果 ...

  9. spring-cloud-hystrix-dasboard服务调用监控

    除了隔离依赖服务的调用以外,hystrix还提供了准实时的调用监控(hystrix dashboard),hystrxi会持续的记录所有通过hyxtrix发起的请求的执行信息,并以统计报表和图形的形式 ...

  10. Linux Centos7.5中的RocketMQ集群部署

    系统环境 Docker > centos7.5 此镜像已经安装了jdk1.8和maven3.6.0 如果你想知道这个基础镜像的具体情况, 参考此文: https://www.cnblogs.co ...