数据结构（分块）：[HZOI 2015]easy seq

【题目描述】

给定一个序列，下标从0开始，分别为a0,a1,a2...an−1，有m个询问，每次给出l和r，求满足ai=aj且l<=i<=j<=r时j−i的最大值

本题强制在线，l和r均进行了加密，解密过程为

l′=min((l+ans)modn,(r+ans)modn)

r′=max((l+ans)modn,(r+ans)modn)

其中ans为上次答案

【输入格式】

第一行，两个正整数n,m

第二行，n个非负整数a0,a1,a2...an−1

以下m行，每行两个数l,r，表示询问

【输出格式】

m行，表示答案

【样例输入】

5 3
0 4 0 0 1
3 2
4 1
0 0

【样例输出】

1
2
0

0<=l,r<n

【来源】

国家集训队2015论文集

　　离线预处理+分块，复杂度有保障。

 #include <cstdio>
 const int N=,M=,B=;
 int n,Q,lmn[N][M],rmx[N][M];
 int dp[M][M],a[N],bel[N];
 int fst[N][M],lst[N][M];
 int st[M],ed[M],tot;
 int h[N],vis[N],tim;
 char c;int x;
 __inline int min(int x,int y){return x>y?y:x;}
 __inline int max(int x,int y){return x<y?y:x;}
 __inline int Read(){
     while(x=getchar()-,x>||x<);
     while(c=getchar()-,c<=&&c>=)x=x*+c;
     return x;
 }
 int main(){
     freopen("easy_seq.in","r",stdin);
     freopen("easy_seq.out","w",stdout);
     scanf("%d%d",&n,&Q);
     for(register int i=;i<=n;i++)a[i]=Read();
     for(register int i=;i<=n;i++)bel[i]=(i-)/B+;
     for(register int i=;i<=n;i++)ed[bel[i]]=i;
     for(register int i=n;i>=;i--)st[bel[i]]=i;

     tot=bel[n];

     for(register int i=n;i>=;i--)fst[a[i]][bel[i]]=i;
     for(register int i=;i<=n;i++)lst[a[i]][bel[i]]=i;

     for(register int i=;i<=n;i++){
         for(register int j=;j<=bel[i];j++)
             lmn[i][j]=fst[a[i]][j];
         for(register int j=bel[i];j<=tot;j++)
             rmx[i][j]=lst[a[i]][j];
     }

     for(register int i=;i<=n;i++){
         for(register int j=bel[i]-;j>=;j--)
             if(!lmn[i][j])lmn[i][j]=lmn[i][j+];
             else lmn[i][j]=min(lmn[i][j],lmn[i][j+]);
         for(register int j=bel[i]+;j<=tot;j++)
             rmx[i][j]=max(rmx[i][j],rmx[i][j-]);
     }

     for(register int j=,i;j<=n;j++){i=bel[j];
         dp[i][i]=max(dp[i][i],lst[a[j]][i]-j);
         dp[i][i]=max(dp[i][i],j-fst[a[j]][i]);
     }

     for(register int k=;k<=tot;k++)
         for(register int i=,j;(j=i+k-)<=tot;i++){
             dp[i][j]=max(dp[i][j-],dp[i+][j]);
             for(register int t=st[i];t<=ed[i];t++)
             dp[i][j]=max(dp[i][j],rmx[t][j]-t);
         }
     int l_,r_,l,r,ans=;
     while(Q--){
         l=Read();r=Read();
         l_=min((l+ans)%n,(r+ans)%n);
         r_=max((l+ans)%n,(r+ans)%n);
         l=l_+;r=r_+;

         if(bel[l]==bel[r]){
             ++tim;ans=;
             for(register int i=l;i<=r;i++){
                 if(vis[a[i]]!=tim){
                     vis[a[i]]=tim;
                     h[a[i]]=i;
                 }
                 else ans=max(ans,i-h[a[i]]);
             }
             printf("%d\n",ans);
             continue;
         }

         if(st[bel[l]]!=l)l_=ed[bel[l]]+;else l_=l;
         if(ed[bel[r]]!=r)r_=st[bel[r]]-;else r_=r;

         ans=dp[bel[l_]][bel[r_]];++tim;
         for(register int i=l;i<l_;i++){
             ans=max(ans,rmx[i][bel[r_]]-i);
             if(vis[a[i]]!=tim){
                 vis[a[i]]=tim;
                 h[a[i]]=i;
             }
             else ans=max(ans,i-h[a[i]]);
         }
         for(register int i=r_+;i<=r;i++){
             if(!lmn[i][bel[l_]])lmn[i][bel[l_]]=n;
             ans=max(ans,i-lmn[i][bel[l_]]);
             if(vis[a[i]]!=tim){
                 vis[a[i]]=tim;
                 h[a[i]]=i;
             }
             else ans=max(ans,i-h[a[i]]);
         }
         printf("%d\n",ans);
     }
     return ;
 }