标签：
注意到$b$不变，考虑用归并树来维护这个$b$序列，对于每个节点有序地维护$b$，同时在归并的时候预处理出每个元素在左右儿子里的排名。
在归并树上额外维护区间内$a\geq b$的个数以及赋值标记。
那么在区间赋值的时候，只需要在根节点的$b$数组中做一个二分，然后往下通过预处理的名次数组转移即可，标记下放时也是如此，每次转移复杂度显然是$O(1)$。
时间复杂度$O((n+m)\log n)$。
#include<cstdio>
#include<algorithm>
const int N=100010,M=262150,E=1768950;
int Case,n,m,A,B,C=~(1<<31),i,a[N],b[N],ans,c,d,x,sum;
int st[M],en[M],v[M],tag[M],pl[E],pr[E],pool[E],cur;
void build(int x,int a,int b){
tag[x]=-1;
if(a==b){
st[x]=cur+1;
pool[++cur]=::b[a];
en[x]=cur;
v[x]=::a[a]>=::b[a];
return;
}
int mid=(a+b)>>1;
build(x<<1,a,mid),build(x<<1|1,mid+1,b);
v[x]=v[x<<1]+v[x<<1|1];
int al=st[x<<1],ar=en[x<<1],bl=st[x<<1|1],br=en[x<<1|1];
st[x]=cur+1;
while(al<=ar&&bl<=br)pool[++cur]=pool[al]<pool[bl]?pool[al++]:pool[bl++];
while(al<=ar)pool[++cur]=pool[al++];
while(bl<=br)pool[++cur]=pool[bl++];
en[x]=cur;
al=st[x<<1],bl=st[x<<1|1];
for(int i=st[x];i<=cur;i++){
while(al<=ar&&pool[al]<=pool[i])al++;
while(bl<=br&&pool[bl]<=pool[i])bl++;
pl[i]=al-1,pr[i]=bl-1;
if(pl[i]<st[x<<1])pl[i]=0;
if(pr[i]<st[x<<1|1])pr[i]=0;
}
}
inline void tag1(int x,int p){v[x]=p?p-st[x]+1:0;tag[x]=p;}
inline void pb(int x){
if(tag[x]<0)return;
int p=tag[x];
tag1(x<<1,pl[p]);
tag1(x<<1|1,pr[p]);
tag[x]=-1;
}
void change(int x,int a,int b,int p){
if(c<=a&&b<=d){tag1(x,p);return;}
pb(x);
int mid=(a+b)>>1;
if(c<=mid)change(x<<1,a,mid,pl[p]);
if(d>mid)change(x<<1|1,mid+1,b,pr[p]);
v[x]=v[x<<1]+v[x<<1|1];
}
void ask(int x,int a,int b){
if(c<=a&&b<=d){ans+=v[x];return;}
pb(x);
int mid=(a+b)>>1;
if(c<=mid)ask(x<<1,a,mid);
if(d>mid)ask(x<<1|1,mid+1,b);
v[x]=v[x<<1]+v[x<<1|1];
}
inline int lower(int x){
int l=st[1],r=en[1],mid,t=0;
while(l<=r)if(pool[mid=(l+r)>>1]<=x)l=(t=mid)+1;else r=mid-1;
return t;
}
inline int rnd(){
A=(36969+(ans>>3))*(A&65535)+(A>>16);
B=(18000+(ans>>3))*(B&65535)+(B>>16);
return(C&((A<<16)+B))%1000000000;
}
int main(){
for(scanf("%d",&Case);Case--;printf("%d\n",sum)){
scanf("%d%d%d%d",&n,&m,&A,&B);
for(i=1;i<=n;i++)scanf("%d",&a[i]);
for(i=1;i<=n;i++)scanf("%d",&b[i]);
ans=cur=sum=0;
build(1,1,n);
for(i=1;i<=m;i++){
c=rnd()%n+1,d=rnd()%n+1,x=rnd()+1;
if(c>d)std::swap(c,d);
if((c+d+x)&1)change(1,1,n,lower(x));
else{
ans=0;
ask(1,1,n);
sum=(1LL*ans*i+sum)%1000000007;
}
}
}
return 0;
}
标签：
原文地址：http://www.cnblogs.com/clrs97/p/5693741.html
