题目描述

给定一个含有n个数的序列a[1],a[2],a[3]……a[n],程序必须回答这样的询问:对于给定的i,j,k,在a[i],a[i+1],a[i+2]……a[j]中第k小的数是多少(1≤k≤j-i+1),并且,你可以改变一些a[i]的值,改变后,程序还能针对改变后的a继续回答上面的问题。你需要编一个这样的程序,从输入文件中读入序列a,然后读入一系列的指令,包括询问指令和修改指令。对于每一个询问指令,你必须输出正确的回答。

输入

第一行有两个正整数n(1≤n≤10000),m(1≤m≤10000)。分别表示序列的长度和指令的个数。第二行有n个数,表示a[1],a[2]……a[n],这些数都小于10^9。接下来的m行描述每条指令,每行的格式是下面两种格式中的一种。 Q i j k 或者 C i t Q i j k (i,j,k是数字,1≤i≤j≤n, 1≤k≤j-i+1)表示询问指令,询问a[i],a[i+1]……a[j]中第k小的数。C i t (1≤i≤n,0≤t≤10^9)表示把a[i]改变成为t。

输出

对于每一次询问,你都需要输出他的答案,每一个输出占单独的一行。

样例输入

5 3
3 2 1 4 7
Q 1 4 3
C 2 6
Q 2 5 3

样例输出

3
6


题解

动态主席树,用树状数组维护。

静态的主席树原理类似于前缀和,查询的时间复杂度为O(nlogn),而修改的时间复杂度应为O(n^2·logn)。

所以需要一种新的方式来实现,很容易想到树状数组。

这相当于n棵动态开点线段树建立在一个树状数组上。

查询和修改时需要按照树状数组的方式来操作,其余的和静态主席树基本一致。

还有,离散化时,需要将每个修改之后的值也加入到数组中。故需要先读入所有询问再操作。

#include <cstdio>
#include <algorithm>
using namespace std;
struct data
{
int num , rank;
}v[40001];
int a[20001] , flag[20001] , A[20001] , B[20001] , K[20001] , val[40001] , root[20001];
int L[40] , R[40] , tot , top , cnt;
int si[4400000] , lp[4400000] , rp[4400000];
char str[5];
bool cmp1(data a , data b)
{
return a.num < b.num;
}
bool cmp2(data a , data b)
{
return a.rank < b.rank;
}
void pushup(int x)
{
si[x] = si[lp[x]] + si[rp[x]];
}
void update(int x , int &y , int l , int r , int p , int a)
{
y = ++tot;
if(l == r)
{
si[y] = si[x] + a;
return;
}
int mid = (l + r) >> 1;
if(p <= mid) rp[y] = rp[x] , update(lp[x] , lp[y] , l , mid , p , a);
else lp[y] = lp[x] , update(rp[x] , rp[y] , mid + 1 , r , p , a);
pushup(y);
}
int query(int a , int b , int l , int r , int p)
{
if(l == r) return val[l];
int i , suml = 0 , sumr = 0 , mid = (l + r) >> 1;
for(i = 1 ; i <= a ; i ++ ) suml += si[lp[L[i]]];
for(i = 1 ; i <= b ; i ++ ) sumr += si[lp[R[i]]];
if(sumr - suml >= p)
{
for(i = 1 ; i <= a ; i ++ ) L[i] = lp[L[i]];
for(i = 1 ; i <= b ; i ++ ) R[i] = lp[R[i]];
return query(a , b , l , mid , p);
}
else
{
for(i = 1 ; i <= a ; i ++ ) L[i] = rp[L[i]];
for(i = 1 ; i <= b ; i ++ ) R[i] = rp[R[i]];
return query(a , b , mid + 1 , r , p - sumr + suml);
}
}
int main()
{
int n , m , i , j , x , y;
scanf("%d%d" , &n , &m);
for(i = 1 ; i <= n ; i ++ )
{
scanf("%d" , &a[i]);
v[++top].num = a[i];
v[top].rank = top;
}
for(i = 1 ; i <= m ; i ++ )
{
scanf("%s%d%d" , str , &A[i] , &B[i]);
if(str[0] == 'Q') scanf("%d" , &K[i]) , flag[i] = 1;
else v[++top].num = B[i] , v[top].rank = top;
}
sort(v + 1 , v + top + 1 , cmp1);
val[0] = -1;
for(i = 1 ; i <= top ; i ++ )
{
if(v[i].num == val[cnt]) v[i].num = cnt;
else val[++cnt] = v[i].num , v[i].num = cnt;
}
sort(v + 1 , v + top + 1 , cmp2);
for(i = 1 ; i <= n ; i ++ )
for(j = i ; j <= n ; j += j & (-j))
update(root[j] , root[j] , 1 , cnt , v[i].num , 1);
top = n;
for(i = 1 ; i <= m ; i ++ )
{
if(flag[i])
{
x = y = 0;
for(j = A[i] - 1 ; j > 0 ; j -= j & (-j))
L[++x] = root[j];
for(j = B[i] ; j > 0 ; j -= j & (-j))
R[++y] = root[j];
printf("%d\n" , query(x , y , 1 , cnt , K[i]));
}
else
{
for(j = A[i] ; j <= n ; j += j & (-j))
update(root[j] , root[j] , 1 , cnt , v[A[i]].num , -1);
top ++ ;
v[A[i]].num = v[top].num;
for(j = A[i] ; j <= n ; j += j & (-j))
update(root[j] , root[j] , 1 , cnt , v[A[i]].num , 1);
}
}
return 0;
}

【bzoj1901】Zju2112 Dynamic Rankings 离散化+主席树+树状数组的更多相关文章

  1. [BZOJ1901]Zju2112 Dynamic Rankings

    [BZOJ1901]Zju2112 Dynamic Rankings 试题描述 给定一个含有n个数的序列a[1],a[2],a[3]--a[n],程序必须回答这样的询问:对于给定的i,j,k,在a[i ...

  2. BZOJ-1901 Zju2112 Dynamic Rankings 函数式线段树 套 树状数组+离线处理

    1901: Zju2112 Dynamic Rankings Time Limit: 10 Sec Memory Limit: 128 MB Submit: 6058 Solved: 2521 [Su ...

  3. BZOJ1901 Zju2112 Dynamic Rankings 主席树

    欢迎访问~原文出处——博客园-zhouzhendong 去博客园看该题解 题目传送门 - BZOJ1901 题意概括 给你一段序列(n个数),让你支持一些操作(共m次), 有两种操作,一种是询问区间第 ...

  4. [luogu2617][bzoj1901][Zju2112]Dynamic Rankings【树套树+树状数组+主席树】

    题目网址 [传送门] 题目大意 请你设计一个数据结构,支持单点修改,区间查询排名k. 感想(以下省略脏话inf个字) 真的强力吹爆洛谷数据,一般的树套树还给我T了一般的点,加强的待修主席树还给我卡了几 ...

  5. BZOJ1901 Zju2112 Dynamic Rankings 【树状数组套主席树】

    题目 给定一个含有n个数的序列a[1],a[2],a[3]--a[n],程序必须回答这样的询问:对于给定的i,j,k,在a[i],a[i+1],a[i+2]--a[j]中第k小的数是多少(1≤k≤j- ...

  6. 【树状数组套权值线段树】bzoj1901 Zju2112 Dynamic Rankings

    谁再管这玩意叫树状数组套主席树我跟谁急 明明就是树状数组的每个结点维护一棵动态开结点的权值线段树而已 好吧,其实只有一个指针,指向该结点的权值线段树的当前结点 每次查询之前,要让指针指向根结点 不同结 ...

  7. bzoj 1901: Zju2112 Dynamic Rankings【整体二分+树状数组||主席树+树状数组】

    整体二分: 对于每一个修改操作,标记为1,并且加一个标记为-1的这个位置原来值,并且对于a数列每个点都当成修改操作 然后整体二分,扫当前操作区间lr,把在值域区间标记为1和-1的操作都在树状数组对应位 ...

  8. BZOJ1901: Zju2112 Dynamic Rankings(整体二分 树状数组)

    Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 9094  Solved: 3808[Submit][Status][Discuss] Descript ...

  9. BZOJ1901——Zju2112 Dynamic Rankings

    1.题目大意:区间第k小,有单点修改 2.分析:这个是树状数组套线段树,也是主席树....为什么主席树这么多QAQ 就是树套树的那种插入什么的,注意啊,一定要动态开内存..不然会爆.. 然后算答案有两 ...

随机推荐

  1. hive中 udf,udaf,udtf

    1.hive中基本操作: DDL,DML 2.hive中函数 User-Defined Functions : UDF(用户自定义函数,简称JDF函数)UDF: 一进一出  upper  lower ...

  2. Python之celery

    一.celery简介 Celery是一个Python开发的异步分布式任务调度模块.celery本身不提供消息服务,使用第三方服务,也就是borker来传递任务,目前支持rebbing, redis, ...

  3. DSP5509之采样定理

    1. 在实际种信号是模拟连续的,但是AD采样确实离散的数字的,根据采样定理,采样频率要是模拟信号的频率2倍以上采样到的值才没问题. 2. 打开工程 unsigned ]; main() { int i ...

  4. 阅读笔记《JavaScript语言精粹》

    阅读笔记<JavaScript语言精粹> 对象 1.检索属性 使用[]和. 2.引用传递 JavaScript的简单数据类型包括数字.字符串.布尔值.null值和undefined值.其它 ...

  5. 【vps搬家】--总结--费元星

    20150310  费元星 稍微玩VPS/服务器比较久的站长手中应该不止一台VPS,我们会有多台机器之间的相互使用.比如可能会遇到的是数据传输,我们传统的做法是先用FTP下载数据A到本地,然后再到本地 ...

  6. 错误码:2003 不能连接到 MySQL 服务器在 (10061)

    今天在ubuntu上安装了mysql服务器,在windows上用客户端软件连接mysql服务器时,出现错误: 错误码: 不能连接到 MySQL 服务器在 () 折腾来折腾去没搞好,防火墙也关了,330 ...

  7. lesson 14 A noble gangster

    lesson 14 A noble gangster there was a ++time++ 时期 times 时期/年代 in times of peace a sum of + money 一笔 ...

  8. 关于javascript的一个小问题,请问有人看出啥问题吗?

    最近学习javascript,有一个问题挺奇怪的,先贴出代码: function binarySearch(){ var arr = [0,1,2,3]; var res = actbinarySea ...

  9. [leetcode-658-Find K Closest Elements]

    Given a sorted array, two integers k and x, find the k closest elements to x in the array. The resul ...

  10. 《剑指Offer》题一~题十

    一.赋值运算符函数 题目:如下为类型CMyString的声明,请为该类型添加赋值运算符函数. class CMyString { public: CMyString(char *pData = nul ...