【例题】【主席树】

1、
NKOJ2664 静态区间第K小数
Time Limit:40000MS Memory Limit:234567K
Total Submit:56 Accepted:39
Case Time Limit:3000MS

Description
给定一个由N个数组成的序列{A1,A2,…,AN}
每次提问(l,r,k)表示提问序列中{Al,..,Ar}中第k小的数的值。

Input
第一行两个数N,Q，表示有N个数，Q次提问。
第二行N个数，为给定的序列。
接下来Q行每行三个数(l,r,k)，表示提问。

Output
对每个提问输出一行，只有一个数，为对应的答案。

Sample Input
7 3
1 5 2 6 3 7 4
2 5 3
4 4 1
1 7 3

Sample Output
5
6
3

Hint
对于40%的数据
1<=N<=1000
1<=Q<=1000
对于60%的数据（包含上述40%的数据），Ai是1到N的一个排列
对于100%的数据
1<=N<=100000
1<=Q<=200000
1<=Ai<=2*109
1<=k<=r-l+1

Source 感谢nodgd提供数据

思路：
先离散化。
前缀和思想，对于n个区间 [1,i] 分别建立n个线段树，维护区间内各数的个数
使用可持久化二叉树的思想减少空间
就变成可持续化线段树了….

注：理论上需先建一颗完整的树，但该题中为空树，不建也没有影响

#include<cstdio>
#include<algorithm>
#include<iostream>
using namespace std;
const int need=100003;
const int deep_=17;

//...............................................
inline int in_()
{
    char t=getchar();int d;
    while(t<'0'|t>'9') t=getchar();
    for(d=0;!(t<'0'|t>'9');t=getchar()) d=(d<<1)+(d<<3)+t-'0';
    return d;
}
inline void out_(int x)
{
    if(x>=10) out_(x/10);
    putchar(x%10+'0');
}
//...............................................
typedef int sdx[(need*(deep_+1)+need)*2];

int root[need],tot;
sdx ls,rs,cnt;//,aa,bb;
int a[need],b[need];

int d;

int add(int pre,int l,int r)
{
    int s=++tot;
    ls[s]=ls[pre],rs[s]=rs[pre],cnt[s]=cnt[pre]+1;//,aa[s]=aa[pre],bb[s]=bb[pre];
    if(l==r) return s;
    int mid=(l+r)>>1;
    if(d<=mid) ls[s]=add(ls[pre],l,(l+r)>>1);
    else rs[s]=add(rs[pre],(l+r)/2+1,r);
    return s;
}

int ask(int x,int y,int l,int r,int k)
{
    if(l==r) return l;
    int cnt_=cnt[ls[y]]-cnt[ls[x]];
    if(cnt_>=k) return ask(ls[x],ls[y],l,(l+r)>>1,k);
    else return ask(rs[x],rs[y],(l+r)/2+1,r,k-cnt_);
}

int main()
{
    int n,m;n=in_(),m=in_();
    for(int i=1;i<=n;i++) b[i]=a[i]=in_();
    sort(b+1,b+1+n);
    int tot_=unique(b+1,b+1+n)-&b[1];
    root[0]=tot=1;//build(1,tot_);
    for(int i=1;i<=n;i++)
    {
        d=lower_bound(b+1,b+1+tot_,a[i])-&b[0];
        root[i]=add(root[i-1],1,tot_);
    }
    for(int i=1,x,y,k;i<=m;i++)
    {
        x=in_(),y=in_(),k=in_();
        d=ask(root[x-1],root[y],1,tot_,k);
        out_(b[d]),putchar(10);
    }
}

2、
NKOJ3677 Couriers
Time Limit:30000MS Memory Limit:265536K
Total Submit:13 Accepted:10

Description
给一个长度为n的序列a。1≤a[i]≤n。
m组询问，每次询问一个区间[l,r]，是否存在一个数在[l,r]中出现的次数大于(r-l+1)/2。如果存在，输出这个数，否则输出0。

Input
第一行两个数n，m。
第二行n个数，a[i]。
接下来m行，每行两个数l,r，表示询问[l,r]这个区间

Output
m行，每行对应一个答案。

Sample Input
7 5
1 1 3 2 3 4 3
1 3
1 4
3 7
1 7
6 6

Sample Output
1
0
3
0
4

Hint
n,m≤500000

Source
poi 2014

思路：
基本同上：线段树依然维护区间内所有数字出现次数

查询时，
若两线段树该节点总数之差<=k一定无解
若该节点左儿子节点总数之差>k 左子树可能有解，右子树一定无解
…

#include<cstdio>
#include<algorithm>
#include<iostream>
using namespace std;
const int need=100003;
const int deep_=17;

//...............................................
inline int in_()
{
    char t=getchar();int d;bool mark=0;
    while(t<'0'|t>'9') {if(t=='-') mark=1;t=getchar();}
    for(d=0;!(t<'0'|t>'9');t=getchar()) d=(d<<1)+(d<<3)+t-'0';
    return mark ? -d : d;
}
inline void out_(int x)
{
    if(x<0) {putchar('-');x=-x;}
    if(x>=10) out_(x/10);
    putchar(x%10+'0');
}
//...............................................
typedef int sdx[(need*(deep_+1)+need)*2];

int root[need],tot;
sdx ls,rs,cnt;
int a[need];
int d,k;

int add(int pre,int l,int r)
{
    int s=++tot;
    ls[s]=ls[pre],rs[s]=rs[pre],cnt[s]=cnt[pre]+1;
    if(l!=r) 
    {
        int mid=(l+r)>>1;
        if(d<=mid) ls[s]=add(ls[pre],l,(l+r)>>1);
        else rs[s]=add(rs[pre],(l+r)/2+1,r);
    }
    return s;
}

int ask(int x,int y,int l,int r)
{
    if(cnt[y]-cnt[x]<=k) return 0;
    if(l==r) return l;
    else
    {
        if(cnt[ls[y]]-cnt[ls[x]]>k) return ask(ls[x],ls[y],l,(l+r)>>1);
        if(cnt[rs[y]]-cnt[rs[x]]>k) return ask(rs[x],rs[y],(l+r)/2+1,r);
        else return 0;
    }
}

int main()
{
    int n,m;n=in_(),m=in_();
    root[0]=tot=1;//build(1,n);
    for(int i=1;i<=n;i++)
    {
        d=a[i]=in_();
        root[i]=add(root[i-1],1,n);
    }
    int temp;
    for(int i=1,x,y;i<=m;i++)
    {
        x=in_(),y=in_();
        k=(y-x+1)>>1;
        temp=ask(root[x-1],root[y],1,n);
        if(temp!=0) out_(temp),putchar(10);
        else putchar('0'),putchar(10);
    }
}