题意:带修改不带插入的区间k大。

裸的可持久化线段树。。由于有修改,要用树状数组维护。其它跟不带修改的可持久化线段树一样。

因为我没有找到网上用指针写的代码。。CLJ写这道题也用的不是可持久化线段树,于是我就没有任何模板可以参照。。就参考网上数组版的自己脑补了一个指针版。。你们就可以看到代码优美度下降了好多。

由于数字范围很大,我们需要把所有输入读进来然后离散化。。不离散化的话就要动态开点(这个我暂时还不会)

这份代码用的空间比较多。。在zju是A不了的,需要空间优化(怎么优化我也不知道)

  1 #include<cstdio>
  2 #include<algorithm>
  3 const int MAXN=10000+5;
  4 const int INF=~0U>>1;
  5 const int BUFFER_SIZE=10000;
  6 const int MAXM=10000+5;
  7 struct Tree{
  8     Tree* pl,*pr;
  9     int l,r,sum;
 10     Tree* set(int _l,int _r,Tree* _pl,Tree* _pr)
 11     {
 12         l=_l;r=_r;pl=_pl;pr=_pr;
 13         sum=pl->sum+pr->sum;
 14         return this;
 15     }
 16     Tree* set(int _l,int _r,int all);
 17     Tree* add(int pos,int ad);
 18 };
 19 Tree* buffer=0,*cur;
 20 inline Tree* get()
 21 {
 22     if(!buffer || cur-buffer==BUFFER_SIZE)
 23         buffer=new Tree[BUFFER_SIZE],cur=buffer;
 24     return cur++;
 25 }
 26 Tree* Tree::set(int _l,int _r,int all)
 27 {
 28     l=_l;r=_r;
 29     if(l+1==r) sum=all;
 30     else
 31     {
 32         int m=(l+r)>>1;
 33         pl=get()->set(l,m,all);
 34         pr=get()->set(m,r,all);
 35         sum=pl->sum+pr->sum;
 36     }
 37     return this;
 38 }
 39 Tree* Tree::add(int pos,int ad)
 40 {
 41     if(l+1==r) return get()->set(l,r,sum+ad);
 42     int m=(l+r)>>1;
 43     if(pos<m) return get()->set(l,r,pl->add(pos,ad),pr);
 44     else return get()->set(l,r,pl,pr->add(pos,ad));
 45 }
 46 inline int lowbit(int x)
 47 {    return x&-x;    }
 48 int n,m;
 49 Tree* root[MAXN];
 50 void update(int x,int pos,int ad)
 51 {
 52     for(;x<=n;x+=lowbit(x))
 53         root[x]=root[x]->add(pos,ad);
 54 }
 55 Tree* ptl[MAXN],*ptr[MAXN];
 56 int lnum,rnum;
 57 inline void get_trees(int x,int y)
 58 {
 59     lnum=rnum=0;
 60     for(;x>0;x-=lowbit(x))
 61         ptl[lnum++]=root[x];
 62     for(;y>0;y-=lowbit(y))
 63         ptr[rnum++]=root[y];
 64 }
 65 inline int get_sum()
 66 {
 67     int cnt1=0,cnt2=0;
 68     for(int i=0;i<lnum;++i) cnt1+=ptl[i]->pl->sum;
 69     for(int i=0;i<rnum;++i) cnt2+=ptr[i]->pl->sum;
 70     return cnt2-cnt1;
 71 }
 72 inline int query(int x,int y,int k)
 73 {
 74     get_trees(x-1,y);
 75     while(ptr[0]->l +1 < ptr[0]->r)
 76     {
 77         int cnt=get_sum();
 78         if(cnt<=k)
 79         {
 80             k-=cnt;
 81             for(int i=0;i<lnum;++i) ptl[i]=ptl[i]->pr;
 82             for(int i=0;i<rnum;++i) ptr[i]=ptr[i]->pr;
 83         }
 84         else
 85         {
 86             for(int i=0;i<lnum;++i) ptl[i]=ptl[i]->pl;
 87             for(int i=0;i<rnum;++i) ptr[i]=ptr[i]->pl;
 88         }
 89     }
 90     return ptr[0]->l;
 91 }
 92 struct Ask{
 93     int kind,a,b,c;
 94 }ask[MAXM];
 95 int w[MAXN],sortw[MAXN+MAXM];
 96 int main()
 97 {
 98     freopen("1.in","r",stdin);
 99     scanf("%d%d",&n,&m);
100     int p;
101     for(int i=0;i<n;++i)
102         scanf("%d",w+i),sortw[i]=w[i];
103     p=n;
104     char cmd[5];
105     for(int i=0;i<m;++i)
106     {
107         scanf("%s",cmd);
108         if(cmd[0]=='Q')
109         {
110             ask[i].kind=0;
111             scanf("%d%d%d",&ask[i].a,&ask[i].b,&ask[i].c);
112         }
113         else
114         {
115             ask[i].kind=1;
116             scanf("%d%d",&ask[i].a,&ask[i].b);
117             sortw[p]=ask[i].b;
118             p++;
119         }
120     }
121     std::sort(sortw,sortw+p);
122     int nn=std::unique(sortw,sortw+p)-sortw;
123     for(int i=0;i<n;++i) w[i]=std::lower_bound(sortw,sortw+nn,w[i])-sortw;
124     root[0]=get()->set(0,nn,0);
125     for(int i=1;i<=n;++i) root[i]=root[0]->add(w[i-1],1);
126     for(int i=1;i<=n;++i) update(i+lowbit(i),w[i-1],1);
127     for(int i=0;i<m;++i)
128     {
129         if(!ask[i].kind) printf("%d\n",sortw[query(ask[i].a,ask[i].b,ask[i].c-1)]);
130         else
131         {
132             int t=std::lower_bound(sortw,sortw+nn,ask[i].b)-sortw;
133             update(ask[i].a,w[ask[i].a-1],-1);
134             w[ask[i].a-1]=t;
135             update(ask[i].a,t,1);
136         }
137     }
138     return 0;
139 }

View Code

转载于:https://www.cnblogs.com/lowsfish/p/4435491.html

BZOJ 1901 Zju2112 Dynamic Rankings 题解相关推荐

  1. Bzoj 1901: Zju2112 Dynamic Rankings 主席树,可持久,树状数组,离散化

    1901: Zju2112 Dynamic Rankings Time Limit: 10 Sec  Memory Limit: 128 MB Submit: 6321  Solved: 2628 [ ...

  2. bzoj 1901: Zju2112 Dynamic Rankings(离线树状数组+主席树)

    1901: Zju2112 Dynamic Rankings Time Limit: 10 Sec  Memory Limit: 128 MB Submit: 8144  Solved: 3378 [ ...

  3. bzoj 1901: Zju2112 Dynamic Rankings

    Time Limit: 10 Sec  Memory Limit: 128 MB Submit: 6245  Solved: 2593 [Submit][Status][Discuss] Descri ...

  4. 【BZOJ】1901: Zju2112 Dynamic Rankings(区间第k小+树状数组套主席树)

    http://www.lydsy.com/JudgeOnline/problem.php?id=1901 首先还是吐槽时间,我在zoj交无限tle啊!!!!!!!!我一直以为是程序错了啊啊啊啊啊啊. ...

  5. BZOJ1901:Zju2112 Dynamic Rankings——题解

    http://www.lydsy.com/JudgeOnline/problem.php?id=1901 Description 给定一个含有n个数的序列a[1],a[2],a[3]--a[n],程序 ...

  6. BZOJ1901: Zju2112 Dynamic Rankings

    题解:很显然的树套树 主席树+树状数组 动态维护区间权值 然后类似于查询区间第K大 复杂度nlog^2n 空间也如此 #include <bits/stdc++.h> const int ...

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

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

  8. Zju2112 Dynamic Rankings(树状数组套可持久化权值线段树)

    Zju2112 Dynamic Rankings description solution code description 给定一个含有n个数的序列a[1],a[2],a[3]--a[n],程序必须 ...

  9. 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- ...

  10. 【BZOJ1901】Zju2112 Dynamic Rankings

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

最新文章

  1. 比赛--建金字塔问题--解题报告
  2. zipaligin不是内部或外部命令,也不是可运行的程序
  3. Linux(基础命令和三剑客)
  4. linux ntp连接失败,linux ntp服务器连接异常
  5. Spring MVC —— 前后台传递JSON
  6. 不能在本地计算机启动 dtc,win10系统dTC服务无法启动的操作技巧
  7. 2005年2月24日(星期四) 中午,晴+煙 - Central Incubator。
  8. ASP.NET MVC 在WebService中Token的使用方法
  9. 实战| JSP Servlet Mysql学生信息管理系统
  10. python文件处理,将DNA序列转换为RNA序列
  11. 人生小故事,看完后想十秒钟,你会获益匪浅
  12. 杰控组态西门子PLC
  13. NoSQL数据库知识
  14. MacOS Ventura 13.0.1 (22A400) 正式版带 OC 0.8.6 and winPE 双分区原版黑苹果镜像
  15. TCP和UDP的区别与联系
  16. UCOSII系统移植详解
  17. 显微镜下的大明内容_显微镜下的大明
  18. Android视频列表自动播放功能
  19. Linux printf Class.this git rebase Persistent Bugger
  20. 财务自由到底是啥感觉啊...

热门文章

  1. java并发编程(1)--线程 可见性 volatile怎么用
  2. 2013年广西计算机二级c语言试题及答案,2013年计算机二级C语言上机模拟试题及答案...
  3. linux创建表空间 没有权限,ORA-12913: 无法创建字典管理的表空间
  4. mysql备份 php_详解MYSQL的备份还原(PHP实现)
  5. 列举微型计算机的主要性能指标,计算机基础考试纲.ppt
  6. mysql查询语句内连接、左连接、右连接以及全连接查询
  7. java 视频格式转换-ffmpeg
  8. 定时任务莫名停止,Spring 定时任务存在 Bug??
  9. k8s的精简版k3s安装
  10. Redis主从配置,哨兵,集群的设计原理