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)

[b]在线做法可以看czl博客[/b]
[b]我的写法是离线莫队……说白了就是把暴力优化一下罢了……[/b]
[b]一般离线暴力不就是枚举起点然后向右跑过去么，那么我们就[/b][b]分块枚举起点，在起点左边的部分就暴力搜过去。
预先离散化，开个桶然后记录值并维护max。[/b]
[b]复杂度n*sqrt(n)[/b]

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

#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\n",ans[i]);
}

View Code