2019牛客暑期多校训练营（第七场）-E Find the median （线段树+离散化 区间为点）

题目链接：https://ac.nowcoder.com/acm/contest/887/E

题意：给出L[i],R[i]，每次添加L[i]...R[i]，求出此时的中位数。

思路：因为添加的数范围为1e9，首先想到要用离散化，这里是用一个点来表示一个区间。

将右区间加一的主要目的是优化处理，将区间最后一个元素并入到前面，自己模拟一下就能明白啦。

然后用线段树维护每个节点所包含的元素个数，用懒惰标记laz表示加的次数，然后每次查询时若左子树的元素个数足够则在左子树查询，否则在右子树查询。

详见代码：

#include<cstdio>
#include<algorithm>
#include<vector>
using namespace std;

typedef long long LL;
const int maxn=4e5+;

struct node{
    int l,r,laz;
    LL sz;
}tr[maxn<<];

int n;
LL X[maxn],Y[maxn],A1,A2,B1,B2,C1,C2,M1,M2;
vector<int> vc;

void build(int v,int l,int r){
    tr[v].l=l,tr[v].r=r;
    if(l==r){
        return;
    }
    int mid=(l+r)>>;
    build(v<<,l,mid);
    build(v<<|,mid+,r);
}

void pushdown(int v){
    tr[v<<].sz+=(vc[tr[v<<].r+]-vc[tr[v<<].l])*tr[v].laz;
    tr[v<<].laz+=tr[v].laz;
    tr[v<<|].sz+=(vc[tr[v<<|].r+]-vc[tr[v<<|].l])*tr[v].laz;
    tr[v<<|].laz+=tr[v].laz;
    tr[v].laz=;
}

void update(int v,int l,int r){
    if(l<=tr[v].l&&r>=tr[v].r){
        tr[v].sz+=(vc[tr[v].r+]-vc[tr[v].l]);
        tr[v].laz+=;
        return;
    }
    if(tr[v].laz) pushdown(v);
    int mid=(tr[v].l+tr[v].r)>>;
    if(l<=mid) update(v<<,l,r);
    if(r>mid) update(v<<|,l,r);
    tr[v].sz=tr[v<<].sz+tr[v<<|].sz;
}

int query(int v,LL k){
    if(tr[v].l==tr[v].r){
        int tmp=tr[v].sz/(vc[tr[v].l+]-vc[tr[v].l]);
        return vc[tr[v].l]+(k-)/tmp;
    }
    if(tr[v].laz) pushdown(v);
    if(k<=tr[v<<].sz) return query(v<<,k);
    else return query(v<<|,k-tr[v<<].sz);
}

int main(){
    scanf("%d",&n);
    scanf("%lld%lld%lld%lld%lld%lld",&X[],&X[],&A1,&B1,&C1,&M1);
    scanf("%lld%lld%lld%lld%lld%lld",&Y[],&Y[],&A2,&B2,&C2,&M2);
    for(int i=;i<=n;++i){
        X[i]=(A1*X[i-]%M1+B1*X[i-]%M1+C1)%M1;
        Y[i]=(A2*Y[i-]%M2+B2*Y[i-]%M2+C2)%M2;
    }
    for(int i=;i<=n;++i){
        ++X[i],++Y[i];
        if(X[i]>Y[i]) swap(X[i],Y[i]);
        vc.push_back(X[i]),vc.push_back(Y[i]+);
    }
    sort(vc.begin(),vc.end());
    vc.erase(unique(vc.begin(),vc.end()),vc.end());
    LL sum=;
    int cnt=vc.size();
    build(,,cnt-);
    for(int i=;i<=n;++i){
        sum+=Y[i]-X[i]+;
        int x,y;
        x=lower_bound(vc.begin(),vc.end(),X[i])-vc.begin();
        y=lower_bound(vc.begin(),vc.end(),Y[i]+)-vc.begin();
        update(,x,y-);
        printf("%d\n",query(,(sum+)>>));
    }
    return ;
}