【回滚莫队】三校联考 10.15 —— ants

题目描述

然而贪玩的 dirty 又开始了他的第三个游戏。
dirty 抓来了 n 只蚂蚁，并且赋予每只蚂蚁不同的编号，编号从 1 到 n。最开始，它们按某个
顺序排成一列。现在 dirty 想要进行 m 场比赛，每场比赛给出 l 和 r ，表示选出从左向右数第 l
只至第 r 只蚂蚁。被选出的蚂蚁需要快速地按编号从小到大排序，之后这些蚂蚁中编号连续的蚂
蚁将围成一个圈。每场比赛结束后，蚂蚁们还需要快速地回到最开始的位置。
按照蚂蚁的审美标准，围成的圈越大美观值就越大。于是 dirty 每次需要找到最大的圈，但由
于比赛多到难以处理，他只需要知道这个圈中蚂蚁的数目。

分析

回滚莫队应该说不是一个新的东西了,在这场考试之前已经见过一次了,但是只写了一道版题.
算法的介绍引用某大佬的博客(传送门)
当一个点被统计,则需要将与它相邻的两个数连上边,最后统计一下答案.
另外可以利用链表,比较快捷的维护出.
首先需要将之前的状态清空,然后将右端点更新至当前询问的位置,并且记录下当前的答案,接着将左区间的点加进去,统计出答案,最后撤销对于左区间的操作,将答案还原到之前的状态.这样做就免去删除操作的时间,将时间复杂度进一步降低.
详见代码.

#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;

const int MAXN=100005;


int n,m,nowans,BLOCK,top;
int ans[MAXN],bl[MAXN],br[MAXN];
int a[MAXN],stk[MAXN],tag[MAXN],T;
bool vis[MAXN];

struct node{
    int l,r,id;
    node(int _l=0,int _r=0,int _id=0):l(_l),r(_r),id(_id){}
}que[MAXN],qs[MAXN];

bool cmp_l(node a,node b){return a.l<b.l;}
bool cmp_r(node a,node b){return a.r<b.r||(a.r==b.r&&a.l<b.l);}

void update(int a){
    if(tag[a]==T) return;
    tag[a]=T,vis[a]=0,bl[a]=a,br[a]=a;
}

void restore(){
    while(top){
        int t=stk[top];
        if(t!=1&&vis[t-1]) br[bl[t-1]]=t-1;
        if(t!=n&&vis[t+1]) bl[br[t+1]]=t+1;
        vis[t]=0,top--;
    }
}

void work(int t,int f){
    update(t),update(t-1),update(t+1);
    int ll=t,rr=t;
    if(t!=1&&vis[t-1]) ll=bl[t-1];
    if(t!=n&&vis[t+1]) rr=br[t+1];
    nowans=max(nowans,rr-ll+1);
    br[ll]=rr,bl[rr]=ll,vis[t]=1;
    if(f) stk[++top]=t;
}

void sol(){
    BLOCK=(int)sqrt(n);
    sort(qs+1,qs+m+1,cmp_l);
    int now=1;
    for(int i=1;i<=n;i+=BLOCK){
        nowans=1;
        int cnt=0;
        while(now<=m&&qs[now].l>=i&&qs[now].l<=min(i+BLOCK-1,n))
            que[++cnt]=qs[now++];
        if(!cnt) continue;
        T++;
        sort(que+1,que+cnt+1,cmp_r);
        int now_r=min(n,i+BLOCK-1);
        for(int j=1;j<=cnt;j++){
            while(now_r<que[j].r) work(a[++now_r],0);
            int t_ans=nowans;
            for(int k=min(n,min(i+BLOCK-1,que[j].r));k>=que[j].l;k--) work(a[k],1);
            ans[que[j].id]=nowans; restore();
            for(int k=min(n,min(i+BLOCK-1,que[j].r));k>=que[j].l;k--) vis[a[k]]=0;
            nowans=t_ans;
        }
    }
}

int main(){
    freopen("ants.in","r",stdin);
    freopen("ants.out","w",stdout);

    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&a[i]);
    for(int i=1;i<=m;i++){
        scanf("%d%d",&qs[i].l,&qs[i].r);
        qs[i].id=i;
    }
    sol();
    for(int i=1;i<=m;i++) printf("%d\n",ans[i]);

    fclose(stdin);
    fclose(stdout);
}