博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
数列分块入门 9 / 蒲公英
阅读量:6079 次
发布时间:2019-06-20

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

区间众数

这题是强制在线的,所以只有分块了,不然的话还可以莫队。。

预处理出每块之间的众数以及每个数出现的位置,查询的时候块内直接查询,不足整块的直接二分l和r,相减就是一个数出现的次数了,然后更新答案即可。

#include 
#define INF 0x3f3f3f3f#define full(a, b) memset(a, b, sizeof a)using namespace std;typedef long long ll;inline int lowbit(int x){ return x & (-x); }inline int read(){ int X = 0, w = 0; char ch = 0; while(!isdigit(ch)) { w |= ch == '-'; ch = getchar(); } while(isdigit(ch)) X = (X << 3) + (X << 1) + (ch ^ 48), ch = getchar(); return w ? -X : X;}inline int gcd(int a, int b){ return b ? gcd(b, a % b) : a; }inline int lcm(int a, int b){ return a / gcd(a, b) * b; }template
inline T max(T x, T y, T z){ return max(max(x, y), z); }template
inline T min(T x, T y, T z){ return min(min(x, y), z); }template
inline A fpow(A x, B p, C lyd){ A ans = 1; for(; p; p >>= 1, x = 1LL * x * x % lyd)if(p & 1)ans = 1LL * x * ans % lyd; return ans;}const int N = 100005;int n, t, id, lt[N], rt[N], a[N], f[3005][3005], freq[N], val[N], pos[N], c, _;map
m;vector
v[N];void pretreatment(int k){ full(freq, 0); int ans = 0, mx = 0; for(int i = lt[k]; i <= n; i ++){ freq[a[i]] ++; if(freq[a[i]] > mx || (freq[a[i]] == mx && val[a[i]] < val[ans])) mx = freq[a[i]], ans = a[i]; f[k][pos[i]] = ans; }}int calc(int x, int l, int r){ int a = (int)(upper_bound(v[x].begin(), v[x].end(), r) - v[x].begin()); int b = (int)(lower_bound(v[x].begin(), v[x].end(), l) - v[x].begin()); return a - b;}int query(int l, int r){ int p = pos[l], q = pos[r], ans = 0, mx = 0; if(p == q){ for(int i = l; i <= r; i ++){ int k = calc(a[i], l, r); if(k > mx || (k == mx && val[a[i]] < val[ans])) mx = k, ans = a[i]; } } else{ ans = f[p + 1][q - 1], mx = calc(ans, l, r); for(int i = l; i <= rt[p]; i ++){ int k = calc(a[i], l, r); if(k > mx || (k == mx && val[a[i]] < val[ans])) mx = k, ans = a[i]; } for(int i = lt[q]; i <= r; i ++){ int k = calc(a[i], l, r); if(k > mx || (k == mx && val[a[i]] < val[ans])) mx = k, ans = a[i]; } } return ans;}int main(){ //freopen("data.txt", "r", stdin); n = read(), _ = read(); for(int i = 1; i <= n; i ++){ a[i] = read(); if(!m[a[i]]) m[a[i]] = ++id, val[id] = a[i]; a[i] = m[a[i]]; v[a[i]].push_back(i); } t = (int)sqrt(_ * ((int)(log2(0.1*n)) + 1)); c = n / t; for(int i = 1; i <= t; i ++){ lt[i] = (i - 1) * c + 1; rt[i] = i * c; } if(rt[t] < n) t ++, lt[t] = rt[t - 1] + 1, rt[t] = n; for(int i = 1; i <= t; i ++){ for(int j = lt[i]; j <= rt[i]; j ++){ pos[j] = i; } } for(int i = 1; i <= t; i ++) pretreatment(i); int x = 0; for(int i = 1; i <= _; i ++){ int l = read(), r = read(); l = (l + x - 1) % n + 1, r = (r + x - 1) % n + 1; if(l > r) swap(l, r); printf("%d\n", val[query(l, r)]); x = val[query(l, r)]; } return 0;}

转载于:https://www.cnblogs.com/onionQAQ/p/10888315.html

你可能感兴趣的文章
redis分布式锁redisson
查看>>
什么样的企业可以称之为初创企业?
查看>>
Python爬虫之BeautifulSoup
查看>>
《HTML 5与CSS 3权威指南(第3版·下册)》——第20章 使用选择器在页面中插入内容...
查看>>
如何判断自己适不适合做程序员?这几个特点了解一下
查看>>
newinstance()和new有什么区别
查看>>
android下载封装类
查看>>
[node] 用 node-webkit 开发桌面应用
查看>>
Nginx访问控制和虚拟主机
查看>>
report widget not working for external users
查看>>
windows phone 摄像头得到图片是旋转90°
查看>>
Linux--sed使用
查看>>
没有显示器的情况下安装和使用树莓派
查看>>
Q85 最大矩形
查看>>
【android】使用handler更新UI
查看>>
mochiweb 源码阅读(十五)
查看>>
前端面试中的常见的算法问题
查看>>
计算机语言的基本理论
查看>>
nodejs流之行读取器例子
查看>>
批量文件重命名工具
查看>>