POJ 2104为例

思想:

利用归并排序的思想:

  • 建树过程和归并排序类似,每个数列都是子树序列的合并与排序。
  • 查询过程,如果所查询区间完全包含在当前区间中,则直接返回当前区间内小于所求数的元素个数,否则递归的对子树进行求解并相加。
  • 使用STL中的merge对子序列进行合并及排序。
  • 时间复杂度O(nlogn+mlog3n)

代码(vector实现):

#include<cstdio>
#include<algorithm>
#include<vector>
#include<iostream>
using namespace std;//[)
const int maxn = 1000010;
vector<int>dat[maxn*3];
int sorted[maxn], a[maxn];
void build(int k, int l, int r)
{
if(r - l == 1) dat[k].push_back(a[l]);
else {
int lch = 2 * k + 1, rch = 2 * k + 2;
build(lch, l, (l + r)/2);
build(rch, (l + r)/2, r);
dat[k].resize(r - l);
merge(dat[lch].begin(), dat[lch].end(), dat[rch].begin(), dat[rch].end(),dat[k].begin());
}
}
int query(int l, int r, int x, int k, int L, int R)
{
if(r <= L|| l >= R) return 0;//不相交
else if(l <= L && R <= r){//完全包含
return lower_bound(dat[k].begin(), dat[k].end(), x) - dat[k].begin();
}
else {
int lch = 2 * k + 1, rch = 2 * k + 2;
int lsum = query(l, r, x, lch, L, (L + R)/2);
int rsum = query(l, r, x, rch, (L + R)/2, R);
return lsum + rsum;
}
}
int main (void)
{
int n, m;scanf("%d%d",&n,&m);
for(int i = 0; i < n; i++){
scanf("%d",&a[i]);
sorted[i] = a[i];
}
sort(sorted, sorted+n);
build(0, 0, n);
int tl, tr, k;
while(m--){
scanf("%d%d%d",&tl,&tr,&k);
int l = 0, r = n;
while(r - l >1){
int mid = l + (r - l)/2;
int c = query(tl-1, tr, sorted[mid], 0, 0 ,n);
if(c <= k - 1 ) l = mid;
else r = mid;
}
printf("%d\n", sorted[l]);
}
return 0;
}//6000+ms

代码(数组实现):

#include<cstdio>
#include<algorithm>
#include<vector>
#include<iostream>
using namespace std;//[)
const int maxn = 1000010;
int dat[20][3*maxn];
int sorted[maxn], a[maxn];
void build(int p, int l, int r)
{
if(r - l == 1) dat[p][l] = a[l];
else {
build(p+1, l, (l + r)/2);
build(p+1, (l + r)/2, r);
merge(dat[p+1] + l, dat[p+1] + (l + r)/2, dat[p+1] + (l + r)/2, dat[p+1] + r, dat[p]+l);
}
}
int query(int l, int r, int x, int p, int L, int R)
{
if(r <= L|| l >= R) return 0;//不相交
else if(l <= L && R <= r){//完全包含
return lower_bound(dat[p]+L, dat[p]+R, x) - (dat[p]+L);
}
else {
int lsum = query(l, r, x, p+1, L, (L + R)/2);
int rsum = query(l, r, x, p+1, (L + R)/2, R);
return lsum + rsum;
}
}
int main (void)
{
int n, m;scanf("%d%d",&n,&m);
for(int i = 0; i < n; i++){
scanf("%d",&a[i]);
sorted[i] = a[i];
}
sort(sorted, sorted+n);
build(0, 0, n);
int tl, tr, k;
while(m--){
scanf("%d%d%d",&tl,&tr,&k);
int l = 0, r = n;
while(r - l >1){
int mid = l + (r - l)/2;
int c = query(tl-1, tr, sorted[mid], 0, 0 ,n);
if(c <= k - 1 ) l = mid;
else r = mid;
}
printf("%d\n", sorted[l]);
}
return 0;
}//2000+ms
  • 数组实现的要比vector快很多。
  • 归并树需要二分求解,但是划分树并不需要。因为划分树是从上到下,每次都用数组记录划分到左子树的元素个数,所以可以直接求得区间第k大数,而归并树是由下到上,每次对子树进行简单的合并和排序,并没有对划分到左子树的元素进行追踪,所以需要二分搜索答案,即线段树+二分。所以在求静态区间第k大时划分树也就比归并树要快。

静态区间第k大(归并树)的更多相关文章

  1. poj2104&&poj2761 (主席树&&划分树)主席树静态区间第k大模板

    K-th Number Time Limit: 20000MS   Memory Limit: 65536K Total Submissions: 43315   Accepted: 14296 Ca ...

  2. 主席树(静态区间第k大)

    前言 如果要求一些数中的第k大值,怎么做? 可以先就这些数离散化,用线段树记录每个数字出现了多少次. ... 那么考虑用类似的方法来求静态区间第k大. 原理 假设现在要有一些数 我们可以对于每个数都建 ...

  3. 可持久化线段树(主席树)——静态区间第k大

    主席树基本操作:静态区间第k大 #include<bits/stdc++.h> using namespace std; typedef long long LL; ,MAXN=2e5+, ...

  4. HDU3473--Minimum Sum(静态区间第k大)

    Minimum Sum Time Limit: 16000/8000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Tota ...

  5. 主席树学习笔记(静态区间第k大)

    题目背景 这是个非常经典的主席树入门题——静态区间第K小 数据已经过加强,请使用主席树.同时请注意常数优化 题目描述 如题,给定N个整数构成的序列,将对于指定的闭区间查询其区间内的第K小值. 输入输出 ...

  6. POJ2104-- K-th Number(主席树静态区间第k大)

    [转载]一篇还算可以的文章,关于可持久化线段树http://finaltheory.info/?p=249 无修改的区间第K大 我们先考虑简化的问题:我们要询问整个区间内的第K大.这样我们对值域建线段 ...

  7. HDU 2665 Kth number(主席树静态区间第K大)题解

    题意:问你区间第k大是谁 思路:主席树就是可持久化线段树,他是由多个历史版本的权值线段树(不是普通线段树)组成的. 具体可以看q学姐的B站视频 代码: #include<cmath> #i ...

  8. 主席树初步学习笔记(可持久化数组?静态区间第k大?)

    我接触 OI也快1年了,然而只写了3篇博客...(而且还是从DP跳到了主席树),不知道我这个机房吊车尾什么时候才能摸到大佬们的脚后跟orz... 前言:主席树这个东西,可以说是一种非常畸形的数据结构( ...

  9. POJ 2104 && POJ 2761 (静态区间第k大,主席树)

    查询区间第K大,而且没有修改. 使用划分树是可以做的. 作为主席树的入门题,感觉太神奇了,Orz /* *********************************************** ...

随机推荐

  1. node入门(一)——安装

    node可以让我们用js写服务器.此外还可以用来前端自动化开发,它找到特定服务要使用的包,然后下载.安装.管理. 首先安装node,进入官网下载需要的node版本,然后一键式傻瓜安装.(我的环境是wi ...

  2. php debug/phpstorm调试

    apache+phpstorm调试php代码,修改php.ini配置文件开启调试,没有以下代码加上即可, [XDebug]zend_extension="C:\php\php-7.0.12- ...

  3. Sqlmap脱库之“你的数据我所见”

    废话不多说,介绍和理论啥的网上一搜一大把,一次只为tuoku. 实战操作: 1.在网上找个index.php?id=1的站点 2.很明显看到了脱下裤子的希望,sqlmap走一把.查看数据库 3.怎么看 ...

  4. Android EditText 输入金额(小数点后两位)

    EditText edit = new EditText(context); InputType.TYPE_NUMBER_FLAG_DECIMAL //小数点型 InputType.TYPE_CLAS ...

  5. XCode的debug断点调试

    debug 流程控制 当你通过 Xcode 的源码编辑器的侧边槽 (或者通过下面的方法) 插入一个断点,程序到达断点时会就会停止运行. 调试条上会出现四个你可以用来控制程序的执行流程的按钮. 从左到右 ...

  6. Vue.js——打包之后资源路径产生问题

    https://blog.csdn.net/qq_30632003/article/details/79353035 https://www.cnblogs.com/diantao/p/7776523 ...

  7. 8086汇编——Introduction(8086内部寄存器,段寄存器,存储器分段)

    8086汇编--Introduction 一.8086CPU的三种工作模式 实模式:只有低20位地址线起作用,仅能寻址第一个1MB的内存空间.MS DOS运行于该模式下. 保护模式:在该模式下,机器可 ...

  8. PHP面向对象考察点

    面向对象三大特性 封装 封装性就是把对象的属性和服务结合成一个独立的相同单位,并尽可能隐蔽对象的内部细节,包含两个含义: 把对象的全部属性和全部服务结合在一起,形成一个不可分割的独立单位(即对象). ...

  9. Android(java)学习笔记199:JNI之JNI概念

    1. JNI是什么? java native interface (java本机接口) 比如方法声明: public final native Class<?>  getClass(): ...

  10. Swift 命名空间形式扩展的实现

    Swift 的 extension 机制很强大,不仅可以针对自定义的类型,还能作用于系统库的类型,甚至基础类型比如 Int.当在对系统库做 extension 的时候,就会涉及到一个命名冲突的问题.O ...