poj2104K-th Number【划分树入门题】

卡了半天攻下的知识点，没有想象中的难啊，某大牛不是说过“主席树是当时我弱不会划分树研究出来的替代品，被一小撮别有用心的人取了这么一个名字”(⊙﹏⊙)b 分明是主席树更难一点好么 ……

说说划分树是个什么玩意，归并排序知道吧，不知道就翻一下大话数据结构，两两排序然后合并，再合并大区间，再合并大大区间……划分树建树的思路和他的差不多，我记得归并排序最开始分几波是任意的，而且这个数字也影响效率，但是划分树拿来一个数组就分两组，找到中间值，小于中间值的放左边，大于的放右边，等于的看数量决定分到哪里，反正肯定是左子树分出去一半，右子树是另一半。刚刚说的这个过程是第二层树。原数组是第一层树。拿出第二层的左、右子树进行递归操作直到子树只有一个数。

建树的时候当然是要给赋值的，要不然建树干嘛。每一层树是一个结构体，结构体里面有三个数组，大小就都是原始数组的长度。三个数组意义分别是：val[i]当前下标的值

sum[i]  记录当前区间内比当前元素小的和（这个有没有看情况，要是只问了第k小的元素是什么，这项就不用要了）

num[i]此下标元素所在区间的当前位置之前进入左子树的个数（什么叫做当前区间？这个问题困扰过我，是指那些分成的当前层的分出的子树，而不是1-n）  

/***************
poj2104
2016.1.28
14768K	1219MS	G++	2059B
***************/
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

const int MAXN = 1000010;
int sorted[MAXN];   //对原集合中元素排序后的值
int val[30][MAXN];  //val记录第k层当前位置的值
int toleft[30][MAXN];  //记录元素所在区间当前位置前的元素进入到左子树的个数
int sum[30][MAXN];  //记录比当前元素小的元素的和
int n;

void build(int l, int r, int d) {
    if (l == r) return ;
    int mid = (l + r) >> 1;
    int same = mid - l + 1;
    for (int i=l; i<=r; i++)
        if (val[d][i] < sorted[mid])
            same--;
    int lp = l, rp = mid+1;
    for (int i=l; i<=r; i++) {
        if (i == l) toleft[d][i] = 0;
        else toleft[d][i] = toleft[d][i-1];

        if (val[d][i] < sorted[mid]) {
            toleft[d][i]++;
            val[d+1][lp++] = val[d][i];
        } else if (val[d][i] > sorted[mid])
            val[d+1][rp++] = val[d][i];
        else {
            if (same) {
                same--;
                toleft[d][i]++;
                val[d+1][lp++] = val[d][i];
            } else val[d+1][rp++] = val[d][i];
        }
    }
    build(l, mid, d+1);
    build(mid+1, r, d+1);
}
int query(int a, int b, int k, int l, int r, int d) {
    if (a == b) return val[d][a];

    int mid = (l + r) >> 1;
    int s, ss, sss;
    if (a == l) {
        s = toleft[d][b];
        ss = 0;
    } else {
        s = toleft[d][b] - toleft[d][a-1];
        ss = toleft[d][a-1];
    }

    if (s >= k) {
        a = l + ss;
        b = l + ss + s - 1;
        return query(a, b, k, l, mid, d+1);
    } else {
        a = mid+1 + a - l - ss;
        b = mid+1 + b - l - toleft[d][b];
        return query(a, b, k-s, mid+1, r, d+1);
    }
}

int main()
{
    //freopen("cin.txt","r",stdin);
     int n, m;
    scanf("%d %d", &n, &m);
    for (int i=1; i<=n; i++) {
        scanf("%d", &sorted[i]);
        val[0][i] = sorted[i];
    }
    sort(sorted+1, sorted+1+n);
    build(1, n, 0);

    int a, b, k;
    while (m--) {
        scanf("%d%d%d", &a, &b, &k);
        printf("%d\n", query(a, b, k, 1, n, 0));
    }
    return 0;
    return 0;
}

对面坐着三个搞数学建模的，实在好奇，刚刚手贱百度了一下，在数学建模吧里有人问ACM和它的取舍，下面的回复除了中立就是推荐前者……加油吧少女~~（捂脸，还有就是今天突发奇想翻了墙，第一次成功登陆CF，但是TC虽然登陆上了，没找到在哪里注册比赛。争取毕业前BC和CF双红吧，加油↖(^ω^)↗