bzoj4241 历史研究

　　莫队算法,只考虑对一个序列加入数字的话，重要度的维护是O(1）的，排完序后的询问，若左端点是在同一块中的话，右端点递增。因为右端点是递增的因此可以O（1）维护，而左端点的话，对于每个询问我们可以暴力插入，询问后还原，对于一次询问来说是O（sqrtn)的

，因此总复杂度O(nsqrtn）

#include<cstdio>
#include<algorithm>
#include<map>
#define N 200010
using namespace std;
int n,m,i,a[N],b[N],c[N],tot,j,k,o;
long long ans,tmp,Ans[N],sum[N];
int cur,lim;
map<int,int> id;
struct g{
    int l,r,id;
}q[N];
bool cmp(g a,g b)
{
    if (a.l/300==b.l/300)
    return a.r<b.r;
    return a.l/300<b.l/300;
}
int main()
{
    scanf("%d%d",&n,&m);
    for (i=1;i<=n;i++)
    scanf("%d",&a[i]),b[i]=a[i]; 
    sort(a+1,a+1+n);
    for (i=1;i<=n;i++)
    {
        if (a[i]!=a[i-1]) tot++;
        id[a[i]]=tot;
    }
    for (i=1;i<=n;i++)
    c[i]=id[b[i]];
    for (i=1;i<=m;i++)
    {
        scanf("%d%d",&q[i].l,&q[i].r);
        q[i].id=i;
    }

    sort(q+1,q+1+m,cmp);
    for (i=1;i<=m;i++)
    {
        for (j=i;j<=m;j++)
        if (q[j].l/300!=q[i].l/300) break;
        k=j-1;cur=(q[i].l/300+1)*300;lim=cur-1;
        ans=0;
        for (j=1;j<=tot;j++)
        sum[j]=0;
        for (j=i;j<=k;j++)
        {
            while (cur<=q[j].r) 
            {
                sum[c[cur]]++;
                ans=max(ans,sum[c[cur]]*b[cur]);
                cur++;
            } 
            tmp=ans;
            for (o=q[j].l;o<=min(q[j].r,lim);o++)
            {
                sum[c[o]]++;
                tmp=max(tmp,sum[c[o]]*b[o]);
            }
            Ans[q[j].id]=tmp;
            for (o=q[j].l;o<=min(q[j].r,lim);o++)
                sum[c[o]]--;
        }
        i=k;
    }
    
    for (i=1;i<=m;i++)
    printf("%lld
",Ans[i]);
}