BZOJ2872 : 优莱卡

时间:2023-12-21 15:38:32

$A$类数据:直接输出$B[l]$,$O(n+m)$。
$B$类数据:ST表求区间最小值,$O(n\log n+m)$。
$C$类数据:可持久线段树求区间$k$小值,$O((n+m)\log n)$。
剩下的数据先用可持久线段树求出$a,b$的实际值,转化为矩形内$B$值的询问。
$D$类数据:即求矩形$B$最小值,对序列分治,对于序列$[l,r]$,用扫描线+线段树处理所有经过$mid$的询问,$O(n\log^2n+m\log n)$。
$E,F$类数据:即求矩形$B$第$k$小值,可持久化线段树套权值线段树即可,$O((n+m)\log^2n)$。

#include<cstdio>
const int N=30010,M=N*18,BUF=30000000,OUT=10000000;
char Buf[BUF],*buf=Buf,Out[OUT],*ou=Out;int Outn[30],Outcnt;
inline void write(int x){
if(!x)*ou++=48;
else{
for(Outcnt=0;x;x/=10)Outn[++Outcnt]=x%10+48;
while(Outcnt)*ou++=Outn[Outcnt--];
}
}
inline void writeln(int x){write(x);*ou++='\n';}
inline void read(int&a){for(a=0;*buf<48;buf++);while(*buf>47)a=a*10+*buf++-48;}
inline int readstr(char*a){
int n=0;
while(*buf<33)buf++;while(*buf>32)*a++=*buf++,n++;
return n;
}
int n,m,e,i,j,l,r,c,d,k,a[N],b[N],last,v[M],sl[M],sr[M],tot,T[N];char type[9];
namespace SubA{
void Main(){
while(m--)read(l),read(r),read(c),read(d),read(k),writeln(b[l]);
}
}
int ins(int x,int a,int b,int c){
int y=++tot;
v[y]=v[x]+1;
if(a==b)return y;
int mid=(a+b)>>1;
if(c<=mid)sl[y]=ins(sl[x],a,mid,c),sr[y]=sr[x];
else sl[y]=sl[x],sr[y]=ins(sr[x],mid+1,b,c);
return y;
}
int kth(int x,int y,int k){
int a=1,b=n,mid,t;
while(a<b){
mid=(a+b)>>1;
t=v[sl[x]]-v[sl[y]];
if(k<=t)b=mid,x=sl[x],y=sl[y];else k-=t,a=mid+1,x=sr[x],y=sr[y];
}
return a;
}
namespace SubB{
int Log[N],f[20][N];
inline int min(int a,int b){return a<b?a:b;}
inline int ask(int l,int r){
int k=Log[r-l+1];
return min(f[k][l],f[k][r-(1<<k)+1]);
}
void Main(){
for(i=2;i<=n;i++)Log[i]=Log[i>>1]+1;
for(i=1;i<=n;i++)f[0][i]=a[i];
for(j=1;j<=Log[n];j++)for(i=1;i+(1<<j)-1<=n;i++)f[j][i]=min(f[j-1][i],f[j-1][i+(1<<(j-1))]);
while(m--){
read(l),read(r),read(c),read(d),read(k);
writeln(ask(l,r));
}
}
}
namespace SubC{
void Main(){
while(m--){
read(l),read(r),read(c),read(d),read(k);
writeln(kth(T[r],T[l-1],c));
}
}
}
namespace SubD{
const int MAXQ=500010;
struct E{int l,r,a,b,ans;}q[MAXQ];
int pool[N],cnt[N],v[MAXQ<<1],val[70000],M;
inline void up(int&a,int b){a>b?(a=b):0;}
inline int ask(int x,int y){
int t=N;
for(x+=M-1,y+=M+1;x^y^1;x>>=1,y>>=1){
if(~x&1)up(t,val[x^1]);
if(y&1)up(t,val[y^1]);
}
return t;
}
inline void change(int x,int y){for(x+=M;x;x>>=1)up(val[x],y);}
inline void clear(int x){for(x+=M;x;x>>=1)val[x]=N;}
void build(int n){
for(M=1;M<n+2;M<<=1);
int i;
for(i=1;i<=n;i++)v[i+M]=a[i];
for(i=M-1;i;i--)v[i]=v[i<<1]+v[i<<1|1];
}
void solve(int l,int r){
if(l>r)return;
int i,j,k,mid=(l+r)>>1;
solve(l,mid-1),solve(mid+1,r);
for(i=mid;i>=l;i--){
change(a[i],b[i]);
for(j=cnt[i];j>cnt[i-1];j--){
k=v[j];
if(q[k].l<=mid&&q[k].r>=mid&&q[k].l>=l&&q[k].r<=r)up(q[k].ans,ask(q[k].a,q[k].b));
}
}
for(i=mid;i>=l;i--)clear(a[i]);
for(i=mid;i<=r;i++){
change(a[i],b[i]);
for(j=cnt[i];j>cnt[i-1];j--){
k=v[j];
if(q[k].l<=mid&&q[k].r>=mid&&q[k].l>=l&&q[k].r<=r)up(q[k].ans,ask(q[k].a,q[k].b));
}
}
for(i=mid;i<=r;i++)clear(a[i]);
}
void Main(){
for(i=1;i<=m;i++){
read(l),read(r),read(c),read(d),read(k);
c=kth(T[r],T[l-1],c);
d=kth(T[r],T[l-1],d);
q[i].l=l,q[i].r=r,q[i].a=c,q[i].b=d,q[i].ans=N;
cnt[l]++,cnt[r]++;
}
for(i=1;i<=n;i++)cnt[i]+=cnt[i-1];
for(i=1;i<=n;i++)pool[i]=cnt[i];
for(i=1;i<=m;i++)v[pool[q[i].l]--]=i,v[pool[q[i].r]--]=i;
for(M=1;M<n+2;M<<=1);
for(i=1;i<=n+M;i++)val[i]=N;
solve(1,n);
for(i=1;i<=m;i++)writeln(q[i].ans);
}
}
namespace SubF{
const int MAXN=9720000;
int sonl[M],sonr[M],id[M],root[N],tot;
int L[MAXN],R[MAXN],V[MAXN],TOT;
int q[N][2],cnt;
int INS(int x,int a,int b,int c){
int y=++TOT;
V[y]=V[x]+1;
if(a==b)return y;
int mid=(a+b)>>1;
if(c<=mid)L[y]=INS(L[x],a,mid,c),R[y]=R[x];else L[y]=L[x],R[y]=INS(R[x],mid+1,b,c);
return y;
}
int ins(int x,int a,int b,int c,int d){
int y=++tot;
id[y]=INS(id[x],0,n,d);
if(a==b)return y;
int mid=(a+b)>>1;
if(c<=mid)sonl[y]=ins(sonl[x],a,mid,c,d),sonr[y]=sonr[x];
else sonl[y]=sonl[x],sonr[y]=ins(sonr[x],mid+1,b,c,d);
return y;
}
void ask(int x,int y,int a,int b,int c,int d){
if(c<=a&&b<=d){
q[++cnt][0]=id[x];
q[cnt][1]=id[y];
return;
}
int mid=(a+b)>>1;
if(c<=mid)ask(sonl[x],sonl[y],a,mid,c,d);
if(d>mid)ask(sonr[x],sonr[y],mid+1,b,c,d);
}
inline int query(int x,int y,int l,int r,int k){
cnt=0;
ask(x,y,1,n,l,r);
int a=0,b=n,mid,t,i;
while(a<b){
mid=(a+b)>>1;
t=0;
for(i=1;i<=cnt;i++)t+=V[L[q[i][0]]]-V[L[q[i][1]]];
if(k<=t){
b=mid;
for(i=1;i<=cnt;i++)q[i][0]=L[q[i][0]],q[i][1]=L[q[i][1]];
}else{
k-=t,a=mid+1;
for(i=1;i<=cnt;i++)q[i][0]=R[q[i][0]],q[i][1]=R[q[i][1]];
}
}
return a;
}
void Main(){
for(i=1;i<=n;i++)root[i]=ins(root[i-1],1,n,a[i],b[i]);
while(m--){
read(l),read(r),read(c),read(d),read(k);
k=(k-1+e*last)%(d-c+1)+1;
c=kth(T[r],T[l-1],c);
d=kth(T[r],T[l-1],d);
last=query(root[r],root[l-1],c,d,k);
writeln(last);
}
}
}
int main(){
fread(Buf,1,BUF,stdin);
readstr(type);
read(n);
for(i=1;i<=n;i++)read(a[i]);
for(i=1;i<=n;i++)read(b[i]);
read(m),read(e);
if(type[0]=='A')SubA::Main();
else if(type[0]=='B')SubB::Main();
else{
for(i=1;i<=n;i++)T[i]=ins(T[i-1],1,n,a[i]);
if(type[0]=='C')SubC::Main();
else if(type[0]=='D')SubD::Main();
else SubF::Main();
}
fwrite(Out,1,ou-Out,stdout);
return 0;
}