题目描述

如题，已知一棵包含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. 记一场与 cookie 的相遇

    简介: cookie 翻译过来为 “小甜点,一种酥性甜饼干,很美味的...”,咳咳,打住!我们这里说的是 “甜点” 文件,它是浏览器储存在用户电脑上的一小段纯文本格式的文件. 由于 http 是一种无 ...

  2. CPU或以太网模块重启DHCP请求

    最近,带了一个实习生,一天: 他:师傅,PLC如何首次分配IP地址啊? 我:不是教过你了吗? 他:我怎么用BOOTP软件,半天没有反应啊? 我:..... 原来,他做实验的CPU已经被分配IP地址了, ...

  3. 转：jsp与servlet的区别与联系

    jsp与servlet的区别与联系 - gsyabc - 博客园https://www.cnblogs.com/sanyouge/p/7325656.html jsp和servlet的区别和联系:1. ...

  4. VS Code保存使用项目Eslint规则格式化代码

    One.文件-首选项-设置-选择项目                  Two.打开右上角JSON设置            Three. 插入以下代码 { "eslint.options& ...

  5. HDU - 1062

    格式错误2遍:没考虑到连续两个空格的情况,遇到空格最后要输出这个空格,因为题目只需要转换单词. 另外,开cin,cout加速要注意读入不能用scanf,printf,puts,getchar这些.ge ...

  6. 如何在sublime安装ctags函数追踪插件

    一.这个插件是需要下载执行软件才能够操作的,下载地址http://prdownloads.sourceforge.net/ctags/ctags58.zip 然后解压到自己喜欢的目录下,我是安装到 D ...

  7. redis缓存雪崩、缓存穿透、数据库和redis数据一致性

    一.缓存雪崩 回顾一下我们为什么要用缓存(Redis):减轻数据库压力或尽可能少的访问数据库. 在前面学习我们都知道Redis不可能把所有的数据都缓存起来(内存昂贵且有限),所以Redis需要对数据设 ...

  8. Shell中字符串的切割、拼接、比较、替换

    [截取] 一.Linux shell 截取字符变量的前8位,有方法如下: expr substr “$a” 1 8 : 二.按指定的字符串截取 第一种方法: ${varible##*string} # ...

  9. mybatis调用oracle存储过程的几个参考例子

    首先写一个存储过程: create or replace procedure p_syn_equipment_20161205 is sqlstr ); begin --清空表 sqlstr := ' ...

  10. 电脑移动后WIFI连接失败解决方法

    1.现象原因 经常会发现将自己的电脑带到不同的地方后连接附近WIFI失败的现象,这是什么原因造成的了,觉得明明之前还有连过这个无线,密码都是正确的,无线连接的图标显示一个大大大的感叹号!  像下面一样 ...