4241: 历史研究

题目大意：区间内出现次数*权值最大的数。

与2724差不多：bzoj2724题解

主要check的时候返回num*a就行了。

卡时卡卡了好几次（雾）

code：

#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<cstring>
#include<vector>
#include<cmath>
#include<algorithm>
#define LL long long
using namespace std;
vector<int> v[100010];
LL n,m,e[100010],L[330],R[330],pre[100010];
struct node1{
    LL a,b,num;
}a[100010];
struct node{
    long long num;
    LL a;
}f[330][330];
bool cmp(node1 a,node1 b){return a.b<b.b;}
LL read()
{
    LL x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
long long check(LL x,LL y,LL c)
{
    LL l=0,r=v[c].size()-1,ans;
    while(l<=r)
    {
        LL mid=(l+r)>>1;
        if(v[c][mid]>=x) ans=mid,r=mid-1;
        else l=mid+1;
    }
    l=0,r=v[c].size()-1;LL ans1;
    while(l<=r)
    {
        LL mid=(l+r)>>1;
        if(v[c][mid]<=y) ans1=mid,l=mid+1;
        else r=mid-1;
    }
    return (long long)(ans1-ans+1)*e[c];
}
node findmax(LL x,LL y,LL X,LL Y)
{
    node ans;ans.num=0;
    for(LL i=x;i<=y;i++)
    {
        LL tmp=check(X,Y,a[i].a);
        if(ans.num<tmp) ans.num=tmp,ans.a=a[i].a;
    }
    return ans;
}
int main()
{
    n=read();m=read(); 
    for(LL i=1;i<=n;i++) a[i].b=read(),a[i].num=i;
    sort(a+1,a+n+1,cmp);
    long long tmp=0;
    for(LL i=1;i<=n;i++)
    {
        if(i==1||a[i].b!=a[i-1].b) tmp++;
        a[a[i].num].a=tmp;
        e[tmp]=a[i].b;
    }
    for(LL i=1;i<=tmp;i++) v[i].push_back(0);
    for(LL i=1;i<=n;i++) v[a[i].a].push_back(i);
    tmp=sqrt(n);
    for(LL i=1;i<=n;i++)
    {
        pre[i]=i/tmp+1;
        if(pre[i]!=pre[i-1]) L[pre[i]]=i,R[pre[i-1]]=i-1;
    }
    R[pre[n]]=n;
    for(LL l=1;l<=pre[n];l++)
        for(LL i=pre[1];i<=pre[n];i++)
        {
            LL j=i+l-1;if(j>pre[n]) break;
            node MAX=findmax(L[i],R[i],L[i],R[j]);
            if(i==j) f[i][j]=MAX;
            else
            {
                tmp=check(L[i],R[j],f[i+1][j].a);
                if(MAX.num>tmp) f[i][j]=MAX;
                else f[i][j].num=tmp,f[i][j].a=f[i+1][j].a;
            }
        }
    node ans;
    for(LL i=1;i<=m;i++)
    {
        LL l,r;l=read();r=read();
        if(l>r) swap(l,r);
        if(pre[l]==pre[r]) ans=findmax(l,r,l,r);
        else
        {
            LL x=R[pre[l]],y=L[pre[r]];
            node tmp1=findmax(l,x,l,r),tmp2=findmax(y,r,l,r);
            if(tmp1.num>tmp2.num) ans=tmp1;
            else ans=tmp2;
            if(pre[l]+1<=pre[r]-1)
            {
                tmp1=f[pre[l]+1][pre[r]-1];
                if(tmp1.num>ans.num) ans=tmp1;
            }
        }
        printf("%lld\n",ans.num);
    }
}