19 hdu多校 Welcome Party //主席树

http://acm.hdu.edu.cn/contests/contest_showproblem.php?cid=857&pid=1005

题意：n个人，有xi和yi属性，将n个人分成两个队伍，分到a队只有xi属性，分到b队只有yi属性，每个队伍得至少一人，求分队后最小的abs(a队xi的最大值-b队yi的最大值)。

思路：主席树，枚举x最大是谁，考虑y最大是谁。枚举到x最大的人的下标为p时，当确定p之后，属性xi大于p.xi的人只能在b队。属性xi小于p.xi的人，考虑yi和p.xi最接近的。（大于等于或者小于）

/*   (o O o)
 *   Author : Rshs
 *   Data : 2019-08-21-14.10
 */
#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define FI first
#define SE second
#define MP make_pair
#define PII pair<int,int>
const LL mod = 1e9+7;
const int MX = 1e5+5;
const int max_n =1e5+5,logn=20;
int root[max_n*logn],ls[max_n*logn],rs[max_n*logn]; //下标为i的树他的左右子树的对应sum的下标
int sum[max_n*logn];//下标为i的区间(上面的区间对应进来)的数字数量
int cnt;//树的数量
struct no{
    LL x,y;
}a[MX];
bool operator<(const no &x,const no &y){
    return x.x<y.x;
}
vector<LL>Y;
void build(int l,int r,int &now){
    now=++cnt;//对应新开的区域
    sum[now]=0;
    if(l==r) return ;
    int mid=(l+r)>>1;
    build(l,mid,ls[now]);
    build(mid+1,r,rs[now]);
}
void update(int l,int r,int &now,int last,int pos){
    now=++cnt;//对应新开的区域
    ls[now]=ls[last];
    rs[now]=rs[last];
    sum[now]=sum[last]+1;
    if(l==r) return;
    int mid=(l+r)>>1;
    if(pos<=mid) update(l,mid,ls[now],ls[last],pos);
    else update(mid+1,r,rs[now],rs[last],pos);
}
int query(int l,int r,int t1,int t2,int num){//区间第k大
    if(l==r) return l;
    int mid=(l+r)>>1;
    int cc=sum[rs[t2]]-sum[rs[t1]]; //第一棵树左儿子减第二颗树左儿子
    if(num<=cc) query(mid+1,r,rs[t1],rs[t2],num);
    else query(l,mid,ls[t1],ls[t2],num-cc); //减去已有的数量
}
LL Query(int L,int R,int l,int r,int t1,int t2){ //[L,R]里面元素个数
    if(l>=L&&r<=R)
        return 1LL*sum[t2]-sum[t1];
    int mid=(l+r)>>1;
    LL re=0;
    if(mid>=L) re+=Query(L,R,l,mid,ls[t1],ls[t2]);
    if(mid<R) re+=Query(L,R,mid+1,r,rs[t1],rs[t2]);
    return re;
}
int main(){
    //out<<LLONG_MAX;
    int T;cin>>T;while(T--){
        cnt=0;Y.clear();
        int n;cin>>n;for(int i=1;i<=n;i++) scanf("%I64d%I64d",&a[i].x,&a[i].y);
        sort(a+1,a+1+n);
        for(int i=1;i<=n;i++) Y.push_back(a[i].y);
        sort(Y.begin(),Y.end());Y.erase(unique(Y.begin(),Y.end()),Y.end());
        int sz=Y.size();
       // build(1,sz,root[0]);
        for(int i=1;i<=n;i++){
            int id=lower_bound(Y.begin(),Y.end(),a[i].y)-Y.begin()+1;
            update(1,sz,root[i],root[i-1],id);
        }
        LL ans=LLONG_MAX;
        for(int i=1;i<=n;i++){
            LL MX=LLONG_MIN/2;
            if(i!=n){
                LL g=Y[query(1,sz,root[i],root[n],1)-1];
                MX=g;
            }
            LL tMX=MX;
            if(i!=1){
                int id=upper_bound(Y.begin(),Y.end(),a[i].x)-Y.begin()-1+1;
                int cc;
                if(id<1)cc=0;
                else cc=Query(1,id,1,sz,root[0],root[i-1]);//小于等于a[i].x
                int dd=Query(1,sz,1,sz,root[0],root[i-1]);
                cc=dd-cc;
                for(int j=cc;j<=cc+1;j++){
                    if(j>=1&&j<=dd){
                        LL g=Y[query(1,sz,root[0],root[i-1],j)-1];
                        if(g>tMX&&abs(g-a[i].x)<abs(a[i].x-MX)) MX=g;
                    }
                }
            }
            ans=min(ans,abs(a[i].x-MX));
        }
        cout<<ans<<'\n';
    }
    return 0;
}
/*
10
4
1 9
7 6
5 5
4 7*/