题解 AT1219 【歴史の研究】

莫队

简单分析：题面含有IOI（惊），可知此题是IOI（数字三角形）难度（逃）。

思路：回滚莫队

当然很多人都是抱着学回滚莫队的目标来看这道题的，所以这里介绍一下回滚莫队。

1、按莫队的思路讲询问排序。

2、查询时枚举每个区间，我们需要保证右端点是保持单调递增的，同时左端点每次在一个块中移动，以此来计算每个询问的值。

3、每一次到下一个块就讲左端点移回右端点，移的过程不需要再像莫队一样一个个移，只需要将莫队中改变的数据清0，然后将右端点赋值到左端点就可以了。

所以其实这道题是回滚莫队板子题。

代码

#include<stdio.h>
#include<cmath>
#include<algorithm>
#define int long long
using namespace std;
struct kkk{
    int x,y,id,flag,ans;
}a[1000001];
struct ggg{
    int a,b,p;
}a1[1000001];
int block,blo[1000001],v[1000001],vis[1000001],cnt[1000001],sum[1000001],aa[1000001];
int Vis[1000001],maxx,maxy,pos;
int n,C,m,l,r;
bool cmp1(ggg a,ggg b){
    return a.a<b.a||(a.a==b.a&&a.p<b.p);
}
bool comp(kkk a,kkk b){
    return a.id<b.id;
}
bool cmp(kkk a,kkk b){
    if(blo[a.x]!=blo[b.x])return a.x<b.x;
    return a.y<b.y;
}
void add(int x){
    vis[v[x]]++;maxx=max(maxx,vis[v[x]]*aa[x]);
}
void re(int x){vis[v[x]]--;}
int check(int l,int r){maxy=0;
    for(int i=l;i<=r;i++)Vis[v[i]]=0;
    for(int i=l;i<=r;i++){
        Vis[v[i]]++;
        maxy=max(maxy,Vis[v[i]]*aa[i]);
    }
    return maxy;
}
int Mo(int pos,int bl){
    maxx=0;int last=0,i=pos;
    for(int j=1;j<=n;j++)vis[j]=0;
    int L=min(block*bl,n);
    int l=L+1,r=L;
    for(;blo[a[i].x]==bl;i++){
        if(blo[a[i].x]==blo[a[i].y]){a[i].ans=check(a[i].x,a[i].y);continue;}
        while(r<a[i].y){add(++r);}
        last=maxx;
        while(l>a[i].x){add(--l);}
        a[i].ans=maxx;
        while(l<L+1)re(l++);
        maxx=last;
    }
    return i;
}
signed main(){
    int num;
    scanf("%lld%lld",&n,&m);block=sqrt(n);
    //离散化
    for(int i=1;i<=n;i++){scanf("%lld",&a1[i].a);aa[i]=a1[i].a;a1[i].p=i,blo[i]=(i-1)/block+1;num=blo[i];}
    sort(a1+1,a1+n+1,cmp1);
    for(int i=1,j=0;i<=n;i++)
    {
        if(i==1||a1[i].a!=a1[i-1].a)j++;
        v[a1[i].p]=j;
    }
    for(int i=1;i<=m;i++){
        scanf("%lld%lld",&a[i].x,&a[i].y);
        a[i].id=i;
    }
    pos=1;
    sort(a+1,a+m+1,cmp);
    for(int i=1;i<=num;i++){
    	pos=Mo(pos,i);
    }
    sort(a+1,a+m+1,comp);
    for(int i=1;i<=m;i++){
        printf("%lld\n",a[i].ans);
    }
    return 0;
}