CRB and Queries

Time Limit: 12000/6000 MS (Java/Others)    Memory Limit: 131072/131072 K (Java/Others)
Total Submission(s): 1602    Accepted Submission(s): 409

Problem Description
There are N boys in CodeLand.
Boy i has his coding skill Ai.
CRB wants to know who has the suitable coding skill.
So you should treat the following two types of queries.
Query 1: 1 l v
The coding skill of Boy l has changed to v.
Query 2: 2 l r k
This is a report query which asks the k-th smallest value of coding skill between Boy l and Boy r(both inclusive).
 
Input
There are multiple test cases. 
The first line contains a single integer N.
Next line contains N space separated integers A1, A2, …, AN, where Ai denotes initial coding skill of Boy i.
Next line contains a single integer Q representing the number of queries.
Next Q lines contain queries which can be any of the two types.
1 ≤ N, Q ≤ 105
1 ≤ Ai, v ≤ 109
1 ≤ l ≤ r ≤ N
1 ≤ k ≤ r – l + 1

 
Output
For each query of type 2, output a single integer corresponding to the answer in a single line.
 
Sample Input
5
1 2 3 4 5
3
2 2 4 2
1 3 6
2 2 4 2
 
Sample Output
3
4
 
Author
KUT(DPRK)
 
Source
 

题目大意:给你一个数列有n个数,q次询问。询问1可以把某个位置改为值b。询问2然后问L - R之间第k大的数是多少。简单来说,就是带修改的动态区间第k大问题。

解题思路:用线段树维护Ci在去重后的所有可能出现在数列中的数当中的大小排名,用Treap维护Ci在数列中的位置关系。其实我觉得这才是内涵。然后对于查询区间L - R的第k大。那么如果在线段树左儿子代表的Treap树中在R位置之前的个数减去左儿子代表的Treap树中的(L-1)之前的个数大于k,那么就可以在左儿子中找位置。

#include<stdio.h>
#include<string.h>
#include<time.h>
#include<algorithm>
using namespace std;
#define mid (L+R)/2
#define lson rt*2,L,mid
#define rson rt*2+1,mid+1,R
const int maxn=1e6;
struct Treap{
int rk;
int coun;
int sz;
int v;
Treap *ch[2];
Treap() {
coun=0;
sz=0;
rk=-maxn; //
v=0;
}
int cmp(int x) const {
if(x==v){
return -1;
}
return x<v ? 0:1;
}
};
Treap *seg[maxn*8];
Treap *null;
int oper[maxn][4];
int a[maxn],b[maxn*3];
int discr(int l,int r,int key){
while(l<=r){
int md=(l+r)/2;
if(key<b[md]){
r=md-1;
}else if(key > b[md]){
l=md+1;
}else{
return md;
}
}
return -1;
}
void update(Treap * &o){
o->sz = o->ch[0]->sz + o->coun + o->ch[1]->sz;
}
void rotate(Treap * &o,int d){
Treap *k=o->ch[d^1]; o->ch[d^1]= k->ch[d]; k->ch[d]=o;
update(o); update(k); o=k;
}
void insert_tp(Treap * &o,int x){
if(o==null){
o = new Treap();
o->ch[0]=o->ch[1]=null;
o->v = x; o->coun = 1; o->rk = rand();
}else{
int d=o->cmp(x);
insert_tp(o->ch[d],x);
if(o->ch[d]->rk > o->rk) rotate(o,d^1);
}
update(o);
}
void insert_seg(int rt,int L,int R,int pos,int x){
insert_tp(seg[rt],x);
if(L==R)
return ;
if(pos<=mid){
insert_seg(lson,pos,x);
}else{
insert_seg(rson,pos,x);
}
}
void free_tp(Treap *&o){
if(o->ch[0]==null&&o->ch[1]==null){
free(o);
return ;
}
if(o->ch[0]!=null){
free_tp(o->ch[0]);
}
if(o->ch[1]!=null){
free_tp(o->ch[1]);
}
free(o);
}
void clean(int rt,int L,int R){
if(L==R){
free_tp(seg[rt]);
return ;
}
clean(lson);
clean(rson);
free_tp(seg[rt]);
}
void del_tp(Treap * &o,int x){
int d=o->cmp(x);
if(d==-1){
if(o->ch[0]==null&&o->ch[1]==null){ //
Treap *pt=o;
o = null;
free(pt);
}
else if(o->ch[0]==null){
Treap *pt=o;
o=o->ch[1];
free(pt);
}else if(o->ch[1]==null){
Treap *pt=o;
o= o->ch[0];
free(pt);
}else{
int d2=( o->ch[0]->rk >o->ch[1]->rk ? 1:0 );
rotate(o,d2);
del_tp(o->ch[d2],x);
}
}else{
del_tp(o->ch[d],x);
}
update(o);
}
void del_seg(int rt,int L,int R,int pos,int x){
del_tp(seg[rt],x);
if(L==R)
return ;
if(pos<=mid){
del_seg(lson,pos,x);
}else{
del_seg(rson,pos,x);
}
}
int select(Treap *o,int x){
if(o==null){
return 0;
}
if(o->v > x) return select(o->ch[0],x);
return o->ch[0]->sz + o->coun +select(o->ch[1],x);
}
int query(int rt,int L,int R,int x,int y,int k){
if(L==R) return L;
int ans=select(seg[rt*2],y)-select(seg[rt*2],x);
if(ans>=k)
return query(lson,x,y,k);
else return query(rson,x,y,k-ans);
}
void init(){
null =new Treap();
null->ch[0]=null->ch[1]=null;
null->sz = null->coun=0;
for(int i=0;i<=maxn*4-1;i++){
seg[i] = null;
}
}
int main(){
// freopen("1007.in","r",stdin);
// freopen("OUT.txt","w",stdout);
int n,Q,typ,x,y,z,nn,mm;
while(scanf("%d",&n)!=EOF){
init();
mm=0;
for(int i=0;i<n;i++){
scanf("%d",&a[i]);
b[mm++]=a[i];
}
scanf("%d",&Q);
for(int i=0;i<Q;i++){
scanf("%d",&typ);
oper[i][0]=typ;
if(typ==2){
scanf("%d%d%d",&x,&y,&z);
oper[i][1]=x;
oper[i][2]=y;
oper[i][3]=z;
}else{
scanf("%d%d",&x,&z);
oper[i][1]=x;
oper[i][2]=z;
b[mm++]=z;
}
}
sort(b,b+mm);
nn=1;
for(int i=1;i<mm;i++){
if(b[i]!=b[i-1]){
b[nn++]=b[i];
}
}
for(int i=0;i<n;i++){
a[i]= discr(0,nn-1,a[i])+1;
insert_seg(1,1,nn,a[i],i+1);
}
for(int i=0;i<Q;i++){
if(oper[i][0]==1){
int kk=0;
del_seg(1,1,nn,a[oper[i][1]-1],oper[i][1]);
kk=discr(0,nn-1,oper[i][2]);
a[oper[i][1]-1]=kk+1;
insert_seg(1,1,nn,a[oper[i][1]-1],oper[i][1]);
}else{
int tmp=query(1,1,nn,oper[i][1]-1,oper[i][2],oper[i][3])-1;
printf("%d\n",b[tmp]);
}
}
clean(1,1,nn);
}
return 0;
} /* 5
1 1 1 1 1
10
1 2 3
2 1 4 3
1 5 6
1 4 2
1 5 3
2 1 3 3
2 1 4 4
2 2 5 2
2 2 5 3
2 2 5 4 */

  

HDU 5412——CRB and Queries——————【线段树套Treap(并没有AC)】的更多相关文章

  1. HDU 5412 CRB and Queries(区间第K大 树套树 按值建树)

    题目链接:http://acm.hdu.edu.cn/showproblem.php? pid=5412 Problem Description There are N boys in CodeLan ...

  2. BZOJ3196 Tyvj1730 二逼平衡树 【树套树】 【线段树套treap】

    BZOJ3196 Tyvj1730 二逼平衡树 Description 您需要写一种数据结构(可参考题目标题),来维护一个有序数列,其中需要提供以下操作: 1.查询k在区间内的排名 2.查询区间内排名 ...

  3. BZOJ - 3196 Tyvj 1730 二逼平衡树 (线段树套treap)

    题目链接 区间线段树套treap,空间复杂度$O(nlogn)$,时间复杂度除了查询区间k大是$O(log^3n)$以外都是$O(log^2n)$的. (据说线段树套线段树.树状数组套线段树也能过?) ...

  4. 初学树套树:线段树套Treap

    前言 树套树是一个十分神奇的算法,种类也有很多:像什么树状数组套主席树.树状数组套值域线段树.\(zkw\)线段树套\(vector\)等等. 不过,像我这么弱,当然只会最经典的 线段树套\(Trea ...

  5. BZOJ2141&洛谷1975 排队 【线段树套treap】

    题目 排排坐,吃果果,生果甜嗦嗦,大家笑呵呵.你一个,我一个,大的分给你,小的留给我,吃完果果唱支歌,大家乐和和. 红星幼儿园的小朋友们排起了长长地队伍,准备吃果果.不过因为小朋友们的身高有所区别,排 ...

  6. hdu 5412 CRB and Queries

    题目连接 http://acm.hdu.edu.cn/showproblem.php?pid=5412 CRB and Queries Description There are $N$ boys i ...

  7. bzoj 3196 && luogu 3380 JoyOI 1730 二逼平衡树 (线段树套Treap)

    链接:https://www.lydsy.com/JudgeOnline/problem.php?id=3196 题面; 3196: Tyvj 1730 二逼平衡树 Time Limit: 10 Se ...

  8. 【bzoj3196】Tyvj 1730 二逼平衡树 线段树套Treap

    题目描述 您需要写一种数据结构(可参考题目标题),来维护一个有序数列,其中需要提供以下操作:1.查询k在区间内的排名2.查询区间内排名为k的值3.修改某一位值上的数值4.查询k在区间内的前驱(前驱定义 ...

  9. bzoj3196 二逼平衡树 树套树(线段树套Treap)

    Tyvj 1730 二逼平衡树 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 4697  Solved: 1798[Submit][Status][D ...

随机推荐

  1. 封闭解(Closed-form solution)、解析解(Analytical solution)、数值解(Numerical solution) 释义

    转俞夕的博客 (侵删) 1 解析解 解析解(Analytical solution) 就是根据严格的公式推导,给出任意的自变量就可以求出其因变量,也就是问题的解,然后可以利用这些公式计算相应的问题.所 ...

  2. UIApplication直接应用

    /************ 当程序载入后执行,应用程序启动入口 *****************************/ - (BOOL)application:(UIApplication *) ...

  3. loj #2143. 「SHOI2017」组合数问题

    #2143. 「SHOI2017」组合数问题   题目描述 组合数 Cnm\mathrm{C}_n^mC​n​m​​ 表示的是从 nnn 个互不相同的物品中选出 mmm 个物品的方案数.举个例子, 从 ...

  4. webstrom 10 注册码

    webStorm : UserName:William ===== LICENSE BEGIN ===== 45550-12042010 00001SzFN0n1bPII7FnAxnt0DDOPJA  ...

  5. web flash推流h264视频取消cbr

    flash as3.0 推视频H264关于视频质量 初次设置代码如下: h264Setting.setProfileLevel(H264Profile.BASELINE, H264Level.LEVE ...

  6. PAT天梯赛L2-005 集合相似度

    题目链接:点击打开链接 给定两个整数集合,它们的相似度定义为:Nc/Nt*100%.其中Nc是两个集合都有的不相等整数的个数,Nt是两个集合一共有的不相等整数的个数.你的任务就是计算任意一对给定集合的 ...

  7. Qt 学习之路 2(18):事件

    Home / Qt 学习之路 2 / Qt 学习之路 2(18):事件 Qt 学习之路 2(18):事件  豆子  2012年9月27日  Qt 学习之路 2  60条评论 事件(event)是由系统 ...

  8. react 的理念

    命名式的编程方式: 命名式的编程方式,我们会有百分之六七十都在进行dom的操作. 1.声名式的开发: react是面向数据开发的,react是根据这个数据自动构建这个网站,可以把数据理解成图纸,rea ...

  9. P4289 [HAOI2008]移动玩具

    传送门 广搜 4*4 的方阵只有 0 和 1 显然可以状态压缩 (如样例的开始状态压缩后就是1111000011100010) 为了加快速度用了双向广搜(顺便学了一下双向广搜) 双向广搜顾名思义 就是 ...

  10. iOS WebDriverAgent 环境搭建

    WebDriverAgent简介 WebDriverAgent是Facebook 在去年的 SeleniumConf 大会上推出了一款新的iOS移动测试框架.当时的推文申明,还只支持模拟器,不过在今年 ...