P4211 [LNOI2014]LCA(离线 + 在线 做法)

P4211 [LNOI2014]LCA

有一棵根节点为111树,有mmm次询问,每次给定l,r,zl, r, zl,r,z,输出∑i=lrdep[lca(i,z)]\sum\limits_{i = l} ^{r} dep[lca(i, z)]i=lrdep[lca(i,z)]

乍一看这题好像无从下手,仔细想想lca(i,z)lca(i, z)lca(i,z)有何性质,不难发现lca(i,z)lca(i, z)lca(i,z)一定是在1−>z1->z1>z的路径上的,

那么我们的答案求解就可以转换为每次对i−>1i->1i>1路径上的点加一,求z−>1z->1z>1上路径上的点的和即为我们要求的答案。

考虑如何优化这一过程:

离线求解:

将每个询问(l,r,z)(l, r, z)(l,r,z)拆分成(1,l−1,z),(1,r,z)(1, l - 1, z), (1, r, z)(1,l1,z),(1,r,z),那么答案就为(1,r,z)−(1,l−1,z)(1, r, z) - (1, l - 1, z)(1,r,z)(1,l1,z)了,

我们将所有拆分后询问按照rrr排个序,从小到大每次插入点iii,如果当前插入的点iii等于每个询问的rrr则统计一次答案,

两点之间点权的信息更新可以考虑树剖,注意答案的加减,最后输出答案即可,整体复杂度O(nlog⁡nlog⁡n)O(n \log n \log n)O(nlognlogn)

#include <bits/stdc++.h>
#define mid (l + r >> 1)
#define lson rt << 1, l, mid
#define rson rt << 1 | 1, mid + 1, r
#define ls rt << 1
#define rs rt << 1 | 1using namespace std;const int N = 1e5 + 10, mod = 201314;int head[N], to[N], nex[N], cnt = 1;int fa[N], dep[N], sz[N], son[N], rk[N], id[N], top[N], tot;int sum[N << 2], lazy[N << 2], ans[N], n, m;struct Res {int u, v, id;bool operator < (const Res &t) {return u < t.u;}
}a[N];void add(int x, int y) {to[cnt] = y;nex[cnt] = head[x];head[x] = cnt++;
}void dfs1(int rt, int f) {fa[rt] = f, dep[rt] = dep[f] + 1, sz[rt] = 1;for (int i = head[rt]; i; i = nex[i]) {if (to[i] == f) {continue;}dfs1(to[i], rt);sz[rt] += sz[to[i]];if (!son[rt] || sz[son[rt]] < sz[to[i]]) {son[rt] = to[i];}}
}void dfs2(int rt, int tp) {top[rt] = tp, rk[++tot] = rt, id[rt] = tot;if (!son[rt]) {return ;}dfs2(son[rt], tp);for (int i = head[rt]; i; i = nex[i]) {if (to[i] == fa[rt] || to[i] == son[rt]) {continue;}dfs2(to[i], to[i]);}
}void push_up(int rt) {sum[rt] = (sum[ls] + sum[rs]) % mod;
}void push_down(int rt, int l, int r) {if (lazy[rt]) {lazy[ls] = (lazy[ls] + lazy[rt]) % mod, lazy[rs] = (lazy[rs] + lazy[rt]) % mod;sum[ls] = (sum[ls] + 1ll * (mid - l + 1) * lazy[rt]) % mod;sum[rs] = (sum[rs] + 1ll * (r - mid) * lazy[rt]) % mod;lazy[rt] = 0;}
}void update(int rt, int l, int r, int L, int R, int v) {if (l >= L && r <= R) {lazy[rt] = (lazy[rt] + 1) % mod;sum[rt] = (sum[rt] + r - l + 1) % mod;return ;}push_down(rt, l, r);if (L <= mid) {update(lson, L, R, v);}if (R > mid) {update(rson, L, R, v);}push_up(rt);
}int query(int rt, int l, int r, int L, int R) {if (l >= L && r <= R) {return sum[rt];}push_down(rt, l, r);int ans = 0;if (L <= mid) {ans = query(lson, L, R);}if (R > mid) {ans = (ans + query(rson, L, R)) % mod;}return ans;
}void update(int rt) {while (rt) {update(1, 1, n, id[top[rt]], id[rt], 1);rt = fa[top[rt]];}
}int query(int rt) {int ans = 0;while (rt) {ans = (ans + query(1, 1, n, id[top[rt]], id[rt])) % mod;rt = fa[top[rt]];}return ans;
}int main() {// freopen("in.txt", "r", stdin);// freopen("out.txt", "w", stdout);// ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);scanf("%d %d", &n, &m);for (int i = 2, fa; i <= n; i++) {scanf("%d", &fa);add(fa + 1, i);}dfs1(1, 0);dfs2(1, 1);for (int i = 1, l, r, u; i <= m; i++) {scanf("%d %d %d", &l, &r, &u);l++, r++, u++;a[i * 2 - 1] = {l - 1, u, i * 2 - 1};a[i * 2] = {r, u, i * 2};}sort(a + 1, a + 1 + 2 * m);int cur = 1;while (cur <= 2 * m && a[cur].u == 0) {cur++;    }for (int i = 1; i <= n; i++) {update(i);while (cur <= 2 * m && a[cur].u == i) {int res = query(a[cur].v);if (a[cur].id & 1) {ans[a[cur].id + 1 >> 1] = (ans[a[cur].id + 1 >> 1] + mod - res) % mod;}else {ans[a[cur].id >> 1] = (ans[a[cur].id >> 1] + res) % mod;}cur++;}}for (int i = 1; i <= m; i++) {printf("%d\n", ans[i]);}return 0;
}

在线求解:

考虑主席树,第iii棵主席树为前i−1i - 1i1棵主席树的信息 + 节点编号为iii的点到根节点权值都+1+1+1的信息。

不难想到,最后我们的答案就是在[l,r][l, r][l,r]区间的主席树上查询在1−>z1->z1>z之间的点的权值和。

主席树上打lazylazylazy,空间按道理应该是O(nlog⁡nlog⁡n)O(n \log n \log n)O(nlognlogn)级别的,支持在线求解时间复杂度也是O(nlog⁡nlog⁡n)O(n \log n \log n)O(nlognlogn)级别的。

#include <bits/stdc++.h>using namespace std;const int N = 5e4 + 10, mod = 201314;int head[N], to[N], nex[N], cnt = 1;int fa[N], son[N], sz[N], top[N], dep[N], rk[N], id[N], tot;int root[N], ls[N * 200], rs[N * 200], sum[N * 200], lazy[N * 200], n, m, num;void add(int x, int y) {to[cnt] = y;nex[cnt] = head[x];head[x] = cnt++;
}void dfs1(int rt, int f) {fa[rt] = f, dep[rt] = dep[f] + 1, sz[rt] = 1;for (int i = head[rt]; i; i = nex[i]) {if (to[i] == f) {continue;}dfs1(to[i], rt);sz[rt] += sz[to[i]];if (!son[rt] || sz[son[rt]] < sz[to[i]]) {son[rt] = to[i];}}
}void dfs2(int rt, int tp) {top[rt] = tp, rk[++tot] = rt, id[rt] = tot;if (!son[rt]) {return ;}dfs2(son[rt], tp);for (int i = head[rt]; i; i = nex[i]) {if (to[i] == fa[rt] || to[i] == son[rt]) {continue;}dfs2(to[i], to[i]);}
}void update(int &rt, int pre, int l, int r, int L, int R, int v) {rt = ++num;ls[rt] = ls[pre], rs[rt] = rs[pre], sum[rt] = sum[pre], lazy[rt] = lazy[pre];sum[rt] = (sum[rt] + 1ll * (min(r, R) - max(l, L) + 1) * v) % mod;if (l >= L && r <= R) {lazy[rt] = (lazy[rt] + v) % mod;return ;}int mid = l + r >> 1;if (L <= mid) {update(ls[rt], ls[pre], l, mid, L, R, v);}if (R > mid) {update(rs[rt], rs[pre], mid + 1, r, L, R, v);}
}int query(int rt1, int rt2, int l, int r, int L, int R) {if (l >= L && r <= R) {return (sum[rt2] - sum[rt1] + mod) % mod;}int mid = l + r >> 1, ans = (1ll * (min(r, R) - max(l, L) + 1) * (lazy[rt2] - lazy[rt1]) % mod + mod) % mod;if (L <= mid) {ans = (ans + query(ls[rt1], ls[rt2], l, mid, L, R)) % mod;}if (R > mid) {ans = (ans + query(rs[rt1], rs[rt2], mid + 1, r, L, R)) % mod;}return ans;
}int query(int l, int r, int rt) {int ans = 0;while (rt) {ans = (ans + query(root[l], root[r], 1, n, id[top[rt]], id[rt])) % mod;rt = fa[top[rt]];}return ans;
}int main() {// freopen("in.txt", "r", stdin);// freopen("out.txt", "w", stdout);scanf("%d %d", &n, &m);for (int i = 2, fa; i <= n; i++) {scanf("%d", &fa);add(fa + 1, i);}dfs1(1, 0);dfs2(1, 1);for (int i = 1; i <= n; i++) {int rt = i;root[i] = root[i - 1];while (rt) {update(root[i], root[i], 1, n, id[top[rt]], id[rt], 1);rt = fa[top[rt]];}}for (int i = 1, l, r, rt; i <= m; i++) {scanf("%d %d %d", &l, &r, &rt);l++, r++, rt++;printf("%d\n", query(l - 1, r, rt));}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/313802.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

.NET框架之“小马过河”

.NET框架之“小马过河”有许多流行的 .NET框架&#xff0c;大家都觉得挺“重”&#xff0c;认为很麻烦&#xff0c;重量级&#xff0c;不如其它“轻量级”框架&#xff0c;从而不愿意使用。面对形形色色的框架发愁&#xff0c;笔者也曾发愁。但我发现只要敢于尝试&#xff0c;这…

深入理解 JVM Class文件格式(三)

** JVM常量池中各数据项类型详解 ** 关于常量池的大概内容&#xff0c; 已经在 深入理解 JVM Class文件格式&#xff08;一&#xff09; 中讲解过了&#xff0c; 这篇文章中还介绍了常量池中的11种数据类型。 本文的任务是详细讲解这11种数据类型&#xff0c; 深度剖析源文件…

#6073. 「2017 山东一轮集训 Day5」距离(树链剖分 + 永久标记主席树)

#6073. 「2017 山东一轮集训 Day5」距离 给定一颗有nnn个节点带边权的树&#xff0c;以及一个排列ppp&#xff0c;path(u,v)path(u, v)path(u,v)为u,vu, vu,v路径上的点集&#xff0c;dist(u,v)dist(u, v)dist(u,v)为u,vu, vu,v之间的最短路的长度。 有mmm次询问&#xff0c;…

ML.NET 示例:搜索引擎结果排名

ML.NET 示例中文版&#xff1a;https://github.com/feiyun0112/machinelearning-samples.zh-cn/edit/master/samples/csharp/getting-started/Ranking_Web英文原版请访问&#xff1a;https://github.com/dotnet/machinelearning-samples/tree/master/samples/csharp/getting-st…

深入理解 JVM Class文件格式(四)

&#xff08;3&#xff09;CONSTANT_Integer_info 一个常量池中的CONSTANT_Integer_info数据项, 可以看做是CONSTANT_Integer类型的一个实例。 它存储的是源文件中出现的int型数据的值。 同样&#xff0c; 作为常量池中的一种数据类型&#xff0c; 它的第一个字节也是一个tag值…

.Net Core中使用Quartz.Net Vue开即用的UI管理

Quartz.NETQuartz.Net 定制UI维护了常用作业添加、删除、修改、停止、启动功能&#xff0c;直接使用cron表达式设置作业执行间隔&#xff0c;有完整的日志记录。Quartz.NET是一个功能齐全的开源作业调度系统&#xff0c;可用于从最小的应用程序到大型企业系统。Quartz.NET是一个…

Ancient Distance(妙啊!!!) [2020牛客暑期多校训练营(第四场)]

Ancient Distance 给定一颗根为111有nnn个节点的树&#xff0c;每次可以选定树上kkk节点当作特殊节点&#xff0c; 定义dis(u)dis(u)dis(u)为&#xff0c;从u−>1u->1u−>1遇上的第一个特殊点的距离&#xff0c;如果遇不上特殊点则dis(u)dis(u)dis(u)无穷大。 有nn…

深入理解 JVM Class文件格式(五)

&#xff08;8&#xff09; CONSTANT_Class_info 常量池中的一个CONSTANT_Class_info&#xff0c; 可以看做是CONSTANT_Class数据类型的一个实例。 他是对类或者接口的符号引用。 它描述的可以是当前类型的信息&#xff0c; 也可以描述对当前类的引用&#xff0c; 还可以描述对…

混沌工程详细介绍——Netflix持续交付实践探寻

内容来源&#xff1a;DevOps案例深度研究 – Netflix的文化与工程实践战队&#xff08;本文只展示部分案例PPT及研究成果&#xff0c;更多细节请关注案例分享活动&#xff0c;及本公众号&#xff09;。本案例内容贡献者&#xff1a;高金梅&#xff0c;李晓莉&#xff0c;潘雄鹰…

P4175 [CTSC2008]网络管理(整体二分)

P4175 [CTSC2008]网络管理 给定一棵有nnn个节点的树&#xff0c;点有点权&#xff0c;有两种操作&#xff1a;① 修改某个点的点权&#xff0c;② 查询两点路径间的点权第kkk大。 给定u,vu, vu,v&#xff0c;选定111号节点为根节点&#xff0c;设inf(x)inf(x)inf(x)表示从根节…

深入理解 JVM Class文件格式(六)

经过前几篇文章&#xff0c; 终于将常量池介绍完了&#xff0c; 之所以花这么大的功夫介绍常量池&#xff0c; 是因为对于理解class文件格式&#xff0c;常量池是必须要了解的&#xff0c; 因为class文件中其他地方&#xff0c;大量引用了常量池中的数据项。 对于还不了解常量池…

远程开发初探 - VS Code Remote Development

如果你是学生&#xff0c;你还在你的 windows 电脑上为各种环境配置头疼的时候&#xff0c;你应该了解一下 Remote Development。如果你喜欢 linux 的开发环境和舒适的 shell&#xff0c;但却不舍得抛弃 windows/macos 图形界面给你带来的用户体验和一些软件的兼容(QQ, 微信), …

深入理解 JVM Class文件格式(七)

本专栏列前面的一系列博客&#xff0c; 对Class文件中的一部分数据项进行了介绍。 本文将会继续介绍class文件中未讲解的信息。 先回顾一下上面一篇文章。 在上一篇博客中&#xff0c; 我们介绍了&#xff1a; this_class 对当前类的描述 super_class 对当前类的超类的描述 in…

P3250 [HNOI2016]网络(整体二分)

P3250 [HNOI2016]网络 给定一棵树&#xff0c;有三种操作&#xff1a; 给定u,v,wu, v, wu,v,w&#xff0c;表示u,vu, vu,v路径上有一个重要度为www的请求&#xff0c;给定ttt&#xff0c;第ttt个发生的请求结束&#xff0c;给定一个xxx&#xff0c;假设xxx发生故障&#xff0…

微信小程序集成腾讯云 IM SDK

1、背景因业务功能需求需要接入IM&#xff08;即时聊天&#xff09;功能&#xff0c;一开始想到的是使用 WebSocket 来实现这个功能&#xff0c;然天意捉弄&#xff08;哈哈&#xff09;服务器版本太低不支持 wx 协议&#xff08;也就不支持 WebSocket了&#xff09;不得不寻找…

深入理解 JVM Class文件格式(八)

在本专栏的第一篇文章 深入理解Java虚拟机到底是什么 中&#xff0c; 我们主要讲解了什么是虚拟机&#xff0c; 这篇博客是对JVM的一个概述。 在随后的几篇文章中&#xff0c;一直在讲解class文件格式。 在今天这篇博客中&#xff0c; 将会继续讲解class文件中的其他信息。 在本…

Function!(计蒜客 - 42386)

Function! fa(x)ax(a>0,a≠1)f_a(x) a ^ x(a > 0, \ a \neq 1)fa​(x)ax(a>0, a​1)&#xff0c;我们要求∑a2n(a∑ban⌊fa−1(b)⌋⌈fb−1(a)⌉)\sum\limits_{a 2} ^{n} \left(a \sum\limits_{b a} ^{n} \lfloor f_a ^{-1}(b) \rfloor \lceil f_b ^{-1}(a) \rce…

深入理解 JVM Class文件格式(九)

经过前八篇关于class文件的博客&#xff0c; 关于class文件格式的内容也基本上讲完了。 本文是关于class文件格式的最后一篇。 在这篇博客中&#xff0c; 将会讲解关于方法的几个属性。 理解这篇博客的内容&#xff0c; 对于理解JVM执行引擎起着重要作用。 关于虚拟机执行引擎有…

MongoDB入门及 c# .netcore客户端MongoDB.Driver2.9.1使用

MongoDB 是一个基于分布式文件存储的数据库。由 C 语言编写。旨在为 WEB 应用提供可扩展的高性能数据存储解决方案。什么场景下使用MongoDBMongoDB虽然是NoSQL(非关系型的数据库)&#xff0c;但是实际使用的时候可以当做关系型数据库来用&#xff0c;mysql等数据库中单表数据量…

#4604. The kth maximum number(整体二分 + 树套树)

#4604. The kth maximum number 给定一个大小不超过51055 \times 10 ^ 55105的矩形区域&#xff0c;有一些点有点权。 每次询问给定x1,y1,x2,y2,kx_1, y_1, x_2, y_2, kx1​,y1​,x2​,y2​,k问以x1,y1x_1, y_1x1​,y1​为右下角&#xff0c;x2,y2x_2, y_2x2​,y2​为左上角的…