历史研究
题解
莫队水题。
一看是序列的查询,就知道又要用那些数据结构。
好吧,题目说得很明确了,我们很明显可以用莫队来维护。
可是,如果直接莫队的话应该只有5pts。
忽然发现序列长度不超过,那么就离散化一遍,貌似可以过15pts了。
我们发现我们这样的算法是的,我们还可以在最后统计答案的地方加上个线段树,中间维护那些重要度时用的线段树来维护,这样就可以卡到,太难看了,不过因为常数的原因,只有40pts。
我们发现,这个时间复杂度理论上是过得去的,于是用各种奇技淫巧来卡常,然后就卡过了。。。
其实,还可以用回滚莫队将那一个log去掉,这样就可以以优美复杂度过掉。
源码
带log
#pragma GCC optimize(2)
#pragma GCC optimize(3)
#pragma GCC optimize("Ofast")
#pragma GCC optimize("inline")
#pragma GCC optimize("-fgcse")
#pragma GCC optimize("-fgcse-lm")
#pragma GCC optimize("-fipa-sra")
#pragma GCC optimize("-ftree-pre")
#pragma GCC optimize("-ftree-vrp")
#pragma GCC optimize("-fpeephole2")
#pragma GCC optimize("-ffast-math")
#pragma GCC optimize("-fsched-spec")
#pragma GCC optimize("unroll-loops")
#pragma GCC optimize("-falign-jumps")
#pragma GCC optimize("-falign-loops")
#pragma GCC optimize("-falign-labels")
#pragma GCC optimize("-fdevirtualize")
#pragma GCC optimize("-fcaller-saves")
#pragma GCC optimize("-fcrossjumping")
#pragma GCC optimize("-fthread-jumps")
#pragma GCC optimize("-funroll-loops")
#pragma GCC optimize("-fwhole-program")
#pragma GCC optimize("-freorder-blocks")
#pragma GCC optimize("-fschedule-insns")
#pragma GCC optimize("inline-functions")
#pragma GCC optimize("-ftree-tail-merge")
#pragma GCC optimize("-fschedule-insns2")
#pragma GCC optimize("-fstrict-aliasing")
#pragma GCC optimize("-fstrict-overflow")
#pragma GCC optimize("-falign-functions")
#pragma GCC optimize("-fcse-skip-blocks")
#pragma GCC optimize("-fcse-follow-jumps")
#pragma GCC optimize("-fsched-interblock")
#pragma GCC optimize("-fpartial-inlining")
#pragma GCC optimize("no-stack-protector")
#pragma GCC optimize("-freorder-functions")
#pragma GCC optimize("-findirect-inlining")
#pragma GCC optimize("-fhoist-adjacent-loads")
#pragma GCC optimize("-frerun-cse-after-loop")
#pragma GCC optimize("inline-small-functions")
#pragma GCC optimize("-finline-small-functions")
#pragma GCC optimize("-ftree-switch-conversion")
#pragma GCC optimize("-foptimize-sibling-calls")
#pragma GCC optimize("-fexpensive-optimizations")
#pragma GCC optimize("-funsafe-loop-optimizations")
#pragma GCC optimize("inline-functions-called-once")
#pragma GCC optimize("-fdelete-null-pointer-checks")
#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
#include<vector>
#include<map>
using namespace std;
#define MAXN 100005
typedef long long LL;
typedef pair<LL,LL> pll;
#define gc() getchar()
template<typename _T>
inline void read(_T &x){
_T f=1;x=0;char s=gc();
while(s>'9'||s<'0'){if(s=='-')f=-1;s=gc();}
while(s>='0'&&s<='9'){x=(x<<3)+(x<<1)+(s^48);s=gc();}
x*=f;
}
int sn;
struct ming{
int l,r,id;
}t[MAXN];
bool cmp(const ming &x,const ming &y){
return (x.l/sn)==(y.l/sn)?x.r<y.r:x.l<y.l;
}
int n,q,x[MAXN],tot,idx[MAXN],L,R;
LL maxx[MAXN<<2],ans[MAXN];
map<int,int> mp;
LL Max(LL a,LL b){
return a>b?a:b;
}
void insert(int rt,int l,int r,int ai,int aw){
if(l==r){
maxx[rt]+=1ll*aw;
return ;
}
int mid=(l+r)>>1;
if(ai<=mid)insert(rt<<1,l,mid,ai,aw);
if(ai>mid)insert(rt<<1|1,mid+1,r,ai,aw);
maxx[rt]=Max(maxx[rt<<1],maxx[rt<<1|1]);
}
signed main(){
read(n);read(q);sn=sqrt(n);
for(int i=1;i<=n;++i){
read(x[i]);
if(!mp[x[i]])mp[x[i]]=++tot,idx[tot]=x[i];
x[i]=mp[x[i]];
}
for(int i=1;i<=q;++i)read(t[i].l),read(t[i].r),t[i].id=i;
sort(t+1,t+q+1,cmp);L=1;R=0;
for(int i=1;i<=q;++i){
int l=t[i].l,r=t[i].r;
while(L>l)--L,insert(1,1,n,x[L],idx[x[L]]);
while(R<r)++R,insert(1,1,n,x[R],idx[x[R]]);
while(L<l)insert(1,1,n,x[L],-idx[x[L]]),++L;
while(R>r)insert(1,1,n,x[R],-idx[x[R]]),--R;
ans[t[i].id]=maxx[1];
}
for(int i=1;i<=q;++i)printf("%lld\n",ans[i]);
return 0;
}
不带log
#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
#include<vector>
#include<map>
using namespace std;
#define MAXN 100005
typedef long long LL;
#define int LL
typedef pair<LL,LL> pll;
#define gc() getchar()
template<typename _T>
inline void read(_T &x){
_T f=1;x=0;char s=gc();
while(s>'9'||s<'0'){if(s=='-')f=-1;s=gc();}
while(s>='0'&&s<='9'){x=(x<<3)+(x<<1)+(s^48);s=gc();}
x*=f;
}
struct ming{
int l,r,b,id;
inline bool operator < (const ming &c){
return b!=c.b?l<c.l:r<c.r;
}
}t[MAXN];
int a[MAXN],x[MAXN],an,sn,ans[MAXN];
int sum[MAXN],n,q;
signed main(){
read(n);read(q);sn=sqrt(n+1);
for(int i=1;i<=n;i++)read(x[i]),a[i]=x[i];
sort(a+1,a+n+1);an=unique(a+1,a+n+1)-a-1;
for(int i=1;i<=n;i++)x[i]=lower_bound(a+1,a+an+1,x[i])-a;
for(int i=1;i<=q;i++){
read(t[i].l);read(t[i].r);
t[i].b=(t[i].l-1)/sn+1;t[i].id=i;
if(t[i].b*sn>=t[i].r){
for(int j=t[i].l;j<=t[i].r;j++)ans[i]=max(ans[i],++sum[x[j]]*a[x[j]]);
for(int j=t[i].l;j<=t[i].r;j++)--sum[x[j]];
}
}
sort(t+1,t+q+1);
for(int L=1,R=1;R<=q;L=R+1){
memset(sum,0,sizeof(sum));
for(R=L;R<q&&t[R+1].b==t[L].b;++R);
int num=0;
for(int i=L,p=t[L].b*sn,r=p;i<=R;i++){
if(t[i].r<=p)continue;
while(r<t[i].r)++r,num=max(num,++sum[x[r]]*a[x[r]]);
int tmp=num;
for(int j=t[i].l;j<=p;j++)tmp=max(tmp,++sum[x[j]]*a[x[j]]);
ans[t[i].id]=tmp;
for(int j=t[i].l;j<=p;j++)--sum[x[j]];
}
}
for(int i=1;i<=q;i++)printf("%lld\n",ans[i]);
return 0;
}