Bzoj1901 Dynamic Ranking
动态区间第k小
离散化后
那么每个点开一棵线段树(主席树)再套一个树状数组在外面
每次询问区间内的树的个数时
相当于进行了一次树状数组求区间和的操作,只是是把树状数组那个点看做主席树,对log棵主席树求区间和
然后每次询问,修改时就是把log棵主席树同时跳到儿子,修改也是log棵
时间复杂度O(nlogn*logn)空间复杂度O(nlogn*logn)
# include <bits/stdc++.h>
# define IL inline
# define RG register
# define Fill(a, b) memset(a, b, sizeof(a))
using namespace std;
typedef long long ll;
const int _(2e4 + 10), __(2e6 + 10);
IL ll Read(){
RG char c = getchar(); RG ll x = 0, z = 1;
for(; c < '0' || c > '9'; c = getchar()) z = c == '-' ? -1 : 1;
for(; c >= '0' && c <= '9'; c = getchar()) x = (x << 1) + (x << 3) + (c ^ 48);
return x * z;
}
int n, m, len, p[_], a[_], ql[_], qr[_], qk[_], qo[_];
int ls[__], rs[__], sz[__], rt[__], num, tmp[2][20];
IL void Modify(RG int &x, RG int l, RG int r, RG int pos, RG int val){
if(!x) x = ++num;
sz[x] += val;
if(l == r) return;
RG int mid = (l + r) >> 1;
if(pos <= mid) Modify(ls[x], l, mid, pos, val);
else Modify(rs[x], mid + 1, r, pos, val);
}
IL void PreModify(RG int x, RG int val){
RG int k = lower_bound(p + 1, p + len + 1, a[x]) - p;
for(RG int i = x; i <= n; i += i & -i) Modify(rt[i], 1, len, k, val);
}
IL int Query(RG int l, RG int r, RG int k){
if(l == r) return l;
RG int mid = (l + r) >> 1, sum = 0;
for(RG int i = 1; i <= tmp[1][0]; i++) sum += sz[ls[tmp[1][i]]];
for(RG int i = 1; i <= tmp[0][0]; i++) sum -= sz[ls[tmp[0][i]]];
if(k <= sum){
for(RG int i = 1; i <= tmp[1][0]; i++) tmp[1][i] = ls[tmp[1][i]];
for(RG int i = 1; i <= tmp[0][0]; i++) tmp[0][i] = ls[tmp[0][i]];
return Query(l, mid, k);
}
else{
for(RG int i = 1; i <= tmp[1][0]; i++) tmp[1][i] = rs[tmp[1][i]];
for(RG int i = 1; i <= tmp[0][0]; i++) tmp[0][i] = rs[tmp[0][i]];
return Query(mid + 1, r, k - sum);
}
}
IL int PreQuery(RG int l, RG int r, RG int k){
Fill(tmp, 0);
for(RG int i = r; i; i -= i & -i) tmp[1][++tmp[1][0]] = rt[i];
for(RG int i = l - 1; i; i -= i & -i) tmp[0][++tmp[0][0]] = rt[i];
return Query(1, len, k);
}
int main(RG int argc, RG char* argv[]){
n = Read(); m = Read();
for(RG int i = 1; i <= n; i++) a[i] = Read(), p[++len] = a[i];
for(RG int i = 1; i <= m; i++){
RG char c; scanf(" %c", &c);
qo[i] = c == 'Q';
if(qo[i]) ql[i] = Read(), qr[i] = Read(), qk[i] = Read();
else ql[i] = qr[i] = Read(), qk[i] = Read(), p[++len] = qk[i];
}
sort(p + 1, p + len + 1); len = unique(p + 1, p + len + 1) - p - 1;
for(RG int i = 1; i <= n; i++) PreModify(i, 1);
for(RG int i = 1; i <= m; i++)
if(qo[i]) printf("%d\n", p[PreQuery(ql[i], qr[i], qk[i])]);
else{
PreModify(ql[i], -1);
a[ql[i]] = qk[i];
PreModify(ql[i], 1);
}
return 0;
}
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 探究高空视频全景AR技术的实现原理
· 理解Rust引用及其生命周期标识(上)
· 浏览器原生「磁吸」效果!Anchor Positioning 锚点定位神器解析
· 没有源码,如何修改代码逻辑?
· 一个奇形怪状的面试题:Bean中的CHM要不要加volatile?
· 分享4款.NET开源、免费、实用的商城系统
· 全程不用写代码,我用AI程序员写了一个飞机大战
· Obsidian + DeepSeek:免费 AI 助力你的知识管理,让你的笔记飞起来!
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· 白话解读 Dapr 1.15:你的「微服务管家」又秀新绝活了