[JOISC 2014 Day1]歷史研究

歷史研究

題解

莫隊水題。

一看是序列的查詢，就知道又要用那些數據結構。

好吧，題目說得很明確了，我們很明顯可以用莫隊來維護。

可是，如果直接莫隊的話應該只有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;
}

謝謝！！！