洛谷树剖模板题 P3384 | 树链剖分
对于以u为根的子树,后代节点的dfn显然比他的dfn大,我们可以记录一下回溯到u的dfn,显然这两个dfn构成了一个连续区间,代表u及u的子树
剩下的就和树剖一样了
#include<cstdio>
#include<algorithm>
#include<cstring>
#define N 100010
typedef long long ll;
using namespace std;
int ecnt,head[N],son[N],fa[N],sz[N],n,m,r,deep[N],pos[N],indx[N],tot,op,back[N],top[N];
ll val[N],P;
struct adj
{
int nxt,v;
}e[*N];
struct node
{
int l,r;
ll sum,lz;
}t[*N];
void add(int u,int v)
{
e[++ecnt].v=v;
e[ecnt].nxt=head[u];
head[u]=ecnt;
e[++ecnt].v=u;
e[ecnt].nxt=head[v];
head[v]=ecnt;
}
void dfs1(int x,int father,int dep)
{
deep[x]=dep,fa[x]=father,sz[x]=;
for (int i=head[x];i;i=e[i].nxt)
{
int v=e[i].v;
if (fa[x]==v) continue;
dfs1(v,x,dep+);
sz[x]+=sz[v];
if (sz[v]>sz[son[x]]) son[x]=v;
}
}
void dfs2(int x,int TOP)
{
top[x]=TOP;
pos[x]=++tot;
indx[tot]=x;
if (son[x]) dfs2(son[x],TOP);
for (int i=head[x];i;i=e[i].nxt)
{
int v=e[i].v;
if (v==fa[x] || v==son[x]) continue;
dfs2(v,v);
}
back[x]=tot;
}
void pushdown(int p)
{
if (t[p].l==t[p].r || !t[p].lz) return;
int w=t[p].lz;
t[p<<].sum=(t[p<<].sum+w*(t[p<<].r-t[p].l+)%P)%P;
t[p<<|].sum=(t[p<<|].sum+w*(t[p<<|].r-t[p<<|].l+)%P)%P;
t[p<<].lz+=w;
t[p<<].lz%=P;
t[p<<|].lz+=w;
t[p<<|].lz%=P;
t[p].lz=;
}
void pushup(int p)
{
t[p].sum=(t[p<<].sum+t[p<<|].sum)%P;
}
void build(int p,int l,int r)
{
t[p].l=l,t[p].r=r,t[p].lz=;
if (l!=r)
{
int mid=l+r>>;
build(p<<,l,mid);
build(p<<|,mid+,r);
pushup(p);
}
else
t[p].sum=val[indx[l]]%P;
}
void modify(int p,int l,int r,int w)
{
if (t[p].l==l && t[p].r==r)
{
t[p].sum+=w*(t[p].r-t[p].l+);
t[p].lz+=w;
return ;
}
pushdown(p);
int mid=t[p].l+t[p].r>>;
if (r<=mid) modify(p<<,l,r,w);
else if (l>mid) modify(p<<|,l,r,w);
else modify(p<<,l,mid,w),modify(p<<|,mid+,r,w);
pushup(p);
}
ll query(int p,int l,int r)
{
if (t[p].l==l && t[p].r==r)
return t[p].sum%P;
int mid=t[p].l+t[p].r>>;
pushdown(p);
if (r<=mid) return query(p<<,l,r)%P;
if (l>mid) return query(p<<|,l,r)%P;
return (query(p<<,l,mid)+query(p<<|,mid+,r))%P;
}
void pathInc(int u,int v,int w)
{
while (top[u]!=top[v])
{
if (deep[top[u]]<deep[top[v]]) swap(u,v);
modify(,pos[top[u]],pos[u],w);
u=fa[top[u]];
}
if (deep[u]>deep[v]) swap(u,v);
modify(,pos[u],pos[v],w);
}
ll pathQuery(int u,int v)
{
ll ret=;
while (top[u]!=top[v])
{
if (deep[top[u]]<deep[top[v]]) swap(u,v);
ret=(ret+query(,pos[top[u]],pos[u]))%P;
u=fa[top[u]];
}
if (deep[u]>deep[v]) swap(u,v);
return (ret+query(,pos[u],pos[v]))%P;
}
int main()
{
scanf("%d%d%d%lld",&n,&m,&r,&P);
for (int i=;i<=n;i++)
scanf("%lld",&val[i]);
for (int i=,u,v;i<n;i++)
scanf("%d%d",&u,&v),add(u,v);
dfs1(r,,);
dfs2(r,r);
build(,,n);
for (int i=,x,y,z;i<=m;i++)
{
scanf("%d",&op);
if (op==)
scanf("%d%d%d",&x,&y,&z),pathInc(x,y,z);
else if (op==)
scanf("%d%d",&x,&y),printf("%lld\n",pathQuery(x,y));
else if (op==)
scanf("%d%d",&x,&z),modify(,pos[x],back[x],z);
else scanf("%d",&x),printf("%lld\n",query(,pos[x],back[x]));
}
return ;
}
洛谷树剖模板题 P3384 | 树链剖分的更多相关文章
- 洛谷 P3373 【模板】线段树 2
洛谷 P3373 [模板]线段树 2 洛谷传送门 题目描述 如题,已知一个数列,你需要进行下面三种操作: 将某区间每一个数乘上 xx 将某区间每一个数加上 xx 求出某区间每一个数的和 输入格式 第一 ...
- 洛谷P3372 【模板】线段树 1
P3372 [模板]线段树 1 153通过 525提交 题目提供者HansBug 标签 难度普及+/提高 提交 讨论 题解 最新讨论 [模板]线段树1(AAAAAAAAA- [模板]线段树1 洛谷 ...
- 洛谷P3373 【模板】线段树 2
P3373 [模板]线段树 2 47通过 186提交 题目提供者HansBug 标签 难度提高+/省选- 提交 讨论 题解 最新讨论 为啥WA(TAT) 题目描述 如题,已知一个数列,你需要进行 ...
- 洛谷 P3373 【模板】线段树 2 解题报告
P3373 [模板]线段树 2 题目描述 如题,已知一个数列,你需要进行下面三种操作: 1.将某区间每一个数乘上\(x\) 2.将某区间每一个数加上\(x\) 3.求出某区间每一个数的和 输入输出格式 ...
- 洛谷 P3372 【模板】线段树 1
P3372 [模板]线段树 1 题目描述 如题,已知一个数列,你需要进行下面两种操作: 1.将某区间每一个数加上x 2.求出某区间每一个数的和 输入输出格式 输入格式: 第一行包含两个整数N.M,分别 ...
- 洛谷—— P3372 【模板】线段树 1
P3372 [模板]线段树 1 题目描述 如题,已知一个数列,你需要进行下面两种操作: 1.将某区间每一个数加上x 2.求出某区间每一个数的和 输入输出格式 输入格式: 第一行包含两个整数N.M,分别 ...
- 洛谷——P3373 【模板】线段树 2&& B 数据结构
P3373 [模板]线段树 2 题目描述 如题,已知一个数列,你需要进行下面三种操作: 1.将某区间每一个数乘上x 2.将某区间每一个数加上x 3.求出某区间每一个数的和 线段树维护区间乘法 1.如何 ...
- 洛谷 P3373 【模板】线段树 2 题解
P3373 [模板]线段树 2 题目描述 如题,已知一个数列,你需要进行下面三种操作: 1.将某区间每一个数乘上x 2.将某区间每一个数加上x 3.求出某区间每一个数的和 输入格式 第一行包含三个整数 ...
- 洛谷P3373 【模板】线段树 2 && P2023 [AHOI2009]维护序列——题解
题目传送: P3373 [模板]线段树 2 P2023 [AHOI2009]维护序列 该题较传统线段树模板相比多了一个区间乘的操作.一提到线段树的区间维护问题,就自然想到了“懒标记”:为了降低时间复 ...
随机推荐
- AMD、CMD、Common规范及对比
https://blog.csdn.net/bluesky1215/article/details/71081780 1.名词解释 AMD:Asynchronous Modules Definiti ...
- 树梅派3B kali2.0 启用SSH进行远程登录
工具/原料 kali 2.0 ssh SSH连接工具(XShell)or PUTTY vi /etc/ssh/sshd_config 将#PasswordAuthentication no的注释去掉, ...
- Windows API窗口绘图程序设计
任务目标 设计一个简单的Windows 窗口程序,在程序窗口内任意位置按下鼠标左键,可绘制范围在10-100之间随机大小的正方形.并且显示的正方形用红色填充. 效果图 小结 程序先是触发鼠标左键点击事 ...
- ES6笔记02-箭头函数
eg1:// ES5 匿名函数 var total = values.reduce(function (a, b) { return a + b; }, 0); // ES6 匿名函数 var tot ...
- 使用MyEclipse/Eclipse修改项目名称报Can't convert argument: null!
报错: java.lang.IllegalArgumentException: Can't convert argument: null! 方法/步骤 报错原因:使用MyEclipse修改项目名 ...
- ECSHOP快递物流单号查询插件
本ECSHOP快递物流单号跟踪插件提供国内外近2000家快递物流订单单号查询服务例如申通快递.顺丰快递.圆通快递.EMS快递.汇通快递.宅急送快递.德邦物流.百世快递.汇通快递.中通快递.天天快递等知 ...
- PLC状态机编程第六篇-优化PLC程序生成
还记得第一篇博客中,我们在状态机中手写上升沿来处理有别于传统的一键启停程序,那个手写的上升沿就是优化手段.stateflow状态机是带事件的,事件本身支持上升沿和下降沿等事件,在这里,如果我们选择用事 ...
- PyQuery网页解析库
from pyquery import PyQuery as pq 字符串初始化: doc = pq(html) URL初始化:doc = pq(url = "···") 文件初始 ...
- realloc函数的用法
realloc(void *__ptr, size_t __size):更改已经配置的内存空间,即更改由malloc()函数分配的内存空间的大小. 如果将分配的内存减少,realloc仅仅是改变索引的 ...
- 陌生又熟悉的数据库之ID增加
当我们设计一张表时,通常为了保证记录的唯一性,会为表增加一个ID字段,生成记录时ID自动加一