裸的区间第k大问题,划分树搞起。

#pragma comment(linker, "/STACK:10240000")
#include <map>
#include <set>
#include <cmath>
#include <ctime>
#include <deque>
#include <queue>
#include <stack>
#include <vector>
#include <cstdio>
#include <string>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>using namespace std;#define X                   first
#define Y                   second
#define pb                  push_back
#define mp                  make_pair
#define all(a)              (a).begin(), (a).end()
#define fillchar(a, x)      memset(a, x, sizeof(a))
#define fillarray(a, b)     memcpy(a, b, sizeof(a))typedef long long ll;
typedef pair<int, int> pii;
typedef unsigned long long ull;#ifndef ONLINE_JUDGE
void RI(vector<int>&a,int n){a.resize(n);for(int i=0;i<n;i++)scanf("%d",&a[i]);}
void RI(){}void RI(int&X){scanf("%d",&X);}template<typename...R>
void RI(int&f,R&...r){RI(f);RI(r...);}void RI(int*p,int*q){int d=p<q?1:-1;
while(p!=q){scanf("%d",p);p+=d;}}void print(){cout<<endl;}template<typename T>
void print(const T t){cout<<t<<endl;}template<typename F,typename...R>
void print(const F f,const R...r){cout<<f<<", ";print(r...);}template<typename T>
void print(T*p, T*q){int d=p<q?1:-1;while(p!=q){cout<<*p<<", ";p+=d;}cout<<endl;}
#endif
template<typename T>bool umax(T&a, const T&b){return b<=a?false:(a=b,true);}
template<typename T>bool umin(T&a, const T&b){return b>=a?false:(a=b,true);}const double PI = acos(-1.0);
const int INF = 1e9 + 7;
const double EPS = 1e-12;/* -------------------------------------------------------------------------------- */const int maxn = 1e5 + 7;/** 过程:快排的过程,通过记录进入左子区间的个数的前缀和来解决区间第k大问题 **/
class PartitionTree {int cnt[20][maxn], val[20][maxn], buf[maxn];int n;void init(int a[], int n) {this->n = n;fillchar(cnt, 0);fillchar(val, 0);fillarray(val[0], a);fillarray(buf, a);sort(buf, buf + n);}void build(int l, int r, int dep) {if (l == r) return ;int m = (l + r) >> 1, c = 0, small = 0;for (int i = l; i <= r; i ++) small += val[dep][i] < buf[m];for (int i = l; i <= r; i ++) {if (c < m - l + 1) {if (val[dep][i] < buf[m] || val[dep][i] == buf[m] && small < m - l + 1) {cnt[dep][i] = 1;val[dep + 1][l + c ++] = val[dep][i];small += val[dep][i] == buf[m];}}else break;}for (int i = l; i <= r; i ++) {if (!cnt[dep][i]) val[dep + 1][l + c ++] = val[dep][i];}build(l, m, dep + 1);build(m + 1, r, dep + 1);}/** 第k小 */int querykth(int L, int R, int k, int l, int r, int dep) {if (k <= 0 || k > R - L + 1) return - 1;if (L == R) return val[dep][L];int m = (l + r) >> 1, cl = cnt[dep][L - 1] - cnt[dep][l - 1], cr = cnt[dep][R] - cnt[dep][l - 1];if (cr - cl >= k) return querykth(l + cl, l + cr - 1, k, l, m, dep + 1);return querykth(m + 1 + L - l - cl, m + R - l + 1 - cr, k - cr + cl, m + 1, r, dep + 1);}
public:void build(int a[], int n) {init(a, n);build(1, n, 0);for (int i = 0; i < 20; i ++) {for (int j = 2; j <= n; j ++) {cnt[i][j] += cnt[i][j - 1];}}}int querykth(int L, int R, int k) { return querykth(L, R, k, 1, n, 0); }
};/** 下标从1开始 */PartitionTree pt;
int a[maxn];int main() {#ifndef ONLINE_JUDGEfreopen("in.txt", "r", stdin);//freopen("out.txt", "w", stdout);
#endif // ONLINE_JUDGEint n, m;while (cin >> n >> m) {for (int i = 1; i <= n; i ++) {scanf("%d", a + i);}pt.build(a, n);int l, r, k;for (int i = 0; i < m; i ++) {scanf("%d%d%d", &l, &r, &k);printf("%d\n", pt.querykth(l, r, k));}}return 0;
}

转载于:https://www.cnblogs.com/jklongint/p/4749192.html

[NBUT 1458 Teemo]区间第k大问题,划分树相关推荐

  1. 区间第K大(划分树)

    题目:http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1175 注意一个很明显的结论:求区间[p,q]的第K大就等于求区间[ ...

  2. P2617 Dynamic Rankings 动态区间第K大【树套树】或【整体二分】

    传送门 动态区间第KKK大问题,单点修改(这里是第kkk小,即是从小到大第kkk个) 这里还有个区间修改,有点类似的 P3332 [ZJOI2013]K大数查询 分析 树套树 树套树,就是用一种树形结 ...

  3. 动态区间第k小:树状数组套权值线段树

    所谓树状数组套权值线段树,就是在树状树组上套权值线段树 (逃) 解析 如何解决静态区间第k小? 使用主席树就ok啦 辣么如何解决动态区间第k小嘞- 我们想想主席树为啥不能解决动态区间第k小 因为如果改 ...

  4. J - [永不止步-2017]_区间第K大-线段树维护

    J - [永不止步-2017]_区间第K大 把线段树的结点的数据域设置为vector类型即可别的操作为区间更新模板 思路就是这样runtime error暂时没改对 #include<bits/ ...

  5. R19436221 区间第k小 主席树

    主席树模板题,记录一下. #include<bits/stdc++.h> using namespace std; const int maxn=2e5+5; int root[maxn] ...

  6. HDU2665 求区间第K大 主席树

    题目链接 http://acm.hdu.edu.cn/showproblem.php?pid=2665 代码: //#include<bits/stdc++.h> #include< ...

  7. hdu4417(Super Mario)—— 二分+划分树

    Super Mario Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Tota ...

  8. POJ 2104 K-th Number(区间第k大数)(平方切割,归并树,划分树)

    题目链接: http://poj.org/problem? id=2104 解题思路: 由于查询的个数m非常大.朴素的求法无法在规定时间内求解. 因此应该选用合理的方式维护数据来做到高效地查询. 假设 ...

  9. 利用划分树求解整数区间内第K大的值

    如何快速求出(在log2n的时间复杂度内)整数区间[x,y]中第k大的值(x<=k<=y)? 其实我刚开始想的是用快排来查找,但是其实这样是不行的,因为会破坏原序列,就算另外一个数组来存储 ...

最新文章

  1. Linux-man命令
  2. bp神经网络训练_数据分析模型6——神经网络基础(人工智能的底层模型)
  3. 中心极限与大数定理律的关系_多元函数的极限、连续性分析
  4. 如何在JavaScript中将字符串的首字母大写?
  5. iOS开发Cocoapods执行命令pod setup,执行失败解决RPC failed; curl 56 LibreSSL SSL_read: SSL_ERROR_SYSCALL, errno 54
  6. 奇迹私服gs服务器端口未能连接请设置参数,奇迹私服架设之各个快捷方式参数...
  7. 树莓派 kali系统默认密码
  8. 设计主导型思维在商业领域的崛起,将对创意职业产生积极的溢出效应
  9. windows删除注册表(通用方法)
  10. Flask框架学习笔记10-jinjia2模板继承实例
  11. 通过mysql修改后台密码_怎么通过修改数据库修改网站后台的管理员密码?
  12. [BZOJ3503][CQOI2014]和谐矩阵
  13. 简单的微信使用技巧,你需要掌握的技巧
  14. win8、server 2012 清除winsxs文件夹
  15. 简单的猜数字游戏[1]
  16. javascript开发HTML5游戏--斗地主(单机模式part3)
  17. Deeplav V3总结
  18. SpringAOP原理+使用
  19. 社交网络分析--python-igraph
  20. Python mac运行pygame一片空白不显示内容-解决方案

热门文章

  1. 【算法小积累】 - 提取非0数最右侧的1
  2. node --- 实践中理解跨域
  3. DOS Network一月项目月报
  4. Python并发编程之concurrent.futures
  5. T24412 Cup#182-3 洞穴之旅
  6. redis-cli使用密码登录
  7. JeeSite 4.0 说说前端的那些事
  8. APNIC IP 库
  9. Android底部导航栏实现(一)之BottomNavigationBar
  10. Markdown文件导出为HTML的小程序