来自FallDream的博客,未经允许,请勿转载,谢谢。


维护一个长度为n的正整数序列a_1,a_2,...,a_n,支持以下两种操作:

1 k,将序列a从小到大排序,输出a_k的值。
2 k,将所有严格大于k的数a_i减去k。
n,m<=100000
 
splay维护序列,每次操作小等于k的数字不变,大于2k的数字顺序不变,(k,2k]的数字减小至少一半,最多减小log次,暴力拿出来重新插入就行了。
复杂度nlog^2n
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#define Rint register int
#define ll long long
#define MN 100000
#define INF 1000000000000000000LL
#define getchar() (*S++)
char B[<<],*S=B;
using namespace std;
inline int read()
{
int x = , f = ; char ch = getchar();
while(ch < '' || ch > ''){ if(ch == '-') f = -; ch = getchar();}
while(ch >= '' && ch <= ''){x = x * + ch - '';ch = getchar();}
return x * f;
}
int c[MN+][],fa[MN+],rt=,tag[MN+],n,m,size[MN+],q[MN+],top=;
ll a[MN+];
inline void update(int x){size[x]=size[c[x][]]+size[c[x][]]+;}
inline void Mark(int x,int v){tag[x]+=v;a[x]+=v;}
inline void pushdown(int x){Mark(c[x][],tag[x]);Mark(c[x][],tag[x]);tag[x]=;} void rotate(int x,int&k)
{
int y=fa[x],z=fa[y],l=c[y][]==x,r=l^;
if(y==k) k=x; else c[z][c[z][]==y]=x;
fa[x]=z;fa[y]=x;fa[c[x][r]]=y;
c[y][l]=c[x][r];c[x][r]=y;
update(y);update(x);
} void Splay(int x,int&k)
{
for(;x!=k;rotate(x,k))
if(fa[x]!=k) rotate((c[fa[fa[x]]][]==fa[x]^c[fa[x]][]==x)?x:fa[x],k);
} void build(int&x,int l,int r,int last)
{
if(l>r) {x=;return;}x=l+r>>;fa[x]=last;
build(c[x][],l,x-,x);
build(c[x][],x+,r,x);
update(x);
} int Ask_After(int x,int v)
{
if(!x) return ;int q;
if(tag[x]) pushdown(x);
if(a[x]>v) return (q=Ask_After(c[x][],v))?q:x;
else return Ask_After(c[x][],v);
} int Ask_Before(int x,int v)
{
if(!x) return ;int q;
if(tag[x]) pushdown(x);
if(a[x]<=v) return (q=Ask_Before(c[x][],v))?q:x;
else return Ask_Before(c[x][],v);
} int Find(int x,int rk)
{
int Sz=size[c[x][]]+;
if(tag[x]) pushdown(x);
if(Sz==rk) return x;
if(Sz<rk) return Find(c[x][],rk-Sz);
else return Find(c[x][],rk);
} void Dfs(int x)
{
if(tag[x]) pushdown(x);
if(c[x][]) Dfs(c[x][]);
q[++top]=x;
if(c[x][]) Dfs(c[x][]);
} void Ins(int&x,int k,int last)
{
if(!x){x=k;c[k][]=c[k][]=;fa[k]=last;size[k]=;return;}
if(tag[x]) pushdown(x);
Ins(c[x][a[k]>=a[x]],k,x);
update(x);
} int main()
{
fread(B,,<<,stdin);
n=read();m=read();
for(Rint i=;i<=n;++i) a[i+]=read();a[]=-INF;a[n+]=INF;
sort(a+,a+n+);build(rt,,n+,);
for(Rint i=;i<=m;++i)
{
int op=read(),v;
if(op==){printf("%d\n",(v=Find(rt,read()+),Splay(v,rt),a[v]));}
if(op==)
{
v=read();int x=Ask_Before(rt,v),y=Ask_After(rt,v<<);
Splay(x,rt);Splay(y,c[rt][]);int z=c[y][];
if(z) {top=;Dfs(z);fa[z]=;c[y][]=;}
update(y);update(x);Mark(y,-v);
if(z) for(Rint j=;j<=top;++j) a[q[j]]-=v,Ins(rt,q[j],),Splay(q[j],rt);
}
}
return ;
}

[bzoj4923]K小值查询的更多相关文章

  1. BZOJ4923 K小值查询(splay)

    容易想到建一棵平衡树,修改时打上标记即可.但是修改会导致平衡树结构被破坏.注意到实际上只有[k+1,2k)这一部分数在平衡树中的位置会被改变,所以对这一部分暴力修改,因为每次都会使其至少减小一半,复杂 ...

  2. [BZ4923][Lydsy1706月赛]K小值查询

    K小值查询 题面 维护一个长度为n的正整数序列a_1,a_2,...,a_n,支持以下两种操作: 1 k,将序列a从小到大排序,输出a_k的值. 2 k,将所有严格大于k的数a_i减去k. Input ...

  3. BZOJ4923:[Lydsy1706月赛]K小值查询(Splay)

    Description 维护一个长度为n的正整数序列a_1,a_2,...,a_n,支持以下两种操作: 1 k,将序列a从小到大排序,输出a_k的值. 2 k,将所有严格大于k的数a_i减去k. In ...

  4. BZOJ 4923: [Lydsy1706月赛]K小值查询 Splay + 思维

    Description 维护一个长度为n的正整数序列a_1,a_2,...,a_n,支持以下两种操作: 1 k,将序列a从小到大排序,输出a_k的值. 2 k,将所有严格大于k的数a_i减去k. In ...

  5. BZOJ4923 [Lydsy1706月赛]K小值查询

    题意 维护一个长度为n的正整数序列a_1,a_2,...,a_n,支持以下两种操作: 1 k,将序列a从小到大排序,输出a_k的值. 2 k,将所有严格大于k的数a_i减去k. \(n \leq 10 ...

  6. 4923: [Lydsy1706月赛]K小值查询 平衡树 非旋转Treap

    国际惯例的题面:这种维护排序序列,严格大于的进行操作的题都很套路......我们按照[0,k],(k,2k],(2k,inf)分类讨论一下就好.显然第一个区间的不会变化,第二个区间的会被平移进第一个区 ...

  7. [BZOJ 4923][Lydsy1706月赛]K小值查询

    传送门 势能分析平衡树,splay或treap都可以 放个指针版的就跑 #include <bits/stdc++.h> using namespace std; #define rep( ...

  8. 【BZOJ】3065: 带插入区间K小值

    http://www.lydsy.com/JudgeOnline/problem.php?id=3065 题意:带插入.修改的区间k小值在线查询.(原序列n<=35000, 询问<=175 ...

  9. bzoj 3065: 带插入区间K小值 替罪羊树 && AC300

    3065: 带插入区间K小值 Time Limit: 60 Sec  Memory Limit: 512 MBSubmit: 1062  Solved: 253[Submit][Status] Des ...

随机推荐

  1. Jmeter入门(01)Jmeter的下载和安装

    一.什么是Jmeter 1.一款优秀的.开源的.免费的.功能测试和性能测试 工具 Jmeter ,使用Java开发的一款优秀的开源免费测试工具,主要用来做功能测试和性能测试(压力测试/负载测试),用J ...

  2. OAuth2.0学习(1-13)oauth2.0 的概念:资源、权限(角色)和scope

    mkk 关于资源的解释 : https://andaily.com/blog/?cat=19 resource用于将系统提供的各类资源进行分组管理, 每一个resource对应一个resource-i ...

  3. apigw鉴权分析(1-2)腾讯开放平台 - 鉴权分析

    一.访问入口 http://wiki.open.qq.com/wiki/%E8%85%BE%E8%AE%AF%E5%BC%80%E6%94%BE%E5%B9%B3%E5%8F%B0%E7%AC%AC% ...

  4. Android TabLayout 在宽屏幕上tab不能平均分配的问题解决

    当TabLayout 在宽屏幕的设备上,如平板横屏的时候,tab的宽度超过一定值后,就不在平均分配宽度,而是居中显示.此时设置 app:tabMode="fixed"或者 top_ ...

  5. golang微信公众号请求获取信息

    初次用golang在公众号中获取信息,记录一下 看了下文档,粗略的写了个demo,如下: func HttpGet(c*gin.Context) { var param GetType if er:= ...

  6. java字符串类型常量拼接与变量拼接的区别

    前言 首先看下下面代码结果是什么? package cn.demo_01; public class StringDemo02 { public static void main(String[] a ...

  7. hue集成hbase出现TSocket read 0 bytes

    解决办法:修改hbase的配置文件 添加以下配置 https://stackoverflow.com/questions/20415493/api-error-tsocket-read-0-bytes ...

  8. C# 打开文件夹和保存文件夹

    OpenFileDialog sfd = new OpenFileDialog(); sfd.ShowDialog(); this.textBox1.Text = sfd.FileName; Save ...

  9. java 反射(Reflection)

    看了很多关于java 反射的文章,自己把所看到的总结一下.对自己,对他人或多或少有帮助吧. Java Reflection是什么? 首先来看看官方文档Oracle里面对Reflection的描述: R ...

  10. python3全栈开发- 元类metaclass(面试必考题)

    一.知识储备 #exec:三个参数 #参数一:字符串形式的命令 #参数二:全局作用域(字典形式),如果不指定,默认为globals() #参数三:局部作用域(字典形式),如果不指定,默认为locals ...