BZOJ4923:[Lydsy1706月赛]K小值查询(Splay)
Description
Input
Output
Sample Input
1 5 6 12
2 5
1 1
1 2
1 3
1 4
Sample Output
1
5
7
Solution
把数排个序然后建$Splay$,每次修改对值域为$[1,k]$中的不管,$[k+1,k\times 2]$中的拆出来改完了再暴力插回去,对于$[k\times 2+1,MAX]$中打标记。我也不知道为什么复杂度是对的。
以后别有事没事把标记下传到$0$点,修改着修改着$0$下标的值就不知道被修改成什么鬼畜的数了。
Code
#include<iostream>
#include<cstdio>
#include<algorithm>
#define N (100009)
using namespace std; int n,m,opt,k,a[N];
int Root,Father[N],Son[N][];
int Val[N],Size[N],Max[N],Add[N]; int Get(int x) {return Son[Father[x]][]==x;} void Pushup(int x)
{
Size[x]=Size[Son[x][]]+Size[Son[x][]]+;
Max[x]=max(Val[x],max(Max[Son[x][]],Max[Son[x][]]));
} int Build(int fa,int l,int r)
{
if (l>r) return ;
int mid=(l+r)>>;
Father[mid]=fa; Val[mid]=a[mid];
Son[mid][]=Build(mid,l,mid-);
Son[mid][]=Build(mid,mid+,r);
Pushup(mid); return mid;
} void Rotate(int x)
{
int wh=Get(x);
int fa=Father[x],fafa=Father[fa];
if (fafa) Son[fafa][Son[fafa][]==fa]=x;
Father[fa]=x; Son[fa][wh]=Son[x][wh^];
if (Son[fa][wh]) Father[Son[fa][wh]]=fa;
Father[x]=fafa; Son[x][wh^]=fa;
Pushup(fa); Pushup(x);
} void Pushdown(int x)
{
if (Add[x]!=)
{
if (Son[x][])
{
Val[Son[x][]]+=Add[x];
Add[Son[x][]]+=Add[x];
Max[Son[x][]]+=Add[x];
}
if (Son[x][])
{
Val[Son[x][]]+=Add[x];
Add[Son[x][]]+=Add[x];
Max[Son[x][]]+=Add[x];
}
Add[x]=;
}
} void Push(int x)
{
if (Father[x]) Push(Father[x]);
Pushdown(x);
} void Splay(int x,int tar)
{
Push(x);
for (int fa; (fa=Father[x])!=tar; Rotate(x))
if (Father[fa]!=tar)
Rotate(Get(fa)==Get(x)?fa:x);
if (!tar) Root=x;
} int Findkth(int x)
{
int now=Root;
while ()
{
Pushdown(now);
if (Size[Son[now][]]>=x) now=Son[now][];
else
{
x-=Size[Son[now][]];
if (x==) {Splay(now,); return Val[now];}
x--; now=Son[now][];
}
}
} int Find(int x)
{
int now=Root,ans=;
while ()
{
Pushdown(now);
if (Max[Son[now][]]>x) now=Son[now][];
else
{
if (Val[now]>x) {Splay(now,); return now;}
now=Son[now][];
}
}
} int Pre(int x)
{
Splay(x,);
x=Son[x][];
while (Son[x][]) x=Son[x][];
return x;
} void Insert(int x)
{
int now=Root,fa=;
while ()
{
Pushdown(now);
if (!now)
{
Father[x]=fa;
Son[fa][Val[fa]<Val[x]]=x;
Max[x]=Val[x]; Size[x]=;
Splay(x,); return;
}
fa=now, now=Son[now][Val[now]<Val[x]];
}
} void DFS(int x,int k)
{
if (!x) return;
Pushdown(x);
DFS(Son[x][],k); DFS(Son[x][],k);
Father[x]=Son[x][]=Son[x][]=Size[x]=Max[x]=;
Val[x]-=k; Insert(x);
} void Update(int k)
{
int x=Pre(Find(k)),y=Find(*k);
Splay(x,); Splay(y,x);
int s=Son[y][]; Son[y][]=Father[Son[y][]]=;
DFS(s,k);
Splay(y,);
x=Pre(Find(*k)),y=n+;
Splay(x,); Splay(y,x);
if (!Son[y][]) return;
Add[Son[y][]]-=k; Val[Son[y][]]-=k; Max[Son[y][]]-=k;
} int main()
{
scanf("%d%d",&n,&m);
for (int i=; i<=n+; ++i)
scanf("%d",&a[i]);
a[]=-2e9; a[n+]=2e9; Max[]=-2e9;
sort(a+,a+n+);
Root=Build(,,n+);
for (int i=; i<=m; ++i)
{
scanf("%d%d",&opt,&k);
if (opt==) printf("%d\n",Findkth(k+));
else Update(k);
}
}
BZOJ4923:[Lydsy1706月赛]K小值查询(Splay)的更多相关文章
- 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 ...
- BZOJ4923 [Lydsy1706月赛]K小值查询
题意 维护一个长度为n的正整数序列a_1,a_2,...,a_n,支持以下两种操作: 1 k,将序列a从小到大排序,输出a_k的值. 2 k,将所有严格大于k的数a_i减去k. \(n \leq 10 ...
- [BZ4923][Lydsy1706月赛]K小值查询
K小值查询 题面 维护一个长度为n的正整数序列a_1,a_2,...,a_n,支持以下两种操作: 1 k,将序列a从小到大排序,输出a_k的值. 2 k,将所有严格大于k的数a_i减去k. Input ...
- 4923: [Lydsy1706月赛]K小值查询 平衡树 非旋转Treap
国际惯例的题面:这种维护排序序列,严格大于的进行操作的题都很套路......我们按照[0,k],(k,2k],(2k,inf)分类讨论一下就好.显然第一个区间的不会变化,第二个区间的会被平移进第一个区 ...
- [BZOJ 4923][Lydsy1706月赛]K小值查询
传送门 势能分析平衡树,splay或treap都可以 放个指针版的就跑 #include <bits/stdc++.h> using namespace std; #define rep( ...
- BZOJ4923 K小值查询(splay)
容易想到建一棵平衡树,修改时打上标记即可.但是修改会导致平衡树结构被破坏.注意到实际上只有[k+1,2k)这一部分数在平衡树中的位置会被改变,所以对这一部分暴力修改,因为每次都会使其至少减小一半,复杂 ...
- [bzoj4923]K小值查询
来自FallDream的博客,未经允许,请勿转载,谢谢. 维护一个长度为n的正整数序列a_1,a_2,...,a_n,支持以下两种操作: 1 k,将序列a从小到大排序,输出a_k的值. 2 k,将所有 ...
- 【BZOJ】3065: 带插入区间K小值
http://www.lydsy.com/JudgeOnline/problem.php?id=3065 题意:带插入.修改的区间k小值在线查询.(原序列n<=35000, 询问<=175 ...
- 「BZOJ3065」带插入区间第K小值 替罪羊树×线段树
题目描述 从前有\(n\)只跳蚤排成一行做早操,每只跳蚤都有自己的一个弹跳力\(a_i\).跳蚤国王看着这些跳蚤国欣欣向荣的情景,感到非常高兴.这时跳蚤国王决定理性愉悦一下,查询区间\(k\)小值.他 ...
随机推荐
- win2d 图片水印
本文告诉大家如何使用 win2d 给图片加上水印. 安装 首先需要使用 Nuget 安装 win2d ,安装参见win10 uwp win2d 如果没有更新 dot net core 那么在运行可能会 ...
- RabbitMQ远程不能访问问题的解决
刚刚安装的RabbitMQ-Server-3.3.5,并且也已经开启了Web管理功能,但是现在存在一个问题: 出于安全的考虑,guest这个默认的用户只能通过http://localhost:1567 ...
- .NET Core类库中读取配置文件
最近在开发基于.NET Core的NuGet包,遇到一个问题:.NET Core中已经没有ConfigurationManager类,在类库中无法像.NET Framework那样读取App.conf ...
- 学习c++的五十条忠告
1.把C++当成一门新的语言学习: 2.看<Thinking In C++>,不要看<C++变成死相>: 3.看<The C++ Programming Language ...
- TF-IDF原理
什么是TF-IDF TF-IDF(Term Frequency-Inverse Document Frequency, 词频-逆文件频率). 是一种用于资讯检索与资讯探勘的常用加权技术.TF-IDF ...
- Java高级工程师需要弄明白的20个知识点
一般的程序员或许只需知道一些JAVA的语法结构,能对数据库数据进行CRUD就可以应付了.但要成为JAVA(高级) 工程师,就要对JAVA做比较深入的研究,需要不断学习进步,以下对高级工程师需要突破的知 ...
- a 链接点击下载
1. 将链接设置为.zip 结尾2.在a元素中添加download 属性,(目前只有chrome.firefox和opera支持) function download(src) { var $a = ...
- 解决myeclipse2017安装后闪退问题
修改myeclipse的配置文件myeclipse.ini: 1.将-vm下面的路径改成自己的jdk的javaw.exe的路径 2.在文档最后加上两条语句: -Dgenuitec.honorDevMo ...
- NOI.AC NOIP2018 全国热身赛 第四场
心路历程 预计得分:\(0 + 100 +100\) 实际得分:\(10 + 100 + 0\) 神TM T3模数为啥是\(1e9 + 9\)啊啊啊啊,而且我也确实是眼瞎...真是血的教训啊.. T2 ...
- js-ES6学习笔记-Class(5)
1.原生构造函数会忽略apply方法传入的this,也就是说,原生构造函数的this无法绑定,导致拿不到内部属性.比如,Array构造函数有一个内部属性[[DefineOwnProperty]],用来 ...