[JOISC 2014 Day1]历史研究

历史研究

题解

莫队水题。

一看是序列的查询,就知道又要用那些数据结构。

好吧,题目说得很明确了,我们很明显可以用莫队来维护。

可是a_{i}\leq 10^{9},如果直接莫队的话应该只有5pts。

忽然发现序列长度不超过10^{5},那么就离散化一遍,貌似可以过15pts了。

我们发现我们这样的算法是O\left(n^{2}\sqrt{n} \right )的,我们还可以在最后统计答案的地方加上个线段树,中间维护那些重要度时用O\left(log_{n} \right )的线段树来维护,这样就可以卡到O\left(n\sqrt{n}log_{n} \right )太难看了,不过因为常数的原因,只有40pts。

我们发现,这个时间复杂度理论上是过得去的,于是用各种奇技淫巧来卡常,然后就卡过了。。。

其实,还可以用回滚莫队将那一个log去掉,这样就可以以O\left(n\sqrt{n} \right )优美复杂度过掉。

源码

带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;
}

谢谢!!!

发布了117 篇原创文章 · 获赞 154 · 访问量 1万+

猜你喜欢

转载自blog.csdn.net/Tan_tan_tann/article/details/103832732
今日推荐