LOJ2874 JOISC2014 历史研究 分块、莫队

传送门

---

看到出现次数自然地考虑莫队。

但是发现如果需要删除并动态维护答案的话，则要用一个堆来维护答案，增加了一个(log)。但是加入操作却没有这个(log)，所以我们考虑避免删除操作。

分块，设(l_i,r_i)表示第(i)个块的左右端点，设(f_{i,j})表示区间([l_i,r_j])的答案，可以枚举(i)然后枚举(j)做到(O(nsqrt{n}))；

接下来将询问离线，对于每一组询问如果左右端点距离(leq 2sqrt{n})则暴力计算答案，否则考虑其覆盖的整块([p,q])，则通过莫队将([l_p,r_q])内所有数的出现次数统计下来，然后暴力将零散块的贡献进入出现次数数组并更新当前答案，最后把出现次数数组还原即可。复杂度(O(nsqrt{n}))。

#include<bits/stdc++.h>
using namespace std;

int read(){
	int a = 0; char c = getchar(); bool f = 0;
	while(!isdigit(c)){f = c == '-'; c = getchar();}
	while(isdigit(c)){
		a = a * 10 + c - 48; c = getchar();
	}
	return f ? -a : a;
}

const int _ = 1e5 + 7 , S = sqrt(_) , T = _ / S + 10;
long long mx[T][T] , ans[_]; int arr[_] , tmp[_] , lsh[_] , N , Q , num;
struct query{
	int l , r , id , tL , tR;
	query(int _l , int _r , int _id){
		l = _l; r = _r; id = _id; tL = (l / S + 1) * S; tR = (r / S) * S - 1;
	}
	friend bool operator <(query A , query B){
		return A.tL < B.tL || A.tL == B.tL && A.tR < B.tR;
	}
}; vector < query > qry;

void add(int x){++tmp[arr[x]];} void del(int x){--tmp[arr[x]];}

signed main(){
#ifndef ONLINE_JUDGE
	freopen("in","r",stdin);
	freopen("out","w",stdout);
#endif
	N = read(); Q = read();
	for(int i = 0 ; i < N ; ++i) arr[i] = lsh[i] = read();
	sort(lsh , lsh + N); num = unique(lsh , lsh + N) - lsh - 1;
	for(int i = 0 ; i < N ; ++i) arr[i] = lower_bound(lsh , lsh + num + 1 , arr[i]) - lsh;

	for(int i = 0 ; i < N ; i += S){
		long long now = 0;
		for(int j = i ; j < N ; ++j){
			now = max(now , 1ll * ++tmp[arr[j]] * lsh[arr[j]]);
			if(j % S == S - 1 || j + 1 == N) mx[i / S][j / S] = now;
		}
		memset(tmp , 0 , sizeof(tmp));
	}
	
	for(int i = 1 ; i <= Q ; ++i){
		int p = read() - 1 , q = read() - 1; long long now = 0;
		if(q - p <= 2 * S){
			for(int k = p ; k <= q ; ++k)
				now = max(now , 1ll * ++tmp[arr[k]] * lsh[arr[k]]);
			for(int k = p ; k <= q ; ++k) tmp[arr[k]] = 0;
			ans[i] = now;
		}
		else qry.push_back(query(p , q , i));
	}

	sort(qry.begin() , qry.end()); int L = 0 , R = -1;
	for(auto t : qry){
		while(R < t.tR) add(++R); while(L > t.tL) add(--L);
		while(R > t.tR) del(R--); while(L < t.tL) del(L++);
		long long now = mx[L / S][R / S];
		for(int i = L - 1 ; i >= t.l ; --i)
			now = max(now , 1ll * ++tmp[arr[i]] * lsh[arr[i]]);
		for(int i = R + 1 ; i <= t.r ; ++i)
			now = max(now , 1ll * ++tmp[arr[i]] * lsh[arr[i]]);
		ans[t.id] = now;
		for(int i = t.l ; i < L ; ++i) --tmp[arr[i]];
		for(int i = t.r ; i > R ; --i) --tmp[arr[i]];
	}
	for(int i = 1 ; i <= Q ; ++i) printf("%lld
" , ans[i]);
	return 0;
}