【大杀器】利用划分树秒杀区间内第k大的数
最近看了一道题,大概就是给出一个序列,不断询问其子区间内第k大的数,下面是个截图
绕了一圈没找到中文版题目,if(你是大佬) then 去看截图;else{我来解释:给出一个整数n,和一个整数m,分别表示序列元素个数和询问次数,然后输入n个数和m个询问,每个询问包含3个数,分别是区间起止点(l和r)和k,求出区间内第k大的数并输出;}这是一道很简单的模板题,怎么解决呢?小编最初想到的是打暴力,正所谓暴力出奇迹,说不定可以成功,反正不会优化,先试试吧,直接把规定区间[l,r]排一次序了,然后在查找一遍第k大的数,但是毫无疑问,绝对会超时,怎样能减少时间复杂度呢?这是就请出了二分。
二分运用了分治的思想,不断将子区间分成两半,直到找到第k大的数,虽然很高效,但是面对这道题的多次询问,二分也只能表示手软,仍然过不了这道题。但是如果换成了线段树的结构,效率则会快很多,线段树看起来也用了分治的思想,把原来的整个序列都大约相等的长度分到左子树和右子树,不断分下去,直到全部分成叶子节点,在逐次确定下一层第k大的数在左子树还是右子树,虽然这种树能成功AC,但是并不是最优的,做题不能只讲求AC,下面就来讲一讲线段树的升级版——划分树。
什么是划分树?
划分树是一种基于线段树的数据结构,也利用了分治的思想,却比线段树高效很多,这是为什么?因为划分树又多了一个性质:在划分时不是随意划分,也不是排序后直接划分(因为这样会破坏原有结构),而是排序后仍保持原来的相对顺序再分到左右子树。
具体实现方法:
整个过程分为建树和查询两个阶段:
//copy来的图
1)建树:首先定义一个数组tree[30][1000]第一个维度表示层数,第二个维度表示这一层第i个数的值,用来表示这棵划分树,然后定义sorted[1000]数组,用来存储排序好的原序列,然后记录每一层前i个数有多少进入了下一层的左子树,存在toleft[30][1000]数组中,在建树中没用,但记录下来对查找时有用,用分治的思想分配左子树和右子树,将不大于中间值mid的数分配到左子树中去,大于中间值的分配到右子树中去,但有时为了左右尽可能个数相等,要把等于中间值的数两边都分配,于是定义same来存储多少等于中间值的数进入左子树,分配完毕后再递归分配左子树和右子树。身为递归,怎么也要有个出口吧,递归到叶子节点时就返回,即if(l==r) return;
2)查询:按照之前存储下的进入下一层左子树个数的数组toleft,可以计算出区间内第k大的数在当前节点的左子树还是右子树,并计算出下一层相应子树的左右边界,然后递归相应子树,同上,递归出口也是到达叶子节点时返回。详见注释……
废话不多说,代码呈上:
#include<iostream>
#include<algorithm>
using namespace std;
int tree[][],sorted[],toleft[][],n,m,ans;//tree和toleft的两个维度分别存储深度和序列,sorted存储的是排序好的序列
void buildtree(int l,int r,int dep)//构建划分树
{
if(l==r) return;//遇见叶子节点就返回
int mid=(l+r)/;//二分
int same=mid-l+;//same最终保存的是和中间值相同元素的个数,以便确定分到哪一区间
for(int i=l;i<=r;i++)
if(tree[dep][i]<sorted[mid]) same--;
int lpos=l;int rpos=mid+;//左指针和右指针,并非常用的指针,是用来保存现在各区间内元素个数
for(int i=l;i<=r;i++)
{
if(tree[dep][i]<sorted[mid])//小于中间值
tree[dep+][lpos++]=tree[dep][i];//分配到左子区间
else if(tree[dep][i]==sorted[mid]&&same>)//等于中间值且相同个数大于0
{
same--;
tree[dep+][lpos++]=tree[dep][i];//分配到右子区间
}
else tree[dep+][rpos++]=tree[dep][i];//剩下的分配到右子区间
toleft[dep][i]=toleft[dep][l-]+lpos-l;//toleft数组记录这一层前i个数有多少个进入下一层的左子区间,查询时有用
}
buildtree(l,mid,dep+);//构建左子区间(左子树)
buildtree(mid+,r,dep+);//构建右子区间(右子树)
}
int search(int L,int R,int l,int r,int dep,int k)//查询第k大的数
{
if(l==r) return tree[dep][l];//查询到符合要求的叶子节点,返回相应的值
int mid=(L+R)/;//L,R为大区间(主要是每个左子树,右子树的边界)
int cnt=toleft[dep][r]-toleft[dep][l-];//求出[l,r]区间内有多少数进入下一层左子区间
if(cnt>=k)//第k大的数对应节点在左子树
{
int newl=L+toleft[dep][l-]-toleft[dep][L-];//求出下一层第k大的数所在区间边界
int newr=newl+cnt-;
return search(L,mid,newl,newr,dep+,k);
}
else//在右子树
{
int newr=r+toleft[dep][R]-toleft[dep][r];//求出下一层第k大的数所在区间边界
int newl=newr-(r-l-cnt);
return search(mid+,R,newl,newr,dep+,k-cnt);
}
}
int main()
{
cin>>n>>m;
for(int i=;i<=n;i++)
{
cin>>tree[][i];
sorted[i]=tree[][i];
}
sort(sorted+,sorted+n+);//为sorted数组排序
buildtree(,n,);//建树
int a,b,c;
for(int i=;i<=m;i++)
{
cin>>a>>b>>c;//输入询问
cout<<search(,n,a,b,,c)<<endl;//查询第k大的数
}
return ;
}
//额~好像忘了改scanf和printf了,没过别怪我……
【大杀器】利用划分树秒杀区间内第k大的数的更多相关文章
- 利用划分树求解整数区间内第K大的值
如何快速求出(在log2n的时间复杂度内)整数区间[x,y]中第k大的值(x<=k<=y)? 其实我刚开始想的是用快排来查找,但是其实这样是不行的,因为会破坏原序列,就算另外一个数组来存储 ...
- zoj 2112 动态区间求第k大
题目大意: 动态单点更新,然后多次询问求区间内第k大 这里单个的主席树不能实现,这里采取的是树状数组套主席树 首先可以想的是将静态主席树先构建好,不去动它,这里空间复杂度就是O(nlogn),这个只要 ...
- HDU 3473 Minimum Sum (划分树求区间第k大带求和)(转)
题意:在区间中找一个数,求出该区间每个数与这个数距离的总和,使其最小 找的数字是中位数(若是偶数个,则中间随便哪个都可)接着找到该区间比此数大的数的总和 区间中位数可以使用划分树,然后在其中记录:每层 ...
- 离群点检测与序列数据异常检测以及异常检测大杀器-iForest
1. 异常检测简介 异常检测,它的任务是发现与大部分其他对象不同的对象,我们称为异常对象.异常检测算法已经广泛应用于电信.互联网和信用卡的诈骗检测.贷款审批.电子商务.网络入侵和天气预报等领域.这些异 ...
- 一文读懂机器学习大杀器XGBoost原理
http://blog.itpub.net/31542119/viewspace-2199549/ XGBoost是boosting算法的其中一种.Boosting算法的思想是将许多弱分类器集成在一起 ...
- [NewLife.XCode]反向工程(自动建表建库大杀器)
NewLife.XCode是一个有10多年历史的开源数据中间件,支持nfx/netstandard,由新生命团队(2002~2019)开发完成并维护至今,以下简称XCode. 整个系列教程会大量结合示 ...
- [csu/coj 1080]划分树求区间前k大数和
题意:从某个区间内最多选择k个数,使得和最大 思路:首先题目给定的数有负数,如果区间前k大出现负数,那么负数不选和更大,于是对于所有最优选择,负数不会出现,所以用0取代负数,问题便转化为区间的前k大数 ...
- 使用docker-compose 大杀器来部署服务 上
使用docker-compose 大杀器来部署服务 上 我们都听过或者用过 docker,然而使用方式却是仅仅用手动的方式,这样去操作 docker 还是很原始. 好吧,可能在小白的眼中噼里啪啦的对着 ...
- zoj2112 树状数组+主席树 区间动第k大
Dynamic Rankings Time Limit: 10000MS Memory Limit: 32768KB 64bit IO Format: %lld & %llu Subm ...
随机推荐
- Go语言的并发和并行
不知道你有没有注意到,这段代码如果我跑在两个goroutines里面的话: package main import ( "fmt" ) func loop(done chan bo ...
- debussy与modelsim的联调设置
前段时间看到网上有人在使用debussy软件对Verilog代码进行调试,而且都称赞其是多么的好用,看着很是馋人,说吧,现在用的是quartus与modelsim的联调,似乎还是可以的,但就是每次稍微 ...
- 【NOIP】提高组2012 同余方程
[算法]扩展欧几里德算法 [题解]学完扩欧就可以随便水了... 转化为不定方程ax-by=1. 因为1且题目保证有解,所以方程有唯一解. 紫书曰:同余方程的一个解其实指的是一个同余等价类. 所以满足x ...
- 关于SQL注入的五大报错注入函数
~全部都以查user()为例子~ 1.floor()id = 1 and (select 1 from (select count(*),concat(version(),floor(rand(0) ...
- 代码合并:Merge、Rebase 的选择
图解 Git 命令 基本用法 上面的四条命令在工作目录.stage 缓存(也叫做索引)和 commit 历史之间复制文件. git add files 把工作目录中的文件加入 stage 缓存 git ...
- python基础===【爬虫】爬虫糗事百科首页图片代码
import requests import re import urllib.request def getHtml(url): page = requests.get(url) html = pa ...
- 终于解决了Linux下运行OCCI程序一直报Error while trying to retrieve text for error ORA-01804错误
终于解决了Linux下运行OCCI程序一直报Error while trying to retrieve text for error ORA-01804错误 http://blog.csdn.net ...
- 【UOJ#169】元旦老人与数列
论文题. 考虑到这题的维护和区间操作是反向的,也就是说无法像V那题快速的合并标记. 我们知道,一个区间的最小值和其他值是可以分开来维护的,因为如果一个区间被整体覆盖,那么最小值始终是最小值. 对于被覆 ...
- UCenter创始人、Discuz!创始人、管理员账号的认知(转)
UCenter创始人.Discuz!创始人.管理员账号的认知 什么是创始人?现在可能还有好多的站长对这个概念有点模糊,今天我给大家屡屡思路,讲讲这个概念性问题,没啥技术含量.已经明白这个概 ...
- Jmeter接口测试示例
如果是Web,需要使用badboy进行录制,今天讲的是接口,因此可以不用录制. (1)新建测试计划 (2)添加http请求默认值 (3)添加http信息头管理器 (4)添加token的正则表达式:&q ...