BZOJ4241历史研究题解

题目连接

很显然可以想到分块，用f[i][j]表示块i到块j的ans，然后发现答案一定是f[i][j]

或者其他在边角出现的数字

我们在记下g[i][j]从开头到块i中的数字j出现的次数

这样就每一次就统计边角数字出现的次数，然后更新答案就好了

（好像莫对也可以做。。。。）

注意每一次查询时不要memset，会T飞

# include<iostream>
# include<cstdio>
# include<cstring>
# include<algorithm>
# include<cmath>
using namespace std;
typedef long long LL;
const int mn = ;
int n,m,siz;
int bl[mn],b[mn],c[mn],a[mn];
LL f[][];//块i到j的答案
int g[][mn];//从1到块i中每个数出现的次数
int tmp[mn],st[mn],top;
void pre(int x)
{
    LL ans=;
    memset(tmp,,sizeof(tmp));
    for(int i=;i<=n;i++)
        g[x][i]=g[x-][i];
    int y=min(n,x*siz);
    for(int i=(x-)*siz+;i<=y;i++)
        g[x][a[i]]++;
    for(int i=(x-)*siz+;i<=n;i++)
    {
        tmp[a[i]]++;
        if(1ll*b[a[i]]*tmp[a[i]]>ans) ans=1ll*b[a[i]]*tmp[a[i]];
        if(i%siz==)
            f[x][i/siz]=ans;
        if(i==n)
            f[x][bl[n]]=ans;
    }
}
LL ask(int x,int y)
{
    LL ans=;
    top=;
    //memset(tmp,0,sizeof(tmp));
    if(bl[x]==bl[y])
    {
        for(int i=x;i<=y;i++)
        {
            if(tmp[a[i]]==) st[++top]=a[i];
            tmp[a[i]]++;
            if(1ll*b[a[i]]*tmp[a[i]]>ans) ans=1ll*b[a[i]]*tmp[a[i]];
        }
        for(int i=top;i>=;i--)
            tmp[st[i]]=;
        return ans;
    }
    else {
        if(bl[y]->=bl[x]+) ans=f[bl[x]+][bl[y]-];
        for(int i=x;i<=bl[x]*siz;i++)
        {
            if(tmp[a[i]]==) st[++top]=a[i],tmp[a[i]]=g[bl[y]-][a[i]]-g[bl[x]][a[i]];
            tmp[a[i]]++;
            if(1ll*b[a[i]]*tmp[a[i]]>ans) ans=1ll*b[a[i]]*tmp[a[i]];
        }
        for(int i=(bl[y]-)*siz+;i<=y;i++)
        {
            if(tmp[a[i]]==) st[++top]=a[i],tmp[a[i]]=g[bl[y]-][a[i]]-g[bl[x]][a[i]];
            tmp[a[i]]++;
            if(1ll*b[a[i]]*tmp[a[i]]>ans) ans=1ll*b[a[i]]*tmp[a[i]];
        }
        for(int i=top;i>=;i--)
            tmp[st[i]]=;
        return ans;
    }
}
int main()
{
    //freopen("4241.in","r",stdin);
    //freopen("4241.out","w",stdout);
    int x,y;
    scanf("%d%d",&n,&m);
    siz=sqrt(n*1.0);
    for(int i=;i<=n;i++)
        scanf("%d",&a[i]),b[i]=a[i];
    sort(b+,b++n);
    for(int i=;i<=n;i++)
        a[i]=lower_bound(b+,b++n,a[i])-b;
    /*for(int i=1;i<=n;i++)
        printf("%d ",a[i]);离散化*/
    for(int i=;i<=n;i++)
        bl[i]=(i-)/siz+;
    for(int i=;i<=bl[n];i++)
        pre(i);
    memset(tmp,,sizeof(tmp));
    for(int i=;i<=m;i++)
    {
        scanf("%d%d",&x,&y);
        printf("%lld\n",ask(x,y));
    }
    return ;
}