2019牛客暑期多校（第七场）E题Find the median

题目链接、

题意：

给你一个空数组，n次操作，每次往里面填入里的每个数，也就是l，l+1，，，，r，然后问你中位数是谁

题解：

考虑树状数组维护区间、

众所周知，如果差分的话，这种区间加一的操作只需要加俩点就行了，但是要求中位数的话，需要树状数组+二分才可以

一般的中位数是树状数组维护原数组，然后二分前缀和，这里并不能这么用，因为不能维护原数组。

所以这里考虑区间影响

询问pos，这个位置的时候，想办法把小于等于pos这个位置的所有操作统计一下，这就算出来了有多少个数比pos小，

又因为知道总个数，所以可以二分这个位置，得到答案

问题是怎么统计pos之前的操作，考虑

如果之前操作的区间如果右端点小于等于pos的话，相当于整个区间被包含，

所以我们可以在离散化之后的l位置减一个L,在r位置加一个R+1，这样每次我们询问如果包含这个区间的话，那么正好是

R-L+1，可以正好统计出来

另一种情况是，区间被pos分割，这样需要统计的是区间左端点到pos这段的和

因为之前我们将l位置减掉了一个L，这一步相当于什么呢，相当于在区间全部减掉了个1，

而我们要求的是

显然我们只需要加上pos+1即可（），有多少个这样的分割区间，就加多少个pos+1

统计有多少个这样的区间，容易得出，在l的位置+1，r的位置-1即可

这样树状数组维护两个数组，二分答案就可以了

 

/*author:revolIA*/
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 2e6+7;
ll Sum[maxn],Tree[maxn],n;
void add(int pos,ll x,ll y){
    while(pos<maxn)Sum[pos] += x,Tree[pos] += y,pos += pos&-pos;
}
ll query(ll A[],int pos,ll ans = 0){
    while(pos)ans += A[pos],pos-=pos&-pos;
    return ans;
}
ll x[maxn],y[maxn],a1,b1,c1,m1,a2,b2,c2,m2;
ll l[maxn],r[maxn];
vector<ll> vec;
int main(){
    scanf("%lld",&n);
    scanf("%lld%lld%lld%lld%lld%lld",&x[1],&x[2],&a1,&b1,&c1,&m1);
    scanf("%lld%lld%lld%lld%lld%lld",&y[1],&y[2],&a2,&b2,&c2,&m2);
    for(int i=1;i<=n;i++){
        if(i>2)x[i]=(a1*x[i-1]+b1*x[i-2]+c1)%m1;
        if(i>2)y[i]=(a2*y[i-1]+b2*y[i-2]+c2)%m2;
        l[i] = min(x[i],y[i])+1;
        r[i] = max(x[i],y[i])+1;
        vec.push_back(l[i]);
        vec.push_back(r[i]);
    }
    sort(vec.begin(),vec.end());
    vec.erase(unique(vec.begin(),vec.end()),vec.end());
    ll cnt = 0;
    for(int i=1;i<=n;i++){
        cnt += r[i]-l[i]+1;
        int L = lower_bound(vec.begin(),vec.end(),l[i])-vec.begin()+1;
        int R = lower_bound(vec.begin(),vec.end(),r[i])-vec.begin()+1;
        add(L,1ll,-l[i]),add(R,-1ll,r[i]+1);
        L = 1,R = 1e9;
        int sul = -1;
        while(L<=R){
            int mid = L+R>>1;
            int pos = upper_bound(vec.begin(),vec.end(),mid)-vec.begin();
            ll Cnt = query(Tree,pos)+query(Sum,pos)*(mid+1);
            if(Cnt*2 >= cnt){
                R = mid-1,sul = mid;
            }else
                L = mid+1;
        }
        printf("%d\n",sul);
    }
    return 0;
}