很久没打代码了,不知道为什么,昨天考岭南文化之前突然开始思考起这个问题来,这个问题据说有很多种方法,划分树什么的,不过对于我现在这种水平还是用熟悉的线段树做比较好。这到题今年8月份的时候曾经做过,那个时候是作为对函数式线段树的一个基础题来做的,那个时候不懂,看着别人的代码对拍然后摸索下来,所以到了昨天我就已经彻底忘了,因为代码不是自己写的。昨天想了很久,终于参透了它的精髓。

首先对于给出的a[1]~a[n]这么多个数离散化,然后建立一个线段树,线段树中的结点对应的区间[L,R]表示的是在[L,R](离散化后的值)中有多少个数。

举个例子,a={1,3,2,6,4,7}  那么假如我对这整个数组建线段树,[3,6]的sum应该为2,因为[3,6]里只有3和5在数组里。

函数式的思想是,假如我建出来从a[1]~a[i]对应的线段树T[i],那么当我要去求数组中的某一段[x,y]在[L,R]里有多少个的时候,实际上就是T[y].[L,R]的和-T[x-1].[L,R]的和,这样我们就可以知道,在数组的a[x]~a[y]里在范围在L-R中的数有多少个。这样的话我们可以二分出第kth值。

即从0~size,如果0~mid里已经有超过k个了,说明第k大在左边,否则第k大在右边,这样查找就可以了。这次我终于写出了一个自己的二分。以前二分都是lower_bound水过,如果要用到的话,我会在l+1=r的时候特判,特别猥琐,看了一篇关于二分的总结才知道,二分的写法有分左闭右闭和左开右开的,不过这些都要慢慢领悟了。

再谈谈怎么建线段树,假如我们已经建了线段树T[i],那么当我们要建线段树T[i+1]的时候,如果再建一次空间会非常大,所以要充分利用T[i]的信息,实际上建立T[i+1]的时候只比T[i]多了一个值,所以我们只需要延着路径重建这些新的结点就好,另外一边的只需要用回历史版本的就好。延着线段树往下会修改到logn个结点,所以建立T[0]时需要开4*n,再加上后来每个一个都要logn个结点,所以总的结点数是(4+logn)*n.  logn=16.多,这也就是为什么我们的结点数要开到20倍,所以这种树对空间的消耗是蛮大的。

想想之前线段树都要套模板,现在可以自己敲出来,这一年里还是有进步的。我还是慢慢努力练习练习吧~

贴一记代码,1700ms,非常慢,权当学习- -0

#include<iostream>
#include<algorithm>
#include<cstring>
#include<string>
#include<cstdio>
#include<cmath>
#define maxn 100000
using namespace std; struct Node
{
int l,r,sum;
Node *lc,*rc;
}N[20*maxn];
int top; Node *root[100050];
int n,m;
int a[maxn+20];
int b[maxn+20];
int bsize; Node *build(int L,int R)
{
if(L==R){
Node *ret=&N[top++];
ret->l=L;ret->r=R;ret->sum=0;
ret->lc=NULL;ret->rc=NULL;
return ret;
}
else{
Node* ret=&N[top++];
int M=(L+R)>>1;
ret->l=L;ret->r=R;
ret->lc=build(L,M);ret->rc=build(M+1,R);
ret->sum=ret->lc->sum+ret->rc->sum;
return ret;
}
} Node *build(int L,int R,int k,Node *p) // 对key值k,和历史版本的p重建树
{
if(L==R){
Node *ret=&N[top++];
ret->l=L;ret->r=R;ret->sum=p->sum+1;
ret->lc=NULL;ret->rc=NULL;
return ret;
}
else{
int M=(L+R)>>1;
Node *ret=&N[top++];
ret->l=L;ret->r=R;
if(k<=M){
ret->lc=build(L,M,k,p->lc); // k值在原树的左边的时候要重建左子树
ret->rc=p->rc; // 右子树不变,用历史版本
ret->sum=ret->lc->sum+ret->rc->sum; // 更新和值
return ret;
}
else{
ret->rc=build(M+1,R,k,p->rc); // k值在原树的右边的时候要重建右子树
ret->lc=p->lc; // 左子树不变,用历史版本
ret->sum=ret->lc->sum+ret->rc->sum; // 更新和值
return ret;
}
}
} int query(Node *p,int L,int R) // 对结点p的[L,R]区间求和值
{
if(p->l==L&&p->r==R){
return p->sum;
}
else{
int M=(p->l+p->r)>>1;
if(R<=M) return query(p->lc,L,R);
else if(L>M) return query(p->rc,L,R);
else return query(p->lc,L,M)+query(p->rc,M+1,R);
}
} int search(int x,int y,int k) //找从x到y上的区间第k大
{
int l=0,r=bsize;int mid,num;
while(l<r){
mid=(l+r)>>1;
num=query(root[y],l,mid)-query(root[x-1],l,mid); // 算出<=mid的有多少个,利用区间可以减的性质
if(k<=num) r=mid;
else{
k-=num;l=mid+1;
}
}
return l;
} int main()
{
while(cin>>n>>m)
{
for(int i=0;i<n;++i){
scanf("%d",&a[i]);
b[i]=a[i];
}
//离散化的过程
sort(b,b+n);bsize=unique(b,b+n)-b;
for(int i=0;i<n;++i){
a[i]=lower_bound(b,b+bsize,a[i])-b;
}
top=0;
root[0]=build(0,bsize);
for(int i=1;i<=n;++i){
root[i]=build(0,bsize,a[i-1],root[i-1]); // 根据a[i-1]的值,以及上一个版本的树建立新的子树
}
int li,ri,ki;
for(int i=0;i<m;++i){
scanf("%d%d%d",&li,&ri,&ki);
printf("%d\n",b[search(li,ri,ki)]);
}
}
return 0;
}

POJ2104 K-th number 函数式线段树的更多相关文章

  1. BZOJ 3123 森林(函数式线段树)

    题目链接:http://61.187.179.132/JudgeOnline/problem.php?id=3123 题意: 思路:总的来说,查询区间第K小利用函数式线段树的减法操作.对于两棵树的合并 ...

  2. 学习笔记--函数式线段树(主席树)(动态维护第K极值(树状数组套主席树))

    函数式线段树..资瓷 区间第K极值查询 似乎不过似乎划分树的效率更优于它,但是如果主席树套树状数组后,可以处理动态的第K极值.即资瓷插入删除,划分树则不同- 那么原理也比较易懂: 建造一棵线段树(权值 ...

  3. POJ- 2104 hdu 2665 (区间第k小 可持久化线段树)

    可持久化线段树 也叫函数式线段树也叫主席树,其主要思想是充分利用历史信息,共用空间 http://blog.sina.com.cn/s/blog_4a0c4e5d0101c8fr.html 这个博客总 ...

  4. BZOJ 3207 花神的嘲讽计划Ⅰ(函数式线段树)

    题目链接:http://61.187.179.132/JudgeOnline/problem.php?id=3207 题意:给出一个数列,若干询问.每个询问查询[L,R]区间内是否存在某个长度为K的子 ...

  5. hdu 5111 树链剖分加函数式线段树

    这题说的是给了两棵树,各有100000 个节点,然后Q个操作Q<=50000; 每个操作L1 R1 L2 R2.因为对于每棵树都有一个与本棵树其他点与众不同的值, 最后问 在树上从L1到R1这条 ...

  6. Codeforces538F A Heap of Heaps(函数式线段树)

    题意:给你一个数组a[n],对于数组每次建立一个完全k叉树,对于每个节点,如果父节点的值比这个节点的值大,那么就是一个违规点,统计出1~n-1完全叉树下的违规点的各自的个数. 一个直觉的思想就是暴力, ...

  7. [Usaco2014 Open Gold ]Cow Optics (树状数组+扫描线/函数式线段树)

    这道题一上手就知道怎么做了= = 直接求出原光路和从目标点出发的光路,求这些光路的交点就行了 然后用树状数组+扫描线或函数式线段树就能过了= = 大量的离散+模拟+二分什么的特别恶心,考试的时候是想到 ...

  8. HDU 1394 Minimum Inversion Number(线段树求最小逆序数对)

    HDU 1394 Minimum Inversion Number(线段树求最小逆序数对) ACM 题目地址:HDU 1394 Minimum Inversion Number 题意:  给一个序列由 ...

  9. BZOJ 3110([Zjoi2013]K大数查询-区间第k大[段修改,在线]-树状数组套函数式线段树)

    3110: [Zjoi2013]K大数查询 Time Limit: 20 Sec   Memory Limit: 512 MB Submit: 418   Solved: 235 [ Submit][ ...

随机推荐

  1. 【学习】Windows PE文件学习(一:导出表)

    今天做了一个读取PE文件导出表的小程序,用来学习. 参考了<Windows PE权威指南>一书. 首先, PE文件的全称是Portable Executable,可移植的可执行的文件,常见 ...

  2. 安装Java EE失败,解决方案

    笔者安装Java EE(版本是java_ee_sdk-7-jdk7-windows-x64-ml.exe)时,遇到错误提示提示"Could not find the required ver ...

  3. HttpWebResponse取不到Cookie?原来是因为被跳转了

    今天做模拟登陆的时候,发现HttpWebResponse的Cookie都为空,但是Fiddler看是有的...后来看见是302状态,才知道请求这个的时候,Response回来已经是跳转了...这样Co ...

  4. html5圆角

    以下是代码片段:#roundCornerI{   font-family: Arial;   border: 5px solid #dedede;   -moz-border-radius-tople ...

  5. 【Qt】QWidget、QDialog、QMainWindow的异同点【转】

    简述 在分享所有基础知识之前,很有必要在这里介绍下常用的窗口-QWidget.QDialog.QMainWindow. 熟悉Qt的同学都应该知道,在新建Qt Widgets项目进行类信息选择时会碰到它 ...

  6. php保存base64数据

    php保存base64数据 if(isset($param['cover_pic']) && !empty($param['cover_pic'])) {
 if (preg_matc ...

  7. c#配置log4net步骤

    1.引入添加log4net.dll引用 2.建立配置文件Log4Net.config(名字自定义).文件内容参考,输出的文件名称可更改 .运行是要放入到相应bin/debug(release) 目录 ...

  8. AJAX请求和跨域请求详解(原生JS、Jquery)

    一.概述 AJAX 是一种在无需重新加载整个网页的情况下,能够更新部分网页的技术. AJAX = 异步 JavaScript 和 XML,是一种用于创建快速动态网页的技术.通过在后台与服务器进行少量数 ...

  9. WPF学习02:Routed Events

    与传统的桌面开发相比,在事件模型上WPF引入了Routed Events,从开发者的角度上,我们获得了两个便利: 1.可以实现事件路由,即向XAML结构中的父元素路由或者是向子元素路由. 2. Rou ...

  10. SQL Server 本地语言版本

    要一些实验是往往喜欢使用英文的Windows 以及SQL Server ,但有时需要使用中文的环境方便理解.中文的SQL Server 不能被安装在英文的Windows 系统上. 根据文档可得知以下兼 ...