标签：now 前缀 int Rankings ll Dynamic include define
本来想用这个题来写一下整体二分。
不过看了很久的树套树神奇数据结构貌似还从来没写过。
来练下手好了。
考虑我们没有修改操作时，区间\(kth\)可以用主席树来做，相当于维护了一个前缀和（这个前缀和的形式是由权值线段树来表达的）。
但是我们有修改操作了，如果我们继续维护这个暴力的前缀和，那么每次修改是\(O（n）\)的。
但是如果我们套上一个树状数组，考虑变成了单点加到区间查询就变得非常容易，这样的代价是我们从\(O(log)\)的一次查询，变成每次要同时跑\(log\)颗树。
是一个\(O(nlog^2n)\)的算法。
Dynamic Rankings
#include<iostream>
#include<algorithm>
#include<cstdio>
#define ll long long
#define N 200005
#define lowbit(x) (x & -x)
ll n,m;
ll a[N],num[N << 1];//数集
struct P{
int opt;
ll l,r,k;
}e[N];
ll head[N];//树状数组对应的线段树的头结点。
struct E{ll ls,rs,val;}T[N * 400];
#define l(x) T[x].ls
#define r(x) T[x].rs
#define v(x) T[x].val
#define mid ((l + r) >> 1)
ll cnt;
inline void add(ll &now,int l,int r,int to,int ki){
if(!now)now = ++cnt;
// std::cout<<now<<" "<<l<<" "<<r<<" "<<to<<" "<<ki<<std::endl;
v(now) += ki;
if(l == r)return;
if(to <= mid)
add(l(now),l,mid,to,ki);
else
add(r(now),mid + 1,r,to,ki);
}
inline void change(int now,int val){
ll k = std::lower_bound(num + 1,num + num[0] + 1,a[now]) - num;
for(int i = now;i <= num[0];i += lowbit(i))
add(head[i],1,num[0],k,val);
}
ll t1[N << 2],t2[N << 2];
inline ll find(int l,int r,int k){
// std::cout<<l<<" "<<r<<" "<<k<<std::endl;
if(l == r)return l;
ll sum = 0;
for(int i = 1;i <= t2[0];++i)
sum += v(l(t2[i]));
for(int i = 1;i <= t1[0];++i)
sum -= v(l(t1[i]));
// std::cout<<sum<<std::endl;
if(sum >= k){
for(int i = 1;i <= t1[0];++i)
t1[i] = l(t1[i]);
for(int i = 1;i <= t2[0];++i)
t2[i] = l(t2[i]);
return find(l,mid,k);
}
else{
for(int i = 1;i <= t1[0];++i)
t1[i] = r(t1[i]);
for(int i = 1;i <= t2[0];++i)
t2[i] = r(t2[i]);
return find(mid + 1,r,k - sum);
}
}
inline ll pre(ll l,ll r,ll k){
t1[0] = 0;
t2[0] = 0;
for(int i = l - 1;i;i -= lowbit(i))
t1[++t1[0]] = head[i];
for(int i = r;i;i -= lowbit(i))
t2[++t2[0]] = head[i];
// for(int i = 1;i <= t1[0];++i)
// std::cout<<t1[i]<<" ";
// puts("");
// for(int i = 1;i <= t2[0];++i)
// std::cout<<t2[i]<<" ";
// puts("");
return find(1,num[0],k);
}
int main(){
scanf("%lld%lld",&n,&m);
for(int i = 1;i <= n;++i)
scanf("%lld",&a[i]),num[++num[0]] = a[i];
for(int i = 1;i <= m;++i){
char a;
while(a != 'Q' && a != 'C')
a = getchar();
if(a == 'Q'){
e[i].opt = 1;
scanf("%lld%lld%lld",&e[i].l,&e[i].r,&e[i].k);
}else{
e[i].opt = 2;
scanf("%lld%lld",&e[i].k,&e[i].l);
num[++num[0]] = e[i].l;
}
a = 'z';
}
std::sort(num + 1,num + num[0] + 1);
num[0] = std::unique(num + 1,num + num[0] + 1) - num - 1;
for(int i = 1;i <= n;++i)
change(i,1);
// for(int i = 1;i <= n;++i)
// std::cout<<head[i]<<std::endl;
for(int i = 1;i <= m;++i){
if(e[i].opt == 1){
std::cout<<num[pre(e[i].l,e[i].r,e[i].k)]<<std::endl;
}else{
change(e[i].k,-1);
a[e[i].k] = e[i].l;
change(e[i].k,1);
}
}
}
标签：now,前缀,int,Rankings,ll,Dynamic,include,define 来源： https://www.cnblogs.com/dixiao/p/14776438.html
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。