2019牛客暑期多校训练营（第七场）E.Find the median （动态中位数）（权值线段树+点代表区间）

题意：Q次操作，每次给你一个区间[l,r],区间内[l,r]的值都出现+1次，问你中位数是多少

题解：点代表区间，找第K大，这道题我搞了好久，问了好多人才明白究竟是啥意思。

          首先可以每次把【l,r】的 r+1, 按照【l,r+1】建树，好处多多，这样就可以让你的结点化为一个左臂右开的区间，并且这个区间可以都依附于左端点的这个节点，可以直接满足加减性，

          区间更新+查找第K大

          比如样例建图如下：

    该图就完美解释代码了

细节：建树时，只用考虑离散化数组-1的长度，因为离散化数组最后一位一定是r+1的值【r+1，+00】；

          还有就是可以在建树时，直接处理区间长度

          对给的区间【l,r】离散化时要注意，因为可能给的l,r两值树的值域没有，得离散化找对应的树中位置。

#include <bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn = 1e6 + 5;
const int INF = 0x3f3f3f3f;
const ll MOD = 1e9 + 7;
ll L[maxn];ll R[maxn];
ll X[maxn]; ll Y[maxn];
vector<ll>v;
ll sum[maxn*4];ll num[maxn*4];ll len[maxn*4];
int add[maxn*4];
void push_down(int p){
	if(add[p]){
		sum[p<<1]+=add[p]*len[p<<1];
		sum[p<<1|1]+=add[p]*len[p<<1|1];
		add[p<<1]+=add[p];
		add[p<<1|1]+=add[p];
		add[p]=0;
	}
}
void build(int l,int r,int p){
	if(l==r){
		sum[p]=0;add[p]=0; 
		len[p]=v[l]-v[l-1];//每个结点区间长度 
	//	cout<<l<<" : "<<r<<" "<<len[p]<<endl;
		return ;
	}
	int mid=(l+r)/2;
	build(l,mid,p<<1);
	build(mid+1,r,p<<1|1);
	len[p]=len[p<<1]+len[p<<1|1];
}
void update(int l,int r,int L,int R,int p){
	if(L<=l&&r<=R) {
		sum[p]+=len[p];//每次区间更新 
		add[p]++;
	//	cout<<"插入时： "<<l<<" : "<<r<<" "<<sum[p]<<endl;
		return ;
	}
	push_down(p);
	int mid=(l+r)/2;
	if(L<=mid) update(l,mid,L,R,p<<1);
	if(R>mid) update(mid+1,r,L,R,p<<1|1);
	sum[p]=sum[p<<1]+sum[p<<1|1];
}
ll query(int l,int r,int p,ll k){
	if(l==r){
		ll cnt=sum[p]/len[p];//每组几个 
		ll tmp=v[l-1]+(k-1)/cnt;
		return tmp;
	}
	push_down(p);
	int mid=(l+r)/2;
	ll Lsum=sum[p<<1];
//	cout<<"查询： "<<p<<" "<<l<<" "<<r<<" "<<Lsum<<endl;
	if(k<=Lsum) return query(l,mid,p<<1,k);
	else return query(mid+1,r,p<<1|1,k-Lsum);
}
int main(){
	int n;   cin>>n;
	ll x1,x2,a1,b1,c1,m1;           ll y1,y2,a2,b2,c2,m2;
	cin>>x1>>x2>>a1>>b1>>c1>>m1;    cin>>y1>>y2>>a2>>b2>>c2>>m2;
	L[1]=min(x1,y1)+1; L[2]=min(x2,y2)+1; X[1]=x1; X[2]=x2;
	R[1]=max(x1,y1)+1; R[2]=max(x2,y2)+1; Y[1]=y1; Y[2]=y2;
	v.push_back(L[1]); v.push_back(R[1]+1);
	v.push_back(L[2]); v.push_back(R[2]+1);
	for(int i=3;i<=n;i++){
		X[i]=(a1*X[i-1]+b1*X[i-2]+c1)%m1;
		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;
		v.push_back(L[i]);
		v.push_back(R[i]+1);
	}
	sort(v.begin(),v.end());
	v.erase(unique(v.begin(),v.end()),v.end());
	for(int i=1;i<=n;i++){
		L[i]=lower_bound(v.begin(),v.end(),L[i])-v.begin()+1;
		R[i]=lower_bound(v.begin(),v.end(),R[i]+1)-v.begin()+1;//满足左臂右开
	}
	int size=v.size(); 	
    build(1,size-1,1);
    ll ps=0;
	for(int i=1;i<=n;i++){
		update(1,size-1,L[i],R[i]-1,1);//更新时按照原先的[l,r]区间 
		ps+=v[R[i]-1]-v[L[i]-1];//直接满足该区间长度 
		ll ans=query(1,size-1,1,(ps+1)/2);
		cout<<ans<<endl;
	}
	return 0;
}