题目链接

Description

You are working for Macrohard company in data structures department. After failing your previous task about key insertion you were asked to write a new data structure that would be able to return quickly k-th order statistics in the array segment.

That is, given an array a[1...n] of different integer numbers, your program must answer a series of questions Q(i, j, k) in the form: "What would be the k-th number in a[i...j] segment, if this segment was sorted?"

For example, consider the array a = (1, 5, 2, 6, 3, 7, 4). Let the question be Q(2, 5, 3). The segment a[2...5] is (5, 2, 6, 3). If we sort this segment, we get (2, 3, 5, 6), the third number is 5, and therefore the answer to the question is 5.

Input

The first line of the input file contains n --- the size of the array, and m --- the number of questions to answer (1 <= n <= 100 000, 1 <= m <= 5 000).

The second line contains n different integer numbers not exceeding 109 by their absolute values --- the array for which the answers should be given.

The following m lines contain question descriptions, each description consists of three numbers: i, j, and k (1 <= i <= j <= n, 1 <= k <= j - i + 1) and represents the question Q(i, j, k).

Output

For each question output the answer to it --- the k-th number in sorted a[i...j] segment.

Sample Input

7 3

1 5 2 6 3 7 4

2 5 3

4 4 1

1 7 3

Sample Output

5

6

3

题意:

给定一个数列a1,a2,a3,···,an和m个三元组表示的查询。对于每个查询(i,j,k),输出ai,ai+1,···,aj的升序排列中第k个数。

分析:

因为查询的个数m很大,朴素的求法没法在规定的时间内求出。

如果x是第k个数,那么一定有:

1.在区间中不超过x的数不少于k个

2.在区间中小于x的数不到k个

因此,如果可以快速求出区间里不超过x的数的个数,就可以通过对x进行二分搜索求出第k个数是多少。

接下来,我们来看一下如何计算在某个区间里不超过x的数的个数。如果不进行预处理,那么就只能遍历一遍所有的元素。

另一方面,如果区间是有序的,那么就可以通过二分搜索法高效地求出不超过x的数的个数了。但是,如果对于每个查询都分别进行一次排序,就完全无法降低复杂度。所以,可以考虑使用平方分割和线段树进行求解。

首先看一下如何使用平方分割来解决这个问题。把数列每b个一组分到哥哥桶里,每个桶内保存有序后的序列。这样,如果要求在某个区间中不超过x的个数,就可以这样求得:

1.对于完全包含在区间内的桶,用二分搜索法进行计算

2.对于所在的桶不完全包含在区间内的元素,逐个检查。

#include<cstdio>
#include<iostream>
#include<string.h>
#include<vector>
#include<algorithm>
using namespace std;
typedef long long ll;
const int B = 1072;
const int maxn = 1e5+5;
int n, m;
int a[maxn], num[maxn];///原序列和排序后的序列
vector<int> buk[maxn/B+1];///每个桶排序后的结果 int main()
{ scanf("%d%d",&n,&m);
for(int i = 0; i < n; i++)
{
scanf("%d", a+i);
buk[i/B].push_back(a[i]);
}
memcpy(num,a,sizeof(int)*n);
sort(num,num+n);
///(n-1)/b <= n/b 当 n%b != 0时候取等号,这时候不是一个完整的桶,所以最后一个桶不sort也没关系
for(int i = n/B;i>=0; i--) sort(buk[i].begin(),buk[i].end());
while(m--)
{
int l,r,k;
scanf("%d%d%d",&l,&r,&k);
int lb = 0, ub = n-1;
while(lb < ub)
{
int md = (lb+ub)>>1;
int x = num[md];
///求[l,r)区间中第k个数
int p = l-1, q = r, c = 0;
///区间两端多出的部分
while(p<q && p%B) if(a[p++] <= x) c++;
while(p<q && q%B) if(a[--q] <= x) c++;
///对每一个桶进行计算
for(int i = p/B, mi = q/B; i < mi; i++)
{
c += upper_bound(buk[i].begin(),buk[i].end(),x)-buk[i].begin();
}
c >= k?ub = md: lb = md+1;
}
printf("%d\n", num[lb]);
} return 0;
}

思路同上,只是计算cnt的时候改用线段树。线段树每个结点保存有序的数组,建树过程是归并排序的完整再现,因此也叫归并树。

更抽象来看,数值大小和位置是两个维度,查询就是询问一个空间内的点数,也就是所谓的区域树了,通过嵌套可以推广到高维。

#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<queue>
#include<vector>
#include<stack>
#include<vector>
#include<map>
#include<set>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;
#define para int d = 0, int l = 0,int r = n
#define lsn d+1, l, mid
#define rsn d+1, mid, r
#define insd ql<=l&&r<=qr
const int maxn = 1e5;
int n, m;
int a[maxn];
const int Log2N_ = 18; int dat[Log2N_][maxn]; void build(para)
{
if(r-l == 1)
{
dat[d][l] = a[l];
}
else
{
int mid = (l+r)>>1;
build(lsn);
build(rsn);
merge(dat[d+1]+l,dat[d+1]+mid,dat[d+1]+mid,dat[d+1]+r,dat[d]+l);
}
} int ql, qr, qval;
int query(para)
{
if(insd)
{
return upper_bound(dat[d]+l,dat[d]+r,qval) - dat[d]-l;
}
else
{
int res = 0;
int mid = (l+r)>>1;
if(ql<mid) res += query(lsn);
if(qr>mid) res += query(rsn);
return res;
}
} int main()
{
scanf("%d%d",&n,&m);
for(int i = 0; i < n; i++)
{
scanf("%d", a+i);
}
build();
while(m--)
{
int k;
scanf("%d%d%d",&ql,&qr,&k);
ql--;
int lb = 0, ub = n-1;
while(lb < ub)
{
int md = (lb+ub)>>1;
qval = dat[0][md];
query() >= k? ub = md:lb = md+1;
}
printf("%d\n", dat[0][lb]);
} return 0;
}

POj 2104 K-th Number (分桶法+线段树)的更多相关文章

  1. POJ 2104 K-th Number(分桶,线段树,主席树)

    一道比较经典的数据结构题.可以用多种方式来做. 一,分桶法(平方分解). 根据数字x的大小和区间内不大于x的数字数量cnt的单调性,可知第k大数kth对应的cnt应该满足cnt≥k, 且kth是满足条 ...

  2. 静态区间第k大(分桶法和平方分割)

    POJ 2104为例 思想: <挑战程序设计竞赛>中介绍的方法. 分桶法:把一排物品或者平面分成桶,每个桶分别维护自己内部的信息,已达到高效计算的目的. 设一共有n个数,每b个分到一个桶里 ...

  3. LeetCode 220. Contains Duplicate III (分桶法)

    Given an array of integers, find out whether there are two distinct indices i and j in the array suc ...

  4. 【POJ 2482】 Stars in Your Window(线段树+离散化+扫描线)

    [POJ 2482] Stars in Your Window(线段树+离散化+扫描线) Time Limit: 1000MS   Memory Limit: 65536K Total Submiss ...

  5. poj 3468:A Simple Problem with Integers(线段树,区间修改求和)

    A Simple Problem with Integers Time Limit: 5000MS   Memory Limit: 131072K Total Submissions: 58269   ...

  6. POJ.3468 A Simple Problem with Integers(线段树 区间更新 区间查询)

    POJ.3468 A Simple Problem with Integers(线段树 区间更新 区间查询) 题意分析 注意一下懒惰标记,数据部分和更新时的数字都要是long long ,别的没什么大 ...

  7. POJ 2104 K-th Number (可持久化线段树)

    题目大意 给一个长度为n的序列,有m个询问,每次询问一个区间里面第k小的数. 解题分析 静态的区间第k大.复习了一下可持久化线段树. 首先对数值离散化,建一颗权值线段树.按照序列的顺序依次插入,每一个 ...

  8. poj 3468 A Simple Problem with Integers(线段树、延迟更新)

    A Simple Problem with Integers Time Limit: 5000MS   Memory Limit: 131072K Total Submissions: 74705   ...

  9. poj 3468 A Simple Problem with Integers(线段树+区间更新+区间求和)

    题目链接:id=3468http://">http://poj.org/problem? id=3468 A Simple Problem with Integers Time Lim ...

随机推荐

  1. Swift-重写(Override)

    子类可以为继承来的实例方法(instance method),类方法(class method),实例属性(instance property),或附属脚本(subscript)提供自己定制的实现(i ...

  2. SQL Server之看懂执行计划

    在SQL Server中,选中一段SQL按Ctrl+L,就可以查看它的执行计划. 上面是一个执行计划的实例,在SQL Server中,执行计划是从右往左看的. SQL Server中,查找数据的方式有 ...

  3. DDL、DML和DCL的比较【引用学习】

    1.DDL       1-1.DDL的概述                DDL(Data Definition Language 数据定义语言)用于操作对象和对象的属性,这种对象包括数据库本身,以 ...

  4. Android基础------高级ul:消息提示

    前言:Android消息提示笔记,刚刚接触Android 1.静态方法Toast 直接调用静态方法 //消息提示(context,"内容",固定时间) Toast.makeText ...

  5. mysql表、函数等被锁住无响应的问题

    场景: 在对表或函数等进行操作的时候,如果出现无法响应的情况(排除外网的网络问题),此时极有可能被某一个线程锁定了(这是函数的情况,表的话可能是被某一个用户锁定了),锁定的原因一般都是死循环出不来,而 ...

  6. Eclipse闪退解决方案

    1. 找到Eclipse目录下的eclipse.exe,右键点击->发送到桌面快捷方式,然后右键点击快捷方式,选择属性,修改“目标”(或target),其中红色代表eclipse的路径,绿色代表 ...

  7. Socket_SSH-3(粘包)

    粘包:两次数据粘到一起了.在Windows中基本看不出来效果. 服务器端的配置: import socket,os,time server=socket.socket() server.bind((' ...

  8. POJ3498:March of the Penguins——题解

    最近的题解的故事背景割. 题目: 描述 在靠近南极的某处,一些企鹅站在许多漂浮的冰块上.由于企鹅是群居动物,所以它们想要聚集到一起,在同一个冰块上.企鹅们不想把自己的身体弄湿,所以它们在冰块之间跳跃, ...

  9. 51NOD 1773:A国的贸易——题解

    http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1773 参考1:FWT讲解 https://www.cnblogs.com ...

  10. python获取当前工作目录

    py文件所在位置/test/pj/hello.py 用户所在位置:/ 用户执行命令python /test/pj/hello.py 1. os.getcwd() 返回的是执行命令的位置 / 2.sys ...