刚开始的朴素想法,也就是老师上课说的部分快排。
但做了以后发现并没那么简单,首先是qsort中,以前写的都有等号,但不确定到底是哪个位置
于是改为无等号。
笨办法,部分需排的由a放入b
#include <stdio.h>
int a[200010];
int b[100000];
int knumber;
void qqsort(int low,int high)
{
int l=low,h=high;
b[0]=b[l];
while(l<h)
{
while( l<h && b[h]>=b[0] )h--;
b[l]=b[h];
while( l<h && b[l]<=b[0] )l++;
b[h]=b[l];
}
b[l]=b[0];
if(l==knumber)
return ;
else if(l>knumber)
qqsort(low ,h-1);
else
qqsort(l+1,high);
}
int main()
{
int i,j,n,m,num;
scanf("%d%d",&n,&m);
for(i=1;i<=n;i++)
scanf("%d",&a[i]);
while(m--)
{
scanf("%d%d%d",&i,&j,&knumber);
for(num=1;num+i-1<=j;num++)
b[num]=a[i+num-1];
qqsort(1,num-1);
printf("%d\n",b[knumber]);
}
return 0;
}
因此需要划分树。
参考:
http://yueyue1105.blog.163.com/blog/static/431117682010716111425892/
http://www.cnblogs.com/pony1993/archive/2012/07/17/2594544.html
划分树是一种基于线段树的数据结构。主要用于快速求出(在log(n)的时间复杂度内)序列区间的第k大值 。
划分树和归并树都是用线段树作为辅助的,原理是基于快排 和归并排序 的。
划分树的建树过程基本就是模拟快排过程,取一个已经排过序的区间中值,然后把小于中值的点放左边,大于的放右边。并且记录d层第i个数之前(包括i)小于中值的放在左边的数。
假设要在区间[l,r]中查找第k大元素,t为当前节点,lch,rch为左右孩子,left,mid为节点t左边界和中间点。
1、sum[r]-sum[l-1]>=k,查找lch[t],区间对应为[ left+sum[l-1] , left+sum[r]-1 ]
2、sum[r]-sum[l-1]<k,查找rch[t],区间对应为[ mid+1+l-left-sum[l-1] , mid+1+r-left-sum[r] ]
模版:
/*
HDU 2665 Kth number
划分树
*/
#include<stdio.h>
#include<iostream>
#include<string.h>
#include<algorithm>
using namespace std;
const int MAXN=100010;
int tree[30][MAXN];//表示每层每个位置的值
int sorted[MAXN];//已经排序的数
int toleft[30][MAXN];//toleft[p][i]表示第i层从1到i有多少个数分入左边
void build(int l,int r,int dep)
{
if(l==r)return;
int mid=(l+r)>>1;
int same=mid-l+1;//表示等于中间值而且被分入左边的个数
for(int i=l;i<=r;i++)
if(tree[dep][i]<sorted[mid])//先假设左边的(mid - l + 1)个数都等于as[mid],然后把实际上小于as[mid]的减去
same--;//为了保持树的平衡,不然若是很多相同的则会全部偏向一边
int lpos=l;
int rpos=mid+1;
for(int i=l;i<=r;i++)
{
if(tree[dep][i]<sorted[mid])//比中间的数小,分入左边
tree[dep+1][lpos++]=tree[dep][i];
else if(tree[dep][i]==sorted[mid]&&same>0)
{
tree[dep+1][lpos++]=tree[dep][i];
same--;
}
else //比中间值大分入右边
tree[dep+1][rpos++]=tree[dep][i];
toleft[dep][i]=toleft[dep][l-1]+lpos-l;//从1到i放左边的个数
}
build(l,mid,dep+1);
build(mid+1,r,dep+1);
}
//查询区间第k大的数,[L,R]是大区间,[l,r]是要查询的小区间
//假设要在区间[l,r]中查找第k大元素,t为当前节点,lch,rch为左右孩子,left,mid为节点t左边界和中间点
//sum[r]-sum[l-1]>=k,查找lch[t],区间对应为[ left+sum[l-1], left+sum[r]-1 ]
//sum[r]-sum[l-1]<k,查找rch[t],区间对应为[ mid+1+l-left-sum[l-1], mid+1+r-left-sum[r] ]
int query(int L,int R,int l,int r,int dep,int k)
{
if(l==r)return tree[dep][l];
int mid=(L+R)>>1;
int cnt=toleft[dep][r]-toleft[dep][l-1];//[l,r]中位于左边的个数
if(cnt>=k)
{
//L+要查询的区间前被放在左边的个数
int newl=L+toleft[dep][l-1]-toleft[dep][L-1];
//左端点加上查询区间会被放在左边的个数
int newr=newl+cnt-1;
return query(L,mid,newl,newr,dep+1,k);
}
else
{
int newr=r+toleft[dep][R]-toleft[dep][r];
int newl=newr-(r-l-cnt);
return query(mid+1,R,newl,newr,dep+1,k-cnt);
}
}
int main()
{
//freopen("in.txt","r",stdin);
//freopen("out.txt","w",stdout);
int T;
int n,m;
int s,t,k;
scanf("%d",&T);
while(T--)
{
scanf("%d%d",&n,&m);
memset(tree,0,sizeof(tree));//这个必须
for(int i=1;i<=n;i++)//从1开始
{
scanf("%d",&tree[0][i]);
sorted[i]=tree[0][i];
}
sort(sorted+1,sorted+n+1);
build(1,n,0);
while(m--)
{
scanf("%d%d%d",&s,&t,&k);
printf("%d\n",query(1,n,s,t,0,k));
}
}
return 0;
}