博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
LOJ2874 JOISC2014 历史研究 分块、莫队
阅读量:4656 次
发布时间:2019-06-09

本文共 2787 字,大约阅读时间需要 9 分钟。


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

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

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

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

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

转载于:https://www.cnblogs.com/Itst/p/11520615.html

你可能感兴趣的文章
ffmpeg的基本命令
查看>>
flex4+fms3.5+cs4开发实时音视频直播及点播详解
查看>>
打开VS2015提示“重新启动处于挂起状态。请在启动Visual Studio”之前重新启动
查看>>
form 中Enctype=multipart/form-data 的作用
查看>>
洛谷P1257 平面上的最接近点对
查看>>
BZOJ1657 [Usaco2006 Mar]Mooo 奶牛的歌声
查看>>
CSS/JavaScript hacks,browserhacks使用
查看>>
SQLServer内部数据库版本获取及匹配
查看>>
【转】<meta>标签用法
查看>>
ajax请求post和get的区别以及get post的选择
查看>>
前端优化之图片优化自动化
查看>>
指向指针的指针
查看>>
(转)VS自带工具:dumpbin的使用
查看>>
HTTP状态码
查看>>
回文数猜想(hd1282)
查看>>
svn 常用操作命令
查看>>
k8s的flannel网络插件配置
查看>>
构造函数
查看>>
并查集模板
查看>>
C++中的static关键字
查看>>