bzoj:4241: 历史研究

Description

IOI国历史研究的第一人——JOI教授，最近获得了一份被认为是古代IOI国的住民写下的日记。JOI教授为了通过这份日记来研究古代IOI国的生活，开始着手调查日记中记载的事件。

日记中记录了连续N天发生的时间，大约每天发生一件。

事件有种类之分。第i天(1<=i<=N)发生的事件的种类用一个整数Xi表示，Xi越大，事件的规模就越大。

JOI教授决定用如下的方法分析这些日记：

1. 选择日记中连续的一些天作为分析的时间段

2. 事件种类t的重要度为t*(这段时间内重要度为t的事件数)

3. 计算出所有事件种类的重要度，输出其中的最大值

现在你被要求制作一个帮助教授分析的程序，每次给出分析的区间，你需要输出重要度的最大值。

Input

第一行两个空格分隔的整数N和Q，表示日记一共记录了N天，询问有Q次。

接下来一行N个空格分隔的整数X1...XN，Xi表示第i天发生的事件的种类

接下来Q行，第i行(1<=i<=Q)有两个空格分隔整数Ai和Bi，表示第i次询问的区间为[Ai,Bi]。

Output

输出Q行，第i行(1<=i<=Q)一个整数，表示第i次询问的最大重要度

Sample Input

5 5
9 8 7 8 9
1 2
3 4
4 4
1 4
2 4

Sample Output

9
8
8
16
16

HINT

1<=N<=10^5

1<=Q<=10^5

1<=Xi<=10^9 (1<=i<=N)

在线做法可以看czl博客

我的写法是离线莫队……说白了就是把暴力优化一下罢了……

一般离线暴力不就是枚举起点然后向右跑过去么，那么我们就分块枚举起点，在起点左边的部分就暴力搜过去。
预先离散化，开个桶然后记录值并维护max。

复杂度n*sqrt(n)

运气#1（截至2016/4/12提交时）

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int read_p,read_ca;
inline int read(){
    read_p=0;read_ca=getchar();
    while(read_ca<'0'||read_ca>'9') read_ca=getchar();
    while(read_ca>='0'&&read_ca<='9') read_p=read_p*10+read_ca-48,read_ca=getchar();
    return read_p;
}
int k;
int n,m,a[100001],pos[100001],id[100001],num=0;
long long ans[100001],tp[100001],an,qq;
struct na{
    int l,r,k,id;
}b[100001];
bool cmp(na a,na b){
    if (a.k==b.k) return a.r<b.r;else
    return a.k<b.k;
}
bool tmp(int x,int y){
    return a[x]<a[y];
}
int main(){
    n=read();m=read();
    int i,s=1,t=1,j,l,o;
    k=200;
    for (i=1;i<=n;i++) a[i]=read(),id[i]=i;
    sort(id+1,id+n+1,tmp);
    for (i=1;i<=n;i++)
    if (a[id[i]]==a[id[i-1]]) pos[id[i]]=num;else pos[id[i]]=++num;
    for (i=1;i<=m;i++) b[i].l=read(),b[i].r=read(),b[i].id=i,b[i].k=((b[i].l-1)/k+1)*k;
    sort(b+1,b+1+m,cmp);
    for (i=k;i<=n;i+=k){
        for (j=0;j<=num;j++) tp[j]=0;
        for (s=t;b[t].k==i;t++);
        l=i;an=0;
        for (;b[s].r<i&&s<t;s++){
            for (j=b[s].l;j<=b[s].r;j++)
            if ((tp[pos[j]]+=a[j])>an) an=tp[pos[j]];
            ans[b[s].id]=an;
            for (j=b[s].l;j<=b[s].r;j++)
            tp[pos[j]]-=a[j];
            an=0;
        }
        for (j=s;j<t;j++){
            for (;l<=b[j].r;l++)
            if ((tp[pos[l]]+=a[l])>an) an=tp[pos[l]];
            qq=an;
            for (o=i-1;o>=b[j].l;o--)
            if ((tp[pos[o]]+=a[o])>qq) qq=tp[pos[o]];
            for (o=b[j].l;o<i;o++) tp[pos[o]]-=a[o];
            ans[b[j].id]=qq;
        }
    }
    for (i=t;i<=m;i++){
        for (j=b[i].l;j<=b[i].r;j++)
        if ((tp[pos[j]]+=a[j])>an) an=tp[pos[j]];
        for (j=b[i].l;j<=b[i].r;j++)
        tp[pos[j]]-=a[j];
        ans[b[i].id]=an;
        an=0;
    }
    for (i=1;i<=m;i++) printf("%lld
",ans[i]);
}

View Code