SPOJ 27000 BGSHOOT - Shoot and kill（ST）

Description
给出n个猎物的出没时间，Q次查询，每次查询一个时间段[L,R]内可以捕到的猎物数量最大值
Input
第一行一整数n表示猎物数量，之后n行每行两个整数l[i],r[i]表示该猎物的出没时间，然后输入一个整数q表示查询数，每次查询输入一个区间L,R表示查询时间段[L,R]内可以捕到的猎物数最大值(1<=n,q<=1e5,1<=l[i],r[i],L,R<=1e9)
Output
对于每次查询，输出时间段[L,R]内可以捕到的猎物数最大值
Sample Input
4
1 2
2 3
4 5
6 7
4
1 5
2 3
4 7
5 7
Sample Output
2
2
1
1
Solution
用一遍前缀和得到每个时间出现的猎物数，然后用ST表维护一下最大值，对于每次查询O(1)输出结果，注意离散化
Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
#define INF 0x3f3f3f3f
#define maxn 444444
int n,q,x[maxn],y[maxn],l[maxn],r[maxn],h[maxn],a[maxn];
int smax[maxn][22];
void ST(int *a,int n)
{
    for(int i=0;i<n;i++)smax[i][0]=i;
    int k=(int)(log(1.0*n)/log(2.0));
    for(int j=1;j<=k;j++)   
        for(int i=0;i+(1<<j)<=n;i++)  
        {
            if(a[smax[i][j-1]]>a[smax[i+(1<<(j-1))][j-1]])smax[i][j]=smax[i][j-1];
            else smax[i][j]=smax[i+(1<<(j-1))][j-1];
        }
}
int query(int l,int r)  
{  
    int k=(int)(log(1.0*(r-l+1))/log(2.0));  
    if(a[smax[l][k]]<a[smax[r-(1<<k)+1][k]])return a[smax[r-(1<<k)+1][k]];
    return a[smax[l][k]];
}
int main()
{
    while(~scanf("%d",&n))
    {
        int m=0;
        for(int i=1;i<=n;i++)
        {
            scanf("%d%d",&x[i],&y[i]);
            h[m++]=x[i],h[m++]=y[i];
        }
        scanf("%d",&q);
        for(int i=1;i<=q;i++)
        {
            scanf("%d%d",&l[i],&r[i]);
            h[m++]=l[i],h[m++]=r[i];
        }
        sort(h,h+m);
        m=unique(h,h+m)-h;
        //for(int i=0;i<m;i++)printf("%d%c",h[i],i==m-1?'\n':' ');
        memset(a,0,sizeof(a));
        for(int i=1;i<=n;i++)
        {
            int L=lower_bound(h,h+m,x[i])-h;
            int R=lower_bound(h,h+m,y[i])-h;
            a[L]++,a[R+1]--;
        } 
        for(int i=1;i<m;i++)a[i]+=a[i-1];
        ST(a,m);
        for(int i=1;i<=q;i++)
        {
            int L=lower_bound(h,h+m,l[i])-h;
            int R=lower_bound(h,h+m,r[i])-h;
            printf("%d\n",query(L,R));
        }
    }
    return 0;
}