hdu 4417(2012 杭州网络赛)

唉，比赛的时候怎么都不会啊！！

题解：http://blog.csdn.net/acm_cxlove/article/details/8014170

划分树+二分答案：划分树可以方便的求解k-number。再利用二分答案，即区间内小于h的个数（最大为r-l+1,最小为0）。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;

#define ls rt<<1
#define rs rt<<1|1
#define lson l,m,ls
#define rson m+1,r,rs

#define MAXN 100010

struct Seg_Tree
{
    int l,r;
}tt[MAXN<<2];
int len ,sorted[MAXN],toLeft[20][MAXN];
int val[20][MAXN];

void build(int d,int l,int r,int rt)
{
    tt[rt].l=l;
    tt[rt].r=r;
    if(tt[rt].l == tt[rt].r) return ;
    int m=(l+r)>>1;
    int lsame = m-l+1;
    for(int i=l;i<=r;i++)
        if(val[d][i] < sorted[m])
            lsame--;
    int lpos =l,rpos=m+1,same=0;
    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[m])
        {
            toLeft[d][i]++;
            val[d+1][lpos++] = val[d][i];
        }
        else if(val[d][i] > sorted[m])
            val[d+1][rpos++] = val[d][i];
        else
        {
            if(same<lsame)
            {
                same++;
                toLeft[d][i]++;
                val[d+1][lpos++] = val[d][i];
            }
            else
                val[d+1][rpos++] = val[d][i];
        }
    }
    build(d+1,lson);
    build(d+1,rson);
}

int query(int k,int d,int l,int r,int rt)
{
    if(l==r) return val[d][l];
    int s,ss;
    if(l == tt[rt].l)
    {
        s = toLeft[d][r];
        ss=0;
    }
    else
    {
        s=toLeft[d][r] - toLeft[d][l-1];
        ss=toLeft[d][l-1];
    }
    if(s >= k)
    {
        int newl = tt[rt].l+ss;
        int newr = tt[rt].l+ss+s-1;
        return query(k,d+1,newl,newr,ls);
    }
    else
    {
        int m = (tt[rt].l + tt[rt].r)>>1;
        int bb = l-tt[rt].l -ss;
        int b= r-l+1-s;
        int newl = m+bb+1;
        int newr = m+bb +b;
        return query(k-s,d+1,newl,newr,rs);
    }
}

int bsearch(int r,int h,int L,int R)//r代表最多有多少个
{
    int l =0;
    while(l<r)
    {
        int m = (l+r)>>1;//二分答案（个数）
        if(query(m,0,L,R,1) >h)
            r=m;
        else
            l=m+1;
    }
    return l;
}

int main()
{
    int t,n,m;
    scanf("%d",&t);
    int cas=1;
    while(t--)
    {
        scanf("%d%d",&n,&m);
        for(int i=1;i<=n;i++)
        {
            scanf("%d",&val[0][i]);
            sorted[i]=val[0][i];
        }
        sort(sorted+1,sorted+n+1);
        build(0,1,n,1);
        printf("Case %d:\n",cas++);
        while(m--)
        {
            int l,r,h;
            scanf("%d%d%d",&l,&r,&h);
            if(query(1,0,l+1,r+1,1) > h)
                puts("0");
            else if(query(r-l+1,0,l+1,r+1,1) <= h)
                printf("%d\n",r-l+1);
            else
            {
                int res = bsearch(r-l+1,h,l+1,r+1);
                while(res != 0 && query(res,0,l+1,r+1,1) >h)
                    res--;
                printf("%d\n",res);
            }
        }
    }
    return 0;
}

原来是将所有的询问读入，按H从小到大排序。然后对于所有的结点也按从小到大排序，然后根据查询的H，将比H小的点加入到线段树，然后就是一个很简单的区间和，，，，，做的题还是太少了。

后来写的离线线段树：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>

using namespace std;

#define mid (l+r)>>1
#define ls rt<<1
#define rs rt<<1|1
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1

#define MAXN 100010

int sum[MAXN<<2],ans[MAXN];

struct Que
{
    int l,r,h,id;
}Q[MAXN];
struct node
{
    int h,id;
}in[MAXN];
int n,m;
bool cmp1(const Que &a,const Que &b)
{
    return a.h<b.h;
}
bool cmp2(const node & a,const node &b)
{
    return a.h<b.h;
}
void pushup(int rt)
{
    sum[rt]=sum[ls]+sum[rs];
}

void build(int l,int r,int rt)
{
    sum[rt]=0;
    if(l==r) return;
    int m = mid;
    build(lson);
    build(rson);
}
void update(int p,int q,int l,int r,int rt)
{
    if(l==r)
    {
        sum[rt]+=q;
        return ;
    }
    int m=mid;
    if(p<=m) update(p,q,lson);
    else update(p,q,rson);
    pushup(rt);
}
int query(int L,int R,int l,int r,int rt)
{
    if(L<=l && r<=R)
        return sum[rt];
    int m=mid;
    int ret=0;
    if(L<=m) ret+=query(L,R,lson);
    if(m<R) ret+=query(L,R,rson);
    return ret;
}

int main()
{
    int t;
    scanf("%d",&t);
    int cas=1;
    while(t--)
    {
        scanf("%d%d",&n,&m);
        for(int i=0;i<n;i++)
        {
            scanf("%d",&in[i].h);
            in[i].id=i+1;
        }
        for(int i=0;i<m;i++)
        {
            scanf("%d%d%d",&Q[i].l,&Q[i].r,&Q[i].h);
            Q[i].id=i+1;
        }
        build(1,MAXN,1);
        sort(Q,Q+m,cmp1);
        sort(in,in+n,cmp2);

        int i=0,j=0;
        while(i<m)
        {
            while(in[j].h<=Q[i].h)
            {
                update(in[j].id,1,1,MAXN,1);
                j++;
            }
            ans[Q[i].id]=query(Q[i].l+1,Q[i].r+1,1,MAXN,1);
            i++;
        }
        printf("Case %d:\n",cas++);
        for(int i=1;i<=m;i++)
            printf("%d\n",ans[i]);

    }
    return 0;
}