题目描述

如题，已知一棵包含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. Keras的一些功能函数

    摘自: https://www.cnblogs.com/Anita9002/p/8136357.html 1.模型的信息提取 # 节点信息提取 config = model.get_config() ...

  2. vue项目的常用配置代码

    { // 针对vue的格式化配置----依赖eslint.prettier.vetur等插件 // 强制单引号 "prettier.singleQuote": true, &quo ...

  3. Ubuntu通过ADB连接手机

    参考 ubuntu14.04 下android studio连接手机 安装 adb $sudo apt install adb $sudo lsusb 得到ID为 Bus Device : ID 12 ...

  4. MySQL中使用group_concat()函数数据字符过长报错的问题解决方法

    最近在办公软件项目,在开发权限指标遇到一个问题:我们系统的一些逻辑处理是用存储过程实现的,但是有一天客户反馈说权限指标分配报错,查了分配的权限数据牵扯到的数据权限基础资源,没有问题.权限指标分配的存储 ...

  5. 微信商户/H5支付申请 被拒原因：网站存在不实内容或不安全信息

    看到这张图片,申请H5支付的朋友估计是崩溃的,被拒的原因是网站存在不实内容或不安全信息. 解决方法有两种: 针对不安全内容,则是微信支付审核团队认为你的网站存在安全漏洞的,需要修复. 针对不实内容,则 ...

  6. Debian 命令行方式配置网络

    一.对于有线网络,如果默认没有安装图形界面,进入了 multi-user.target中时,是没有使用NetworkManager管理网络的,此时需要手动配置才能上网 首先得到网卡名称:ip addr ...

  7. 学习笔记——单片机简介 & 点亮LED & 流水灯 & 电路基础【更新Ing】

    视频地址:https://www.bilibili.com/video/av10765766 超详细!!!!!! 单片机内部三大资源 [资源:单片机可提供使用的东西] FLASH 可以重复擦写 断电后 ...

  8. go语言使用xpath

    1.导包 gopm get -g -v github.com/lestrrat-go/libxml2 2.使用示例 func ExampleHTML() { res, err := http.Get( ...

  9. gym 102059A 树链剖分后odt维护区间

    题意 一棵树 多次修改,每次修改一个点到根的所有边的颜色,并询问现在有哪些颜色染了恰好$m$条边 题解: 稍加思考可以知道,从某个点到根节点的颜色数,均摊复杂度很低,因此,可以考虑珂朵莉树维护重链剖分 ...

  10. lambda-基于谓词筛选值序列

    此方法通过使用延迟执行实现. 即时返回值为一个对象,该对象存储执行操作所需的所有信息. 只有通过直接调用对象的 GetEnumerator 方法或使用 Visual C# 中的 foreach(或 V ...