Luogu P5048 [Ynoi2019模拟赛]Yuno loves sqrt technology III 分块

这才是真正的$N\sqrt{N}$吧$qwq$

---

记录每个数$vl$出现的位置$s[vl]$，和每个数$a[i]=vl$是第几个$vl$,记为$P[i]$，然后预处理出块$[i,j]$区间的答案$f[i][j]$;

对于$[l,r]$,现将$ans$设为$[l,r]$中整块的答案；对于散块，将散块中的每个数$a[i]=vl$,有对应的$P[i]$，我们用$s[vl]$检查第$P[i]+ans+1$个$vl$对应的位置是否在$[l,r]$，如果在的话就更新答案$ans=ans+1$，并重复上述过程，直到$P[i]+ans+1$不在$[l,r]$中。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<cctype>
#include<cstdlib>
#include<vector>
#include<queue>
#include<map>
#include<set>
#define ull unsigned long long
#define ll long long
#define R register int
#define pause (for(R i=1;i<=10000000000;++i))
#define OUT freopen("out.out","w",stdout);
using namespace std;
namespace Fread {
    static char B[<<],*S=B,*D=B;
    #define getchar() (S==D&&(D=(S=B)+fread(B,1,1<<15,stdin),S==D)?EOF:*S++)
    inline int g() {
        R ret=,fix=; register char ch; while(!isdigit(ch=getchar())) fix=ch=='-'?-:fix;
        do ret=ret*+(ch^); while(isdigit(ch=getchar())); return ret*fix;
    } inline bool isempty(const char& ch) {return ch<=||ch>=;}
    inline void gs(char* s) {register char ch; while(isempty(ch=getchar())); do *s++=ch; while(!isempty(ch=getchar()));}
}using Fread::g; using Fread::gs;
const int N=,M=;
int n,m,T,tot,lst; int cnt[N],pos[N],P[N],a[N],b[N],f[M][M];
vector <int> s[N];
#define pb push_back
signed main() {
    n=g(),m=g(),T=pow(n,/2.0); for(R i=;i<=n;++i) a[i]=g(); memcpy(b,a,sizeof(int)*(n+));
    sort(b+,b+n+); tot=unique(b+,b+n+)-b-;
    for(R i=;i<=n;++i) a[i]=lower_bound(b+,b+tot+,a[i])-b;
    for(R i=;i<=n;++i) s[a[i]].pb(i);
    for(R i=;i<=n;++i) P[i]=++cnt[a[i]]-; memset(cnt,,sizeof(cnt));
    for(R i=;i<=n;++i) pos[i]=(i-)/T+;
    for(R p=,L=pos[n],mx=;p<=L;++p,mx=,memset(cnt,,sizeof(cnt))) for(R t=p;t<=L;++t) {
        for(R i=(t-)*T+,LL=min(t*T,n);i<=LL;++i) mx=max(++cnt[a[i]],mx);
        f[p][t]=mx;
    } for(R i=,l,r;i<=m;++i) { R ans=;
        l=g()^lst,r=g()^lst;
        if(pos[l]==pos[r]) {
            for(R i=l;i<=r;++i) ans=max(ans,++cnt[a[i]]); for(R i=l;i<=r;++i) cnt[a[i]]=;
        } else { R p=pos[l]+,q=pos[r]-; ans=f[p][q];
            for(R i=l,L=pos[l]*T;i<=L;++i) { R tmp=P[i];
                while((tmp+ans)<s[a[i]].size()&&s[a[i]][tmp+ans]<=r) ++ans;
            } for(R i=q*T+;i<=r;++i) { R tmp=P[i];
                while((tmp-ans)>=&&s[a[i]][tmp-ans]>=l) ++ans;
            }
        } printf("%d\n",lst=ans);
    }
}

---

2019.07.03