[洛谷P4168][Violet]蒲公英
题目大意:有$n(n\leqslant4\times10^4)$个数,$m(m\leqslant5\times10^4)$个询问,每次问区间$[l,r]$内的众数,若相同输出最小的,强制在线。
题解:先离散化,分块,设块大小为$S$,可以在$O(\dfrac n S n)$的复杂度内预处理出每两个块间的众数。考虑询问,发现询问中的众数要么是整块的那一个众数,要么就是非整块内出现过的数,可以用主席树查询区间数出现个数,比较一下即可,一次查询复杂度$O(2S\log_2 n)$。$S$开的比$\sqrt n$小一点就行了
卡点:主席树查询时左端点没有减一,数组开小
C++ Code:
#include <algorithm>
#include <cctype>
#include <cstdio>
#include <cstring>
namespace __IO {
namespace R {
int x, ch;
inline int read() {
while (isspace(ch = getchar()));
for (x = ch & 15; isdigit(ch = getchar());) x = x * 10 + (ch & 15);
return x;
}
}
}
using __IO::R::read; #define maxn 40010
const int BSZ = 100, Bnum = maxn / BSZ + 5; namespace SgT {
#define N (maxn * 20)
int lc[N], rc[N], V[N], idx; void insert(int &rt, const int l, const int r, const int num) {
lc[++idx] = lc[rt], rc[idx] = rc[rt], V[idx] = V[rt] + 1, rt = idx;
if (l == r) return ;
const int mid = l + r >> 1;
if (num <= mid) insert(lc[rt], l, mid, num);
else insert(rc[rt], mid + 1, r, num);
}
int query(const int L, const int R, const int l, const int r, const int num) {
if (l == r) return V[R] - V[L];
const int mid = l + r >> 1;
if (num <= mid) return query(lc[L], lc[R], l, mid, num);
else return query(rc[L], rc[R], mid + 1, r, num);
} #undef N
} int n, m, ans, anscnt;
int rt[maxn];
int v[maxn], s[maxn], bl[maxn];
int L[Bnum], R[Bnum], cnt[maxn];
int Max[Bnum][Bnum]; int main() {
n = read(), m = read();
for (int i = 1; i <= n; i++) {
v[i] = s[i] = read(); bl[i] = i / BSZ + 1;
}
int tot = (std::sort(v + 1, v + n + 1), std::unique(v + 1, v + n + 1) - v - 1);
for (int i = 1; i <= n; i++) {
SgT::insert(rt[i] = rt[i - 1], 1, tot, s[i] = std::lower_bound(v + 1, v + tot + 1, s[i]) - v);
} const int B = bl[n];
for (int i = 1; i <= B; i++) L[i] = (i - 1) * BSZ, R[i] = L[i] + BSZ - 1;
L[1] = 1, R[B] = n;
for (int l = 1, r, M; l <= B; l++) {
r = l, M = 0;
memset(cnt, 0, sizeof cnt);
for (int i = L[l]; i <= n; i++) {
cnt[s[i]]++;
if (cnt[s[i]] > cnt[M] || (cnt[s[i]] == cnt[M] && s[i] < M)) M = s[i];
if (i == R[r]) Max[l][r] = M, r++;
}
} while (m --> 0) {
int l = (read() + ans - 1) % n + 1, r = (read() + ans - 1) % n + 1;
if (l > r) std::swap(l, r);
const int lb = bl[l], rb = bl[r];
ans = anscnt = 0;
#define check(x) {\
const int tmp = SgT::query(rt[l - 1], rt[r], 1, tot, x); \
if (tmp > anscnt || (tmp == anscnt && x < ans)) ans = x, anscnt = tmp;\
}
if (lb == rb) {
for (register int i = l; i <= r; i++) check(s[i]);
} else {
for (register int i = l; i <= R[lb]; i++) check(s[i]);
if (lb + 1 <= rb - 1) check(Max[lb + 1][rb - 1]);
for (register int i = L[rb]; i <= r; i++) check(s[i]);
}
#undef check
printf("%d\n", ans = v[ans]);
}
return 0;
}
[洛谷P4168][Violet]蒲公英的更多相关文章
- 洛谷 P4168 [Violet]蒲公英 解题报告
P4168 [Violet]蒲公英 题目背景 亲爱的哥哥: 你在那个城市里面过得好吗? 我在家里面最近很开心呢.昨天晚上奶奶给我讲了那个叫「绝望」的大坏蛋的故事的说!它把人们的房子和田地搞坏,还有好多 ...
- 洛谷 P4168 [Violet] 蒲公英
历尽千辛万苦终于AC了这道题目... 我们考虑1个区间\([l,r]\), 被其完整包含的块的区间为\([L,R]\) 那么众数的来源? 1.\([l,L)\)或\((R,r]\)中出现的数字 2.\ ...
- 「分块系列」「洛谷P4168 [Violet]」蒲公英 解题报告
蒲公英 Description 我们把所有的蒲公英看成一个长度为\(n\)的序列(\(a_1,a_2,...a_n\)),其中\(a_i\)为一个正整数,表示第i棵蒲公英的种类的编号. 每次询问一个区 ...
- 洛谷P1445 [Violet] 樱花 (数学)
洛谷P1445 [Violet] 樱花 题目背景 我很愤怒 题目描述 求方程 1/X+1/Y=1/(N!) 的正整数解的组数,其中N≤10^6. 解的组数,应模1e9+7. 输入输出格式 输入格式: ...
- 洛谷P4168 蒲公英 分块处理区间众数模板
题面. 许久以前我还不怎么去机房的时候,一位大佬好像一直在做这道题,他称这道题目为"大分块". 其实这道题目的思想不只可以用于处理区间众数,还可以处理很多区间数值相关问题. 让我们 ...
- BZOJ2721或洛谷1445 [Violet]樱花
BZOJ原题链接 洛谷原题链接 其实推导很简单,只不过我太菜了想不到...又双叒叕去看题解 简单写下推导过程. 原方程:\[\dfrac{1}{x} + \dfrac{1}{y} = \dfrac{1 ...
- 【题解】洛谷P1445 [Violet]樱花 (推导+约数和)
洛谷P1445:https://www.luogu.org/problemnew/show/P1445 推导过程 1/x+1/y=1/n! 设y=n!+k(k∈N∗) 1/x+1/(n!+k)=1 ...
- 洛谷P4168 蒲公英 [Violet] 分块
题解:分块+离散化 解题报告: 一个分块典型题呢qwq还是挺妙的毕竟是道黑题 然,然后发现忘记放链接了先放链接QAQ 有两三种解法,都港下qwq 第一个是O(n5/3)的复杂度,谢总说不够优秀没有港, ...
- 【洛谷 P4168】[Violet]蒲公英(分块)
题目链接 题目大意:给定\(n\)个数和\(m\)个求区间众数的询问,强制在线 这题我\(debug\)了整整一个下午啊..-_- 从14:30~16:45终于\(debug\)出来了,\(debug ...
随机推荐
- Consul初体验
Preface Today I'm gonna implement a consul in my environment to discover service of MySQL da ...
- linux部署MantisBT(二)部署php
二.部署php 1.下载php安装包 http://php.net/downloads.php 2.安装php yum install libxml2 yum install libxml2-deve ...
- 使用jenkins构建一个自由风格的项目
一.创建一个freestyle风格的构建项目 二.输入任务名称和选择任务类型 三.配置项目 3.1:选择代码托管 3.2:到gitlab上去配置deploy key 3.3:到jenkins服务器去生 ...
- <cassert>
文件名: <cassert> (assert.h) 这是一个C语言的诊断库,assert.h文件中定义了一个可作为标准调试工具的宏函数: assert ; 下面介绍这个宏函数:asser ...
- Aizu - 2249
注意先保证距离最短,再来判断价格 邻接矩阵回朝内存 ,要用邻接表的 #include<bits/stdc++.h> using namespace std; #define inf 0x ...
- 【python 3.6】从网站抓图并存放到本地路径
#!/usr/bin/python # -*- coding: UTF-8 -*- _author_ = 'BH8ANK' import urllib.request import re import ...
- django启动创建用户失败
a django应用启动 b 访问127.0.0.1:8000,报错信息如下,原因为没有这个用户需要创建下用户 c 创建用户过程中报错原因是因为添加了app需要告诉django,这个 模型发生了改变, ...
- 贵州省未来二十年的投资机会的探讨1>
贵州的股市 1.000540.SZ 中天金融 2.000589.SZ 黔轮胎A 3.000733.SZ 振华科技 4.000851.SZ 高鸿股份 5.000920.SZ 南方汇通 6.002025. ...
- Python3 Tkinter-Entry
1.创建 from tkinter import * root=Tk() t1=Entry(root) t1.pack() root.mainloop() 2.绑定变量 from tkinter im ...
- 阿里校招内推C++岗位编程题第一题 空格最少的字符串
给定一个字符串S和有效单词的字典D,请确定可以插入到S中的最小空格数,使得最终的字符串完全由D中的有效单词组成.并输出解. 如果没有解则应该输出n/a 例如: 输入: S = “ilikealibab ...