K-th Number (划分树学习)

刚开始的朴素想法，也就是老师上课说的部分快排。

但做了以后发现并没那么简单，首先是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;
}