2019牛客暑期多校训练营（第七场） E 线段树+离散化

题目传送门

题意：按照一定的公式给出若干个$<l,r>$，每次往一个序列中加上l到r的数字，并输出中位数。

思路：需要将每个$区间$离散化，比如把$[1,2]$变成$[1,3)$，也就是$[1,2)$和$[2,3)$，这样做才能完整的表达区间，否则如[2,2]这样的区间就会出现问题。

　　所以我们将每一个$[l,r+1)$离散化，设$（x）$代表x离散化后的数字，每次更新的时候，右区间应该是$（r+1）-1$。比如原来只有一个区间$[2,3]$，我们表示成$[2,4)$,离散化后我们更新的区间是$[1,1]$,加的数字的个数是$ve[2]-ve[1]$，ve是原数组。

　　查询和更新思路一样。

#include<bits/stdc++.h>
#define clr(a,b) memset(a,b,sizeof(a))
#define pb(a)  push_back(a)
#define rep(i,x,n) for(int i=x;i<=n;++i)
#define dep(i,n,x) for(int i=n;i>=x;--i)
using namespace std;
typedef long long ll;
int n,m;
const int inf=0x3f3f3f3f;
const int maxn=4e5+;
ll x[maxn],y[maxn],A1,B1,C1,M1,A2,B2,C2,M2;
vector<ll >ve;
ll sz[maxn<<],lazy[maxn<<];
void add(int o,int l,int r,ll f){
    sz[o]+=(ve[r+]-ve[l])*f;
    lazy[o]+=f;
}
void pushdown(int o,int l,int r){
    if(!lazy[o]||l==r)return;
    int mid=(l+r)>>;
    add(o<<,l,mid,lazy[o]);
    add(o<<|,mid+,r,lazy[o]);
    lazy[o]=;
}
void update(int o,int l,int r,int ql,int qr){
    if(ql<=l&&r<=qr){
        add(o,l,r,);
        return;
    }
    int mid=(l+r)>>;
    pushdown(o,l,r);
    if(ql<=mid)update(o<<,l,mid,ql,qr);
    if(mid<qr)update(o<<|,mid+,r,ql,qr);
    sz[o]=sz[o<<]+sz[o<<|];
}

ll query(int o,int l,int r,int num){
    if(l==r){
        ll tot=sz[o]/(ve[l+]-ve[l]);
        return ve[l]+(num-)/tot;
    }
    pushdown(o,l,r);
    int mid=(l+r)>>;
    if(sz[o<<]>=num)return query(o<<,l,mid,num);
    return query(o<<|,mid+,r,num-sz[o<<]);
}

int main(){
    int n;
    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);
       rep(i,,n){
        x[i]=(A1*x[i-]+B1*x[i-]+C1)%M1;
        y[i]=(A2*y[i-]+B2*y[i-]+C2)%M2;
    }

    rep(i,,n){
        x[i]++,y[i]++;
        if(x[i]>y[i]) swap(x[i],y[i]);
        ve.push_back(x[i]); ve.push_back(y[i]+);
    }
    sort(ve.begin(),ve.end());
    ve.erase(unique(ve.begin(),ve.end()),ve.end());
    int cnt=ve.size();
    rep(i,,n){

          x[i]=lower_bound(ve.begin(),ve.end(),x[i])-ve.begin();
        y[i]=lower_bound(ve.begin(),ve.end(),y[i]+)-ve.begin();
        update(,,cnt,x[i],y[i]-);
        printf("%lld\n",query(,,cnt,(sz[]+)/));
    }
}