P4175 [CTSC2008]网络管理(树状数组套主席树)

P4175 [CTSC2008]网络管理

每次询问两点间的第kkk大,考虑建立nnn棵主席树,每棵主席树记录的是从根节点到当前节点的状态,

由于主席树是一种类似前缀和的数据结构,可以使用差分的方式得到(u,v)(u, v)(u,v)点对之间的主席树为u+v−lca(u,v)−fa(lca(u,v))u + v - lca(u, v) - fa(lca(u, v))u+vlca(u,v)fa(lca(u,v))

即这一段路径上的信息可用这四棵主席树加减得到,但是这样只能处理静态问题,这里还增加了修改操作。

利用dfsdfsdfs序在树状数组上建立主席树,当这个点进栈的时候,把这个点的信息在树状数组上的主席树上进行+1+1+1,当这个点出栈的时候将这个点的信息在树状数组上的主席树上进行−1-11,所以对某个点查询的时候我们查询树状数组上的主席树,其所有值加起来也就是从根节点到当前节点的信息,然而我们把修改操作优化成了log⁡n×log⁡n\log n \times \log nlogn×logn,整体复杂度也就是nlog⁡n×log⁡nn \log n \times \log nnlogn×logn

#include <bits/stdc++.h>using namespace std;const int N = 8e4 + 10, maxn = 100000000;int a[N], n, m;int head[N], to[N << 1], nex[N << 1], cnt = 1;int son[N], sz[N], l[N], r[N], dep[N], top[N], fa[N], tot;int root[N], sum[N << 8], ls[N << 8], rs[N << 8], num;int A[N], B[N], cnt1, cnt2;inline int lowbit(int x) {return x & (-x);
}inline void add(int x, int y) {to[cnt] = y;nex[cnt] = head[x];head[x] = cnt++;
}void update(int &rt, int l, int r, int x, int v) {if (!rt) {rt = ++num;}sum[rt] += v;if (l == r) {return ;}int mid = l + r >> 1;if (x <= mid) {update(ls[rt], l, mid, x, v);}else {update(rs[rt], mid + 1, r, x, v);}
}void dfs1(int rt, int f) {fa[rt] = f, l[rt] = ++tot, sz[rt] = 1, dep[rt] = dep[f] + 1;for (int i = l[rt]; i <= n; i += lowbit(i)) {update(root[i], 0, maxn, a[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];}}r[rt] = tot;for (int i = r[rt] + 1; i <= n; i += lowbit(i)) {update(root[i], 0, maxn, a[rt], -1);}
}void dfs2(int rt, int tp) {top[rt] = tp;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]);}
}int lca(int x, int y) {while (top[x] != top[y]) {if (dep[top[x]] < dep[top[y]]) {swap(x, y);}x = fa[top[x]];}return dep[x] < dep[y] ? x : y;
}int query(int l, int r, int k) {if (l == r) {return l;}int cur = 0, mid = l + r >> 1;for (int i = 1; i <= cnt1; i++) {cur += sum[ls[A[i]]];}for (int i = 1; i <= cnt2; i++) {cur -= sum[ls[B[i]]];}if (cur >= k) {for (int i = 1; i <= cnt1; i++) {A[i] = ls[A[i]];}for (int i = 1; i <= cnt2; i++) {B[i] = ls[B[i]];}return query(l, mid, k);}else {for (int i = 1; i <= cnt1; i++) {A[i] = rs[A[i]];}for (int i = 1; i <= cnt2; i++) {B[i] = rs[B[i]];}return query(mid + 1, r, k - cur);}
}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 = 1; i <= n; i++) {scanf("%d", &a[i]);}for (int i = 1; i < n; i++) {int x, y;scanf("%d %d", &x, &y);add(x, y);add(y, x);}dfs1(1, 0);dfs2(1, 1);for (int i = 1; i <= m; i++) {int k, u, v;scanf("%d %d %d", &k, &u, &v);if (!k) {for (int i = l[u]; i <= n; i += lowbit(i)) {update(root[i], 0, maxn, a[u], -1);}for (int i = r[u] + 1; i <= n; i += lowbit(i)) {update(root[i], 0, maxn, a[u], 1);}a[u] = v;for (int i = l[u]; i <= n; i += lowbit(i)) {update(root[i], 0, maxn, a[u], 1);}for (int i = r[u] + 1; i <= n; i += lowbit(i)) {update(root[i], 0, maxn, a[u], -1);}}else {int f = lca(u, v), ff = fa[f], total = dep[u] + dep[v] - dep[f] - dep[ff];if (total < k) {puts("invalid request!");}else {cnt1 = cnt2 = 0;for (int i = l[u]; i; i -= lowbit(i)) {A[++cnt1] = root[i];}for (int i = l[v]; i; i -= lowbit(i)) {A[++cnt1] = root[i];}for (int i = l[f]; i; i -= lowbit(i)) {B[++cnt2] = root[i];}for (int i = l[ff]; i; i -= lowbit(i)) {B[++cnt2] = root[i];}k = total - k + 1;printf("%d\n", query(0, maxn, k));}}}return 0;
}

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

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

相关文章

Exceptionless使用介绍

一.Exceptionless介绍Exceptionless专注于.net平台提供实时错误和日志报告。主要包括&#xff1a;错误通知、智能分组异常、详细错误报告堆栈跟踪、支持离线、UI查看重要错误和确定优先级、仪表板上的统计信息和趋势、对异常标记为已修复&#xff0c;监视回归、将事件标记为关键…

BIO、NIO 和 AIO

BIO、NIO 和 AIO的区别 https://www.cnblogs.com/barrywxx/p/8430790.html 三种 IO 的用法与原理 netty http://ifeve.com/netty5-user-guide/

#279. [SYZOI Round1] 滑稽♂树(树状数组套主席树)

#279. [SYZOI Round1] 滑稽♂树 子树上的问题&#xff0c;考虑dfsdfsdfs序&#xff0c;第kkk大&#xff0c;可以用主席树嘛&#xff0c;支持修改&#xff0c;那就树状数组上套主席树&#xff0c;参考P4175 [CTSC2008]网络管理&#xff08;树状数组套主席树&#xff09; #incl…

Net Core DocXCore 实现word模板导出

实际工作中&#xff0c;往往有这样的需求&#xff0c;需要导出word&#xff0c;还有各种各样的样式&#xff0c;于是有了word模板导出。实现以下几个需求&#xff1a;1、表单导出2、表格导出3、表单表格混合导出4、实际用例测试解决方案&#xff1a;实现是基于NET Core 2.1 &am…

主席树有关的一些题目(持续更新)

主席树 模板 P3919 【模板】可持久化线段树 1&#xff08;可持久化数组&#xff09; #include <bits/stdc.h>using namespace std;const int N 1e6 10;int root[N], n, m;int ls[N * 25], rs[N * 25], value[N * 25], tot;void build(int &rt, int l, int r) {r…

密钥分离,.Net程序猿不再背锅

引言互联网每隔一段时间就会爆出 【某程序猿在代码托管平台上传了公司机密配置信息&#xff0c;导致公司核心数据被黑客获取或修改】&#xff0c; 一茬又一茬背锅侠层出不穷。拯救自我软件工程理论早以加粗字体给出 经典原则&#xff1a;Never store production passwords or o…

Java ASM与Javassit

操作java字节码的工具有两个比较流行&#xff0c;一个是ASM&#xff0c;一个是Javassit 。 ASM &#xff1a;直接操作字节码指令&#xff0c;执行效率高&#xff0c;要是使用者掌握Java类字节码文件格式及指令&#xff0c;对使用者的要求比较高。 Javassit 提供了更高级的API…

Chrome让人失望,是时候转到Firefox或Edge?

我们曾报道过谷歌正在移除 Chrome 的关闭其他选项卡&#xff0c;这引起很多用户的讨论&#xff0c;很多用户表现出对 Chrome 浏览器的失望。而在国外&#xff0c;也有很多人认为谷歌这项修改使 Chrome 变得越来越臃肿。还有前段时间&#xff0c;谷歌在 Chrome 76 中隐藏了 “ht…

ICPC 徐州 H Yuuki and a problem (树状数组套主席树)

Yuuki and a problem 先不管第一问的修改操作&#xff0c;考虑如何达到第二问的查询操作&#xff0c; 题目要我们给出一个区间[l,r][l, r][l,r]中&#xff0c;不能通过权值得到的最小的数字是什么&#xff0c; 假设我们已经可以得到[1,x][1, x][1,x]之间的数了&#xff0c;且…

Java序列化、反序列化

目录 1&#xff1a;什么是序列化、反序列化&#xff1f; 2&#xff1a;序列化的用途&#xff1f; 3&#xff1a;序列化的n种方式 1&#xff1a;什么是序列化、反序列化&#xff1f; 把对象转换为字节序列的过程称为对象的序列化把字节序列转换为对象的过程中称为对象的反序列…

我认真写下9段如翔一般的代码,只为等你来品鉴

溪源 | 长沙.NET技术社区开篇我们总是很容易就能写出满足某个特定功能的代码&#xff0c;却很难写出优雅代码。又最欣赏那些优雅的代码&#xff0c;因为优雅代码更能体现一个开发者的积累。就像写一篇散文&#xff0c;有的就像初学者不得其门而入&#xff0c;遣词造句都非常困难…

P2596 [ZJOI2006]书架(fhq treap)

P2596 [ZJOI2006]书架 我们用fhq treap来完成这一题 对于一个新插入的节点我们取权值为其索引值&#xff0c;其所记录的valuevaluevalue是其当前索引所在位置。 操作一&#xff1a;把索引为valuevaluevalue的点放到平衡树前面&#xff0c;分别别得到三颗子树x,y,zx, y, zx,y…

「标签管理」用数据管理思维去管理你的日常电子化资料、文件、笔记等

最近一时兴起&#xff0c;研究了一些文件管理的方法论和笔记管理类的知识&#xff0c;自己想到一些不错的方法及落地方案&#xff0c;可能对一部分朋友来说&#xff0c;这些方法和工具的落地会有一些共鸣&#xff0c;故简单给大家做一下分享。跨界应用&#xff1a;数据表结构应…

Java之JMS

一&#xff1a;JMS简介   JMS即Java消息服务&#xff08;Java Message Service&#xff09;&#xff0c;是一个Java平台定义的关于面向消息中间件的API&#xff0c;用于在两个应用程序之间&#xff0c;或分布式系统中发送消息&#xff0c;进行异步通信。JMS是一个与具体平台无…

程序员过关斩将--更加优雅的Token认证方式JWT

点击上方“蓝字”带你去看小星星菜菜&#xff0c;上次你讲的cookie和session认证方式&#xff0c;我这次面试果然遇到了结果怎么样&#xff1f;结果面试官问我还有没有更好的方式&#xff1f;看来你又挂了别说了&#xff0c;伤心呀。到底还有没有更好的方式呢&#xff1f;你猜&…

Java之JMX

JMX&#xff0c;即Java Management Extensions&#xff0c;Java管理扩展。是一个为应用程序、设备、系统等植入管理功能的框架。   JMX提供了一种简单、基础的方法&#xff0c;用来管理应用、设置、服务等资源。由于JMX是动态的&#xff0c;你可以使用JMX技术来监控和管理处于…

平衡树(模板 and 题目)记录

平衡树 替罪羊树 #include <bits/stdc.h>using namespace std;const double alpha 0.725; const int N 2e6 10;struct Spgtree {int ls[N], rs[N], val[N], num[N], fac[N], sz[N], sum[N], cnt, root;int top, stk[N];void update(int rt) {fac[rt] fac[ls[rt]] …

[NewLife.XCode]分表分库(百亿级大数据存储)

NewLife.XCode是一个有15年历史的开源数据中间件&#xff0c;支持netcore/net45/net40&#xff0c;由新生命团队(2002~2019)开发完成并维护至今&#xff0c;以下简称XCode。整个系列教程会大量结合示例代码和运行日志来进行深入分析&#xff0c;蕴含多年开发经验于其中&#xf…

Java -- 泛型

1、什么是泛型&#xff1f; 泛型&#xff08;Generics&#xff09;是把类型参数化&#xff0c;运用于类、接口、方法中&#xff0c;在调用时传入具体的类型。 泛型就是参数化类型 适用于多种数据类型执行相同的代码泛型的类型在使用时指定泛型归根到底就是“模板” 优点&…

HDU 6703 array(主席树 + set)

array 给一个全排列&#xff0c;接下来有两种操作&#xff1a; 一、把pospospos位置上的值10,000,00010,000,00010,000,000。 二、查询[1,r][1, r][1,r]区间&#xff0c;没有出现的且≥k\geq k≥k的最小值是多少。 考虑用主席树 set 求解&#xff0c; 先建立一颗主席树&a…