洛谷 P4168 [Violet]蒲公英

https://www.luogu.org/problemnew/show/P4168

分块大暴力

错误记录：66,82行lans=ans.b

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<tr1/unordered_map>
using namespace std;
using namespace tr1;
struct D
{
    int a,b;
    D(){}
    D(int x,int y):a(x),b(y){}
};
bool operator<(const D &a,const D &b)    {return a.a<b.a||(a.a==b.a&&a.b>b.b);}
int n,c,m,sz=200,sz1;
int a[40100];
int be[40100],st[202],ed[202];
int nnm[202][40100];
int ttt[40100];
D tt2[202][202];
unordered_map<int,int> ma;
int tt[40100];
int main()
{
    int i,j,k,l,r,lans=0;D ans;
    scanf("%d%d",&n,&m);
    //sz=sqrt(n);
    sz1=(n-1)/sz+1;
    for(i=1;i<=n;i++)    scanf("%d",&a[i]),be[i]=(i-1)/sz+1,tt[++tt[0]]=a[i];
    for(i=1;i<sz1;i++)    st[i]=(i-1)*sz+1,ed[i]=i*sz;
    st[sz1]=(sz1-1)*sz+1,ed[sz1]=n;
    sort(tt+1,tt+tt[0]+1);tt[0]=unique(tt+1,tt+tt[0]+1)-tt-1;
    for(i=1;i<=tt[0];i++)    ma[tt[i]]=i;
    for(i=1;i<=n;i++)    a[i]=ma[a[i]];
    for(i=1;i<=sz1;i++)
    {
        memcpy(nnm[i],nnm[i-1],sizeof(nnm[i]));
        for(j=st[i];j<=ed[i];j++)    nnm[i][a[j]]++;
    }
    for(i=1;i<=sz1;i++)
    {
        memset(ttt,0,sizeof(ttt));ans=D(0,0);
        for(j=i;j<=sz1;j++)
        {
            for(k=st[j];k<=ed[j];k++)
            {
                ttt[a[k]]++;
                ans=max(ans,D(ttt[a[k]],a[k]));
            }
            tt2[i][j]=ans;
        }
    }
    memset(ttt,0,sizeof(ttt));
    while(m--)
    {
        scanf("%d%d",&l,&r);l=(l+lans-1)%n+1;r=(r+lans-1)%n+1;
        if(l>r)    swap(l,r);
        if(be[l]+1>=be[r])
        {
            ans=D(0,0);
            for(i=l;i<=r;i++)
            {
                ttt[a[i]]++;
                ans=max(ans,D(ttt[a[i]],a[i]));
            }
            printf("%d
",tt[ans.b]);lans=tt[ans.b];
            for(i=l;i<=r;i++)    ttt[a[i]]--;
        }
        else
        {
            ans=tt2[be[l]+1][be[r]-1];
            for(i=l;i<=ed[be[l]];i++)
            {
                ttt[a[i]]++;
                ans=max(ans,D(ttt[a[i]]+nnm[be[r]-1][a[i]]-nnm[be[l]][a[i]],a[i]));
            }
            for(i=st[be[r]];i<=r;i++)
            {
                ttt[a[i]]++;
                ans=max(ans,D(ttt[a[i]]+nnm[be[r]-1][a[i]]-nnm[be[l]][a[i]],a[i]));
            }
            printf("%d
",tt[ans.b]);lans=tt[ans.b];
            for(i=l;i<=ed[be[l]];i++)    ttt[a[i]]--;
            for(i=st[be[r]];i<=r;i++)    ttt[a[i]]--;
        }
    }
    return 0;
}

---

额外：

loj #6285. 数列分块入门 9

https://loj.ac/problem/6285

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<tr1/unordered_map>
using namespace std;
using namespace tr1;
struct D
{
    int a,b;
    D(){}
    D(int x,int y):a(x),b(y){}
};
bool operator<(const D &a,const D &b)    {return a.a<b.a||(a.a==b.a&&a.b>b.b);}
int n,c,m,sz=300,sz1;
int a[100100];
int be[100100],st[340],ed[340];
int nnm[340][100100];
int ttt[100100];
D tt2[340][340];
unordered_map<int,int> ma;
int tt[100100];
int main()
{
    int i,j,k,l,r,lans=0;D ans;
    scanf("%d",&n);m=n;
    //sz=sqrt(n);
    sz1=(n-1)/sz+1;
    for(i=1;i<=n;i++)    scanf("%d",&a[i]),be[i]=(i-1)/sz+1,tt[++tt[0]]=a[i];
    for(i=1;i<sz1;i++)    st[i]=(i-1)*sz+1,ed[i]=i*sz;
    st[sz1]=(sz1-1)*sz+1,ed[sz1]=n;
    sort(tt+1,tt+tt[0]+1);tt[0]=unique(tt+1,tt+tt[0]+1)-tt-1;
    for(i=1;i<=tt[0];i++)    ma[tt[i]]=i;
    for(i=1;i<=n;i++)    a[i]=ma[a[i]];
    for(i=1;i<=sz1;i++)
    {
        memcpy(nnm[i],nnm[i-1],sizeof(nnm[i]));
        for(j=st[i];j<=ed[i];j++)    nnm[i][a[j]]++;
    }
    for(i=1;i<=sz1;i++)
    {
        memset(ttt,0,sizeof(ttt));ans=D(0,0);
        for(j=i;j<=sz1;j++)
        {
            for(k=st[j];k<=ed[j];k++)
            {
                ttt[a[k]]++;
                ans=max(ans,D(ttt[a[k]],a[k]));
            }
            tt2[i][j]=ans;
        }
    }
    memset(ttt,0,sizeof(ttt));
    while(m--)
    {
        scanf("%d%d",&l,&r);//l=(l+lans-1)%n+1;r=(r+lans-1)%n+1;
        if(l>r)    swap(l,r);
        if(be[l]+1>=be[r])
        {
            ans=D(0,0);
            for(i=l;i<=r;i++)
            {
                ttt[a[i]]++;
                ans=max(ans,D(ttt[a[i]],a[i]));
            }
            printf("%d
",tt[ans.b]);lans=tt[ans.b];
            for(i=l;i<=r;i++)    ttt[a[i]]--;
        }
        else
        {
            ans=tt2[be[l]+1][be[r]-1];
            for(i=l;i<=ed[be[l]];i++)
            {
                ttt[a[i]]++;
                ans=max(ans,D(ttt[a[i]]+nnm[be[r]-1][a[i]]-nnm[be[l]][a[i]],a[i]));
            }
            for(i=st[be[r]];i<=r;i++)
            {
                ttt[a[i]]++;
                ans=max(ans,D(ttt[a[i]]+nnm[be[r]-1][a[i]]-nnm[be[l]][a[i]],a[i]));
            }
            printf("%d
",tt[ans.b]);lans=tt[ans.b];
            for(i=l;i<=ed[be[l]];i++)    ttt[a[i]]--;
            for(i=st[be[r]];i<=r;i++)    ttt[a[i]]--;
        }
    }
    return 0;
}

---