整体二分初识

菊苣是怎么一天学会CDQ和二分的啊我哭了（还不是自己太懒了）

和CDQ分治相似，也有按照时间将修改与询问排序达到分治的操作。

以后补详细理解（因为目前还没啥理解。。。。

K-th Number

POJ链接：http://poj.org/problem?id=2104

居然因为结构体构造函数死循环了。。。构造函数还是不要重名吧，尽量用ll,rr啥的

#include<map>
#include<set>
#include<queue>
#include<cmath>
#include<stack>
#include<ctime>
#include<cstdio>
#include<vector>
#include<string>
#include<sstream>
#include<cstdlib>
#include<cstring>
#include<cassert>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int MAXN=1e5+10;
#define pi  acos(-1.0)
#define INF  0x3f3f3f3f
#define mod   1000000009
#define endll printf("\n")
#define s1(a) scanf("%d",&a)
#define lowbit(x)  ((x)&(-x))
#define s2(a,b) scanf("%d %d",&a,&b)
#define mem(a,b) memset(a,b,sizeof(a))
#define s3(a,b,c) scanf("%d %d %d",&a,&b,&c)
int m,maxn,tot;
struct IN
{
   int l,r,v,pos,op;
   IN(){}
   IN(int ll,int rr,int vv,int pp,int oo)
   {l=ll,r=rr,v=vv,pos=pp,op=oo;}//又是一个教训。。。。
}q[MAXN<<1],q1[MAXN<<1],q2[MAXN<<1];
int sum[MAXN],ans[MAXN];
void fix(int x,int v){for(;x<=maxn;x+=lowbit(x)) sum[x]+=v;}
int fund(int x){int re=0;for(;x;x-=lowbit(x)) re+=sum[x];return re;}
void solve(int l,int r,int L,int R)//分别表示答案可能存在的区间，答案在此区间的操作区间
{
    if(l>r||L>R) return;
    if(l==r)//此时若此区间还有查询操作，那么此区间即为这些操作的答案
    {
        for(int i=L;i<=R;i++)
            if(q[i].op) ans[q[i].pos]=l;
        return;
    }
    int m=(l+r)>>1,cnt1=0,cnt2=0;
    for(int i=L;i<=R;i++)//把原区间内的操作分到左右两区间内
    {
        if(q[i].op)//如果是查询操作
        {
            int val=fund(q[i].r)-fund(q[i].l-1);
            if(val>=q[i].v) q1[++cnt1]=q[i];//说明左区间的数字数量已经大于k了，那么第k大数字一定在左区间
            else q[i].v-=val,q2[++cnt2]=q[i];//否则就在右区间，这时要减去左区间对他的贡献
        }
        else//如果是修改操作
        {
            if(q[i].v<=m) q1[++cnt1]=q[i],fix(q[i].pos,1);//只对位于左区间的修改操作进行修改
            else q2[++cnt2]=q[i];//将修改操作丢到右区间
        }
    }
    for(int i=1;i<=cnt1;i++)
        if(!q1[i].op) fix(q1[i].pos,-1);//
    for(int i=1;i<=cnt1;i++) q[L+i-1]=q1[i];
    for(int i=1;i<=cnt2;i++) q[L+cnt1+i-1]=q2[i];
    solve(l,m,L,L+cnt1-1);solve(m+1,r,L+cnt1,R);
}
int main()
{
    s2(maxn,m);tot=0;
    for(int i=1,x;i<=maxn;i++)
        s1(x),q[++tot]=IN(0,0,x,i,0);
    for(int i=1,l,r,k;i<=m;i++)
        s3(l,r,k),q[++tot]=IN(l,r,k,i,1);
    solve(-INF,INF,1,tot);
    for(int i=1;i<=m;i++) printf("%d\n",ans[i]);
    return 0;
}