BZOJ 1146: [CTSC2008]网络管理Network( 树链剖分 + 树状数组套主席树 )

树链剖分完就成了一道主席树裸题了, 每次树链剖分找出相应区间然后用BIT+(可持久化)权值线段树就可以完成计数. 但是空间问题很严重....在修改时不必要的就不要新建, 直接修改原来的..详见代码. 时间复杂度O(N*log^3(N))

----------------------------------------------------------------------------

#include<cstdio>
#include<cstring>
#include<cctype>
#include<iostream>
#include<algorithm>
using namespace std;
#define h(v) (lower_bound(H, H + n, v) - H + 1)
const int maxn = 80009;
inline int read() {
char c = getchar();
for(; !isdigit(c); c = getchar());
int ret = 0;
for(; isdigit(c); c = getchar())
ret = ret * 10 + c - '0';
return ret;
}
int N, Q, w[maxn], seq[maxn], H[maxn << 1], n;
int top[maxn], fa[maxn], sz[maxn], ch[maxn], dep[maxn], Top;
int Id[maxn], Idn;
struct edge {
int to;
edge* next;
} E[maxn << 1], *Pt = E, *head[maxn];
inline void AddEdge(int u, int v) {
Pt->to = v;
Pt->next = head[u];
head[u] = Pt++;
}
struct O {
int t, a, b;
inline void Read() {
t = read(), a = read() - 1, b = read();
if(t) b--;
}
} o[maxn];
void Init() {
N = read(), Q = read();
for(n = 0; n < N; n++)
H[n] = w[n] = read();
for(int i = 1; i < N; i++) {
int u = read() - 1, v = read() - 1;
AddEdge(u, v);
AddEdge(v, u);
}
for(int i = 0; i < Q; i++) {
o[i].Read();
if(!o[i].t)
H[n++] = o[i].b;
}
sort(H, H + n);
n = unique(H, H + n) - H;
}
struct Node {
Node *lc, *rc;
int v;
Node() : v(0) {
}
} pool[5000000], *pt, *Null, *Root[maxn], *V[maxn];
void Init_sgt() {
pt = pool;
Null = pt++;
Null->lc = Null->rc = Null;
Null->v = 0;
}
int Pos, Val;
Node* Modify(Node* t, int l, int r) {
if(t->v + Val == 0)
return Null;
if(t->v == Val) { 
if(l == r) return t;
int m = (l + r) >> 1;
if(Pos <= m) {
t->lc = Modify(t->lc, l, m);
t->rc = Null;
} else {
t->lc = Null;
t->rc = Modify(t->rc, m + 1, r);
}
return t;
} else {
Node* o = pt++;
o->v = t->v + Val;
int m = (l + r) >> 1;
if(Pos <= m) {
o->lc = Modify(t->lc, l, m);
o->rc = t->rc;
} else {
o->lc = t->lc;
o->rc = Modify(t->rc, m + 1, r);
}
return o;
}
}
Node* Add(Node* t, int l, int r) {
Node* o = pt++;
o->v = t->v + Val;
if(l != r) {
int m = (l + r) >> 1;
if(Pos <= m) {
o->lc = Add(t->lc, l, m);
o->rc = t->rc;
} else {
o->lc = t->lc;
o->rc = Add(t->rc, m + 1, r);
}
}
return o;
}
void DFS(int x) {
sz[x] = 1;
ch[x] = -1;
for(edge* e = head[x]; e; e = e->next) if(e->to != fa[x]) {
dep[e->to] = dep[x] + 1;
fa[e->to] = x;
DFS(e->to);
sz[x] += sz[e->to];
if(ch[x] == -1 || sz[ch[x]] < sz[e->to])
ch[x] = e->to;
}
}
void dfs(int x) {
seq[Id[x] = ++Idn] = x;
top[x] = Top;
if(~ch[x]) dfs(ch[x]);
for(edge* e = head[x]; e; e = e->next)
if(e->to != fa[x] && e->to != ch[x]) dfs(Top = e->to);
}
void Init_slpf() {
fa[0] = -1;
dep[0] = 0;
DFS(0);
Idn = 0;
dfs(0);
}
void Build() {
Init_sgt();
Root[0] = V[0] = Null;
Val = 1;
for(int i = 1; i <= N; i++) {
Pos = h(w[seq[i]]);
Root[i] = Add(Root[i - 1], 1, n);
V[i] = Null;
}
}
Node *L[500], *R[500];
void Work() {
int Ln, Rn;
for(int i = 0; i < Q; i++) {
if(o[i].t) {
int x = o[i].a, y = o[i].b, size = 0;
Ln = Rn = 0;
for(; top[x] != top[y]; x = fa[top[x]]) {
if(dep[top[x]] < dep[top[y]]) swap(x, y);
size += dep[x] - dep[top[x]] + 1;
L[Ln++] = Root[Id[top[x]] - 1];
R[Rn++] = Root[Id[x]];
for(int l = Id[top[x]] - 1; l; l -= l & -l)
L[Ln++] = V[l];
for(int r = Id[x]; r; r -= r & -r)
R[Rn++] = V[r];
}
if(dep[x] > dep[y]) swap(x, y);
size += dep[y] - dep[x] + 1;
L[Ln++] = Root[Id[x] - 1];
R[Rn++] = Root[Id[y]];
for(int l = Id[x] - 1; l; l -= l & -l)
L[Ln++] = V[l];
for(int r = Id[y]; r; r -= r & -r)
R[Rn++] = V[r];
if(size < o[i].t) {
puts("invalid request!");
continue;
} else 
o[i].t = size - o[i].t + 1;
int l = 1, r = n;
while(l < r) {
int m = (l + r) >> 1, cnt = 0;
for(int j = 0; j < Ln; j++)
cnt -= L[j]->lc->v;
for(int j = 0; j < Rn; j++)
cnt += R[j]->lc->v;
if(cnt >= o[i].t) {
for(int j = 0; j < Ln; j++)
L[j] = L[j]->lc;
for(int j = 0; j < Rn; j++)
R[j] = R[j]->lc;
r = m;
} else {
o[i].t -= cnt;
for(int j = 0; j < Ln; j++)
L[j] = L[j]->rc;
for(int j = 0; j < Rn; j++)
R[j] = R[j]->rc;
l = m + 1;
}
}
printf("%d\n", H[l - 1]);
} else {
Pos = h(w[seq[Id[o[i].a]]]), Val = -1;
for(int x = Id[o[i].a]; x <= N; x += x & -x)
V[x] = Add(V[x], 1, n);
Pos = h(w[seq[Id[o[i].a]]] = o[i].b), Val = 1;
for(int x = Id[o[i].a]; x <= N; x += x & -x)
V[x] = Add(V[x], 1, n);
}
}
}
int main() {
Init();
Init_slpf();
Build();
Work();
return 0;
}

 

----------------------------------------------------------------------------

1146: [CTSC2008]网络管理Network

Time Limit: 50 Sec  Memory Limit: 162 MB
Submit: 2943  Solved: 859
[Submit][Status][Discuss]

Description

M公司是一个非常庞大的跨国公司，在许多国家都设有它的下属分支机构或部门。为了让分布在世界各地的N个部门之间协同工作，公司搭建了一个连接整个公司的通信网络。该网络的结构由N个路由器和N-1条高速光缆组成。每个部门都有一个专属的路由器，部门局域网内的所有机器都联向这个路由器，然后再通过这个通信子网与其他部门进行通信联络。该网络结构保证网络中的任意两个路由器之间都存在一条直接或间接路径以进行通信。 高速光缆的数据传输速度非常快，以至于利用光缆传输的延迟时间可以忽略。但是由于路由器老化，在这些路由器上进行数据交换会带来很大的延迟。而两个路由器之间的通信延迟时间则与这两个路由器通信路径上所有路由器中最大的交换延迟时间有关。作为M公司网络部门的一名实习员工，现在要求你编写一个简单的程序来监视公司的网络状况。该程序能够随时更新网络状况的变化信息（路由器数据交换延迟时间的变化），并且根据询问给出两个路由器通信路径上延迟第k大的路由器的延迟时间。【任务】 你的程序从输入文件中读入N个路由器和N-1条光缆的连接信息，每个路由器初始的数据交换延迟时间Ti，以及Q条询问（或状态改变）的信息。并依次处理这Q条询问信息，它们可能是： 1. 由于更新了设备，或者设备出现新的故障，使得某个路由器的数据交换延迟时间发生了变化。 2. 查询某两个路由器a和b之间的路径上延迟第k大的路由器的延迟时间。

Input

第一行为两个整数N和Q，分别表示路由器总数和询问的总数。第二行有N个整数，第i个数表示编号为i的路由器初始的数据延迟时间Ti。紧接着N-1行，每行包含两个整数x和y。表示有一条光缆连接路由器x和路由器y。紧接着是Q行，每行三个整数k、a、b。如果k=0，则表示路由器a的状态发生了变化，它的数据交换延迟时间由Ta变为b。如果k>0，则表示询问a到b的路径上所经过的所有路由器（包括a和b）中延迟第k大的路由器的延迟时间。注意a可以等于b，此时路径上只有一个路由器。

Output

对于每一个第二种询问（k>0），输出一行。包含一个整数为相应的延迟时间。如果路径上的路由器不足k个，则输出信息“invalid request!”（全部小写不包含引号，两个单词之间有一个空格）。

Sample Input

5 5
5 1 2 3 4
3 1
2 1
4 3
5 3
2 4 5
0 1 2
2 2 3
2 1 4
3 3 5

Sample Output

3
2
2
invalid request!

HINT

10% 测试数据满足N<=8000,Q<=3000,

40% 测试数据满足所有询问中1<=K<=5 。即路由器的延迟时间不会发生变化。

100% 测试数据满足N,Q<=80000,任意一个路由器在任何时刻都满足延迟时间小于10^8。对于所有询问满足0<=K<=N 。

Source

 

转载于:https://www.cnblogs.com/JSZX11556/p/5128213.html

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

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

相关文章

struts2遍历select

struts2遍历select

2019独角兽企业重金招聘Python工程师标准>>> <s:bean name"org.apache.struts2.util.Counter" id"counter"> <s:param name"first" value"1" /> <s:param name"last" value"pageCount"…
阅读更多...
Oracle 彻底 kill session

Oracle 彻底 kill session

为什么80%的码农都做不了架构师&#xff1f;>>> --*************************** -- Oracle 彻底 kill session --*************************** kill session 是DBA经常碰到的事情之一。如果kill 掉了不该kill 的session&#xff0c;则具有破坏性&#xff0c;因此尽…
阅读更多...
WebService学习总结——调用第三方提供的webService服务

WebService学习总结——调用第三方提供的webService服务

互联网上面有很多的免费webService服务&#xff0c;我们可以调用这些免费的WebService服务&#xff0c;将一些其他网站的内容信息集成到我们的Web应用中显示&#xff0c;下面就以获取天气预报数据。气象中心的管理系统将收集的天气信息并将数据暴露出来(通过WebService Server)…
阅读更多...
【017】◀▶ C#学习(九) - ADO.NET

【017】◀▶ C#学习(九) - ADO.NET

《C#入门经典&#xff08;中文第四版&#xff09;》在程序中访问数据库学习笔记 --------------------------------------------------------------------------------------------------------- ●● 目录&#xff1a; A0 ………… System.Data.SqlClient 命名空间 A1 ……………
阅读更多...
【转】主流PHP框架间的比较（Zend Framework，CakePHP，CodeIgnit...

【转】主流PHP框架间的比较（Zend Framework，CakePHP，CodeIgnit...

2019独角兽企业重金招聘Python工程师标准>>> 【转】主流PHP框架间的比较&#xff08;Zend Framework&#xff0c;CakePHP&#xff0c;CodeIgniter&#xff0c;Symfony&#xff0c;ThinkPHP&#xff0c;FleaPHP&#xff09; 2011年08月14日 星期日 12:51 转载自 leol…
阅读更多...
如何利用业务时间提升自我

如何利用业务时间提升自我

转载于:https://www.cnblogs.com/helloIT/articles/5140273.html
阅读更多...
ping得通外网，上得了QQ，游戏，却打不开网页。

ping得通外网，上得了QQ，游戏，却打不开网页。

操作系统&#xff1a;win7_64bit 下午帮我小艾优化电脑&#xff0c;不知道把什么插件给卸载了。 她回来以上就上不了网&#xff0c;我首先ping网关&#xff0c;DNS&#xff0c;外网&#xff0c;都通。 登QQ正常。 我的思路&#xff1a; 先把Socket套字结相关的卸载掉就是&#…
阅读更多...
WebGIS中一种根据网格索引判断点面关系的方法

WebGIS中一种根据网格索引判断点面关系的方法

文章版权由作者李晓晖和博客园共有&#xff0c;若转载请于明显处标明出处&#xff1a;http://www.cnblogs.com/naaoveGIS/。 1.背景 判断点面关系的算法有很多&#xff0c;在我之前的博文中有一篇专门对其进行了描述&#xff1a;判断点是否落在面中的Oracle存储过程描述。其中提…
阅读更多...
为Ubuntu Server 安装图形桌面环境

为Ubuntu Server 安装图形桌面环境

From: http://www.oschina.net/question/12_10835 第一步&#xff1a;安装桌面环境 首先需要确保您已经在/etc/apt/sources.list 文件中启用了Universe和Multiverse软件库&#xff0c;接下来执行如下命令&#xff1a; sudo apt-get update sudo apt-get install ubuntu-desk…
阅读更多...
iSensor APP 之 摄像头调试 OV3640 OV2640 MT9d112

iSensor APP 之 摄像头调试 OV3640 OV2640 MT9d112

iSensor app 非常适合调试各种摄像头&#xff0c;已测试通过的sensor有&#xff1a; l OV7670、OV7725、OV9650、OV9655、OV9653、OV5642、OV5640 l MT9T001、MT9M001、MT9P031 今天又测试了上述几款新的型号&#xff0c;之前并未用过。 OV3640 OV2640 MT9d112 先看OV2640&am…
阅读更多...
查询优化器内核剖析第四篇：从一个实例看执行计划

查询优化器内核剖析第四篇：从一个实例看执行计划

查询优化器内核剖析第四篇&#xff1a;从一个实例看执行计划系列文章索引&#xff1a; 查询优化器内核剖析第一篇 查询优化器内核剖析第二篇&#xff1a;产生候选执行计划&执行计划成本估算 查询优化器内核剖析第三篇&#xff1a;查询的执行与计划的缓存 &…
阅读更多...
Makefile中自定义函数的调用

Makefile中自定义函数的调用

From: http://www.cnblogs.com/MyEyes/archive/2012/01/12/2320589.html 自己学习脚印&#xff0c;不喜勿喷&#xff0c;谢谢 ~ Makefile中函数定义&#xff1a; external/genext2fs/Config.mk # $(1): src directory # $(2): output file # $(3): label (if any) # $(4): …
阅读更多...
python selenium ---键盘事件

python selenium ---键盘事件

转自&#xff1a;http://www.cnblogs.com/fnng/p/3258946.html 本节重点&#xff1a; l 键盘按键用法l 键盘组合键用法l send_keys() 输入中文运行报错问题键盘按键键用法&#xff1a; #codingutf-8 from selenium import webdriver from selenium.webdriver.common.keys impor…
阅读更多...
Eclipse中使用自己的makefile管理工程

Eclipse中使用自己的makefile管理工程

From: http://blog.csdn.net/whz_zb/article/details/7101164 我用的环境&#xff1a; &#xff08;1&#xff09;系统&#xff1a; Ubuntu 11.04 &#xff08;2&#xff09;Eclipse CDT 1. 在建立工程的时候选择makefile工程。 2.将makefile文件和你写好的工程文件&#x…
阅读更多...
Maven实战(Maven+Nexus建立私服【Linux系统】)

Maven实战(Maven+Nexus建立私服【Linux系统】)

准备工作 下载及配置Maven3&#xff1a;http://www.cnblogs.com/leefreeman/archive/2013/03/05/2944519.html 下载Nexus&#xff1a;http://nexus.sonatype.org/downloads/ 安装配置Nexus Nexus提供了两种安装方式&#xff0c;一种是内嵌Jetty的bundle&#xff0c;只要你有JRE…
阅读更多...
青少年如何使用 Python 开始游戏开发

青少年如何使用 Python 开始游戏开发

From: http://www.oschina.net/translate/beginning-game-programming-for-teens-with-python 这是一篇由教程团队成员Julian Meyer发表的文章,一个12岁的python开发人员。你可以在Google和Twitter上找到他。 你可曾想过如何创建视频游戏吗&#xff1f;它不像你想象的那么复杂…
阅读更多...
Vue3中的混入（mixins）

Vue3中的混入（mixins）

本文主要介绍Vue3中的混入&#xff08;mixins&#xff09;。 目录 一、在普通写法中使用混入&#xff1a;二、在setup写法中使用混入&#xff1a; 混入是Vue中一种用于在组件中共享可复用功能的特性。在Vue 3中&#xff0c;混入的使用方式有所改变。 一、在普通写法中使用混入…
阅读更多...
RTSP协议

RTSP协议

From: http://www.cnblogs.com/Jimmly/archive/2009/07/27/1531999.html RTSP是由Real network 和Netscape共同提出的如何有效地在IP网络上传输流媒体数据的应用层协议。 实时流协议&#xff08;RTSP&#xff09;建立并控制一个或几个时间同步的连续流媒体&#xff0c;如音频和…
阅读更多...
CSS3 动画 思维导图

CSS3 动画 思维导图

思维导图在新窗口打开浏览 转载于:https://www.cnblogs.com/weaming/p/5193044.html
阅读更多...
python+PyQT+Eric安装配置

python+PyQT+Eric安装配置

From: http://www.cnblogs.com/lhj588/archive/2011/10/03/2198472.html 一、大纲内容&#xff1a; 1、预备PC环境&#xff1a; 2、预备安装程序: 2、1、下载Python3.2 2、2、下载PyQt4 2、3、下载Eric5 3、安装配置步骤&#xff1a; 3、1、安装Pyhon3.2 3、2、安装PyQt4 3、3、…
阅读更多...
推荐文章
最新文章

Copyright @ 2022~2023