HDU 4417 Super Mario（划分树）#by zh

比赛的时候看到这题有很多的查询操作，感觉就是道线段树题，果断用线段树写了，然后悲剧的TLE了……之后问过了的人才知道这题要用划分树写，然后二分第K大数才能过，不过我当时根本不知道划分树是神马，回来之后在网上找了个模版，套模版才A掉了这题，不过现在还是不明白这个数据结构的具体实现，晚上要研究一下

#include <iostream>
#include <cstdio>
#include <string.h>
#include <algorithm>
using namespace std;
#define M 100001
#define LL(x) (x<<1)
#define RR(x) (x<<1|1)
struct Seg_Tree
{
    int left,right;
    int mid()
    {
        return (left + right) >> 1;
    }
} tt[M<<2];
int len;
int sorted[M];
int toLeft[20][M];
int val[20][M];

void build(int l,int r,int d,int idx)
{
    tt[idx].left = l;
    tt[idx].right = r;
    if(tt[idx].left == tt[idx].right)	return ;
    int mid = tt[idx].mid();
    int lsame = mid - l + 1;//lsame表示和val_mid相等且分到左边的
    for(int i = l ; i <= r ; i ++)
    {
        if(val[d][i] < sorted[mid])
        {
            lsame --;//先假设左边的数(mid - l + 1)个都等于val_mid,然后把实际上小于val_mid的减去
        }
    }
    int lpos = l;
    int rpos = mid+1;
    int same = 0;
    for(int i = l ; i <= r ; i ++)
    {
        if(i == l)
        {
            toLeft[d][i] = 0;//toLeft[i]表示[ tt[idx].left , i ]区域里有多少个数分到左边
        }
        else
        {
            toLeft[d][i] = toLeft[d][i-1];
        }
        if(val[d][i] < sorted[mid])
        {
            toLeft[d][i] ++;
            val[d+1][lpos++] = val[d][i];
        }
        else if(val[d][i] > sorted[mid])
        {
            val[d+1][rpos++] = val[d][i];
        }
        else
        {
            if(same < lsame)  //有lsame的数是分到左边的
            {
                same ++;
                toLeft[d][i] ++;
                val[d+1][lpos++] = val[d][i];
            }
            else
            {
                val[d+1][rpos++] = val[d][i];
            }
        }
    }
    build(l,mid,d+1,LL(idx));
    build(mid+1,r,d+1,RR(idx));
}

int query(int l,int r,int k,int d,int idx)
{
    if(l == r)
    {
        return val[d][l];
    }
    int s;//s表示[ l , r ]有多少个分到左边
    int ss;//ss表示 [tt[idx].left , l-1 ]有多少个分到左边
    if(l == tt[idx].left)
    {
        s = toLeft[d][r];
        ss = 0;
    }
    else
    {
        s = toLeft[d][r] - toLeft[d][l-1];
        ss = toLeft[d][l-1];
    }
    if(s >= k)  //有多于k个分到左边,显然去左儿子区间找第k个
    {
        int newl = tt[idx].left + ss;
        int newr = tt[idx].left + ss + s - 1;//计算出新的映射区间
        return query(newl,newr,k,d+1,LL(idx));
    }
    else
    {
        int mid = tt[idx].mid();
        int bb = l - tt[idx].left - ss;//bb表示 [tt[idx].left , l-1 ]有多少个分到右边
        int b = r - l + 1 - s;//b表示 [l , r]有多少个分到右边
        int newl = mid + bb + 1;
        int newr = mid + bb + b;
        return query(newl,newr,k-s,d+1,RR(idx));
    }
}
int main()
{
    //freopen("input.txt","r",stdin);
    int n,m,t;
    scanf("%d",&t);
    for(int cas=1; cas<=t; cas++)
    {
        scanf("%d %d",&n,&m);
        printf("Case %d:\n",cas);
        for(int i = 1 ; i <= n ;  i++)
        {
            scanf("%d",&val[0][i]);
            sorted[i] = val[0][i];
        }
        sort(sorted+1,sorted+n+1);
        build(1,n,0,1);
        while(m--)
        {
            int a,b,h;
            scanf("%d %d %d",&a,&b,&h);
            int low=1,high=b-a+1,mid;
            while(low<=high)
            {
                mid=(high+low)>>1;
                if(query(a+1,b+1,mid,0,1)<=h)
                    low=mid+1;
                else high=mid-1;
            }
            printf("%d\n",low-1);
        }
    }
    return 0;
}