bzoj4534: 基础排序算法练习题

传送门

 

 

策爷的论文题啊……题解在这儿

我只想知道为什么这题的弱化版会出现在我们今天的%你赛里……

题意：给你一堆操作$(l,r)$，表示将区间$(l,r)$按升序排序。以及$q$个询问，每次询问一个数列能否在上述操作之后变为有序

首先，我们要知道如果只有一次询问的话，我们要怎么乱搞。首先，排序的交换次数下界等于逆序对个数，我们只要对于$(l,r)$中的每一个满足$a[i]>a[i+1]$的元素，交换$i$和$i+1$，那么交换次数就能刚好为逆序对个数，即$O(n^2)$

所以对于每一个操作，我们都可以在区间中消逆序对。考虑怎么快速找到所有的逆序对。我们可以把所有的满足$a[i]>a[i+1]$的$i$给丢进set里，那么每次只要二分找到在这个区间中的$i$，然后交换，删掉原来的逆序对并加入新的逆序对就可以了。时间复杂度为$O((n^2+m)logn)$，优于$O(nmlogn)$代码如下

void SORT(rint l,rint r){
    rint i;
    while((i=*s.upper_bound(l-))>=l&&i<r){
        s.erase(i),swap(a[i],a[i+]);
        if(a[i]<a[i+]&&a[i+]>a[i+])s.insert(i+);
        if(a[i-]>a[i]&&a[i-]<a[i+])s.insert(i-);
    }
}

然后就是神仙操作了……我们先考虑什么样的初始序列可以被排好序，我们从$01$序列开始考虑，如果一个$111000$的初始序列可以通过一堆操作之后变得有序，那么$101010$的序列必然也可以（感性理解一下发现很显然），那么我们就说$101010$比$111000$更优。总的来说，一个序列中$0$越靠左，$1$越靠右，那么这个序列越接近升序，这个序列就越优。设$pos(a,i)$表示数列$a$中第$i$个$1$所在的位置。对于两个数列来说，如果$a$比$b$更优，当且仅当对于所有$i$都有$pos(a,i)\geq pos(b,i)$。

那么我们只要能得到最差的能被排好序的初始序列，那么所有比它更优的一定能够排好序。

于是我们可以设原数列为$A$，且是一个$n$的排列（如果不是的话我们可以把它给离散，并把元素相同的按下表排列，可以发现这样并不会影响结果）。考虑如何从$01$序列拓展到这个$n$排列，对于$A$和某个数值$k$，定义$01$序列$B_k$，$B_k[i]$为$1$当且仅当$A[i]\geq k$。

于是我们发现，如果对于所有的$k=2,3,...n$，$B_k$都能在这些操作下变得有序，那么原排列也可以

直接做太慢了，我们设原排列为$1,2,...,n$，然后对于$m$次操作，我们按$i$从$m$到$1$的顺序每一次都用上面的方法使区间$[l_i,r_i]$降序排序（变为最差的）

最后得到的序列中，小于$k$的位置用$0$表示，其他位置用$1$表示，得到的$01$序列为$C_k$，且$C_k$只要更改一个位置上的值就能变为$C_{k-1}$

于是只要所有的$B_k$都比对应的$C_k$优秀，那么原排列就能在若干次操作后变为有序，即$B_k$的任意一个前缀和都不能大于$C_k$的前缀和

如果我们把$B_k$中的$1$取相反数，并把两个前缀和相加，那么只要每一个前缀和都大于等于$0$即可。设$posb[i]$表示元素$i$在$B_k$中的位置，$posc[i]$表示元素$i$在$C_k$中的位置。当求完$k=i$要求$k=i+1$的时候，$posb[i+1]$处的元素变为$-1$，那么它以及后面所有数的前缀和都会减一，$posc[i]$的变化同理。于是我们只要维护好前缀和，并判断最小的前缀和是否小于$0$即可。区间赋值和区间取$min$用线段树就可以了

 //minamoto
 #include<bits/stdc++.h>
 #define rint register int
 #define min(x,y) ((x)<(y)?(x):(y))
 #define ls (p<<1)
 #define rs (p<<1|1)
 #define ppd(p,v) (mn[p]+=(v),tag[p]+=(v))
 #define mem(a) (memset(a,0,sizeof(a)))
 using namespace std;
 #define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
 char buf[<<],*p1=buf,*p2=buf;
 int read(){
     int res,f=;char ch;
     while((ch=getc())>''||ch<'')(ch=='-')&&(f=-);
     for(res=ch-'';(ch=getc())>=''&&ch<='';res=res*+ch-'');
     return res;
 }
 const int N=,M=1e6+;
 int n,m,q,c[N],posc[N],b[N],posb[N],mn[N<<],tag[N<<];
 struct ques{int l,r;}p[M];set<int>s;
 struct node{
     int val,id;
     inline bool operator <(const node &b)const{return val==b.val?id<b.id:val<b.val;}
 }f[N];
 void SORT(rint l,rint r){
     rint i;
     while((i=*s.upper_bound(l-))>=l&&i<r){
         s.erase(i),swap(c[i],c[i+]);
         if(c[i+]<c[i+]&&c[i]>c[i+])s.insert(i+);
         if(c[i-]<c[i]&&c[i-]>c[i+])s.insert(i-);
     }
 }
 inline void pd(int p){if(tag[p])ppd(ls,tag[p]),ppd(rs,tag[p]),tag[p]=;}
 void update(int p,int l,int r,int ql,int qr,int v){
     if(ql<=l&&qr>=r)return (void)(ppd(p,v));
     int mid=(l+r)>>;pd(p);
     if(ql<=mid)update(ls,l,mid,ql,qr,v);
     if(qr>mid)update(rs,mid+,r,ql,qr,v);
     mn[p]=min(mn[ls],mn[rs]);
 }
 int main(){
 //    freopen("testdata.in","r",stdin);
     n=read(),m=read(),q=read();
     for(int i=;i<=m;p[i].l=read(),p[i].r=read(),++i);
     for(int i=;i<=n;c[i]=i,++i);
     for(int i=;i<n;s.insert(i),++i);
     for(int i=m;i;SORT(p[i].l,p[i].r),--i);
     for(int i=;i<=n;posc[c[i]]=i,++i);
     while(q--){
         for(int i=;i<=n;f[i].val=read(),f[i].id=i,++i);
         sort(f+,f++n);
         for(int i=;i<=n;b[f[i].id]=i,++i);
         for(int i=;i<=n;posb[b[i]]=i,++i);
         mem(mn),mem(tag);
         bool flag=;
         for(int i=;i<=n;++i){
             update(,,n,posb[i-],n,);
             update(,,n,posc[i-],n,-);
             if(mn[]<){flag=;break;}
         }
         puts(flag?"AC":"WA");
     }
     return ;
 }