LGP3709题解
题目大意
简化后为区间众数出现次数,简化前为【数据删除】
吐槽
为什么题解只有一篇分块,剩下的全是莫队?
这题不是蒲公英?这和算导例题有何区别???
为什么现在的人都喜欢去看题解而不注重思维???
莫队之前也胡过区间众数莫队,由于太菜胡出来了一个回滚莫队。(毕竟暴力思路这题很难删)
题解
因为某些原因代码是从 P5048 搬过来然后把强制在线去掉后加了个负号(
首先我们对序列分块,然后求出左右端点同时也是某个块的左右端点的区间的众数出现次数。(语文差轻喷)
然后对于一个询问,把询问区间分为整块和边角块后,我们很容易能够知道答案一定在 O(整块答案) 级别。所以我们枚举边角块的 \(O(\sqrt n)\) 个数,看它在区间中出现了多少次就行了。
但是数出现次数的方法有点儿妙。
我们开 \(n\) 个 vector 记录每一种数出现的下标。
假如一个边角块的数可能成为答案(设其为 \(k\)),那么这个数后面的第整块答案个 \(k\) 一定也在区间内。
而很容易知道更新答案的复杂度是 \(O(\sqrt n)\) 的。
code:
#include<algorithm>
#include<cstdio>
#include<vector>
#include<cmath>
const int N=805,M=5e5+5;
int n,m,p,ans,siz,len,a[M],t[M],pre[M],cnt[M],lsh[M],ANS[N][N];
int L[N],R[N],pos[M];
bool vis[M];
std::vector<int>id[M];
inline int min(const int&a,const int&b){
return a>b?b:a;
}
inline int Query(int l,int r){
int i,j,q=pos[l],p=pos[r];
if(q+1>=p||q==p){
for(i=l;i<=r;++i){
if(++cnt[a[i]]>ans)ans=cnt[a[i]];
}
for(i=l;i<=r;++i)--cnt[a[i]];
return ans;
}
ans=ANS[q+1][p-1];
for(i=l;i<=R[q];++i){
for(j=pre[i]+ans;j<id[a[i]].size()&&id[a[i]][j]<=r;++j)++ans;
}
for(i=L[p];i<=r;++i){
for(j=pre[i]-ans;j>=0&&id[a[i]][j]>=l;--j)++ans;
}
return ans;
}
signed main(){
register int i,j,k,x=0,ql,qr;
scanf("%d%d",&n,&m);p=ceil(n/sqrt(m));
for(i=1;i<=n;++i){
scanf("%d",a+i);lsh[i]=a[i];
pos[i]=(i-1)/p+1;
}
siz=pos[n];
for(i=1;i<=siz;++i)L[i]=R[i-1]+1,R[i]=i*p;
R[siz]=n;
std::sort(lsh+1,lsh+n+1);len=std::unique(lsh+1,lsh+n+1)-lsh-1;
for(i=1;i<=n;++i){
a[i]=std::lower_bound(lsh+1,lsh+len+1,a[i])-lsh;
id[a[i]].push_back(i);pre[i]=id[a[i]].size()-1;
}
for(i=1;i<=siz;++i){
for(j=i;j<=siz;++j){
for(k=L[j];k<=R[j];++k){
if(++cnt[a[k]]>ans)ans=cnt[a[k]];
}
ANS[i][j]=ans;
}
for(j=L[i];j<=n;++j)--cnt[a[j]];ans=0;
}
for(i=1;i<=m;++i){
scanf("%d%d",&ql,&qr);ans=len=0;
printf("%d\n",-Query(ql,qr));
}
}
LGP3709题解的更多相关文章
- 2016 华南师大ACM校赛 SCNUCPC 非官方题解
我要举报本次校赛出题人的消极出题!!! 官方题解请戳:http://3.scnuacm2015.sinaapp.com/?p=89(其实就是一堆代码没有题解) A. 树链剖分数据结构板题 题目大意:我 ...
- noip2016十连测题解
以下代码为了阅读方便,省去以下头文件: #include <iostream> #include <stdio.h> #include <math.h> #incl ...
- BZOJ-2561-最小生成树 题解(最小割)
2561: 最小生成树(题解) Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1628 Solved: 786 传送门:http://www.lyd ...
- Codeforces Round #353 (Div. 2) ABCDE 题解 python
Problems # Name A Infinite Sequence standard input/output 1 s, 256 MB x3509 B Restoring P ...
- 哈尔滨理工大学ACM全国邀请赛(网络同步赛)题解
题目链接 提交连接:http://acm-software.hrbust.edu.cn/problemset.php?page=5 1470-1482 只做出来四道比较水的题目,还需要加强中等题的训练 ...
- 2016ACM青岛区域赛题解
A.Relic Discovery_hdu5982 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Jav ...
- poj1399 hoj1037 Direct Visibility 题解 (宽搜)
http://poj.org/problem?id=1399 http://acm.hit.edu.cn/hoj/problem/view?id=1037 题意: 在一个最多200*200的minec ...
- 网络流n题 题解
学会了网络流,就经常闲的没事儿刷网络流--于是乎来一发题解. 1. COGS2093 花园的守护之神 题意:给定一个带权无向图,问至少删除多少条边才能使得s-t最短路的长度变长. 用Dijkstra或 ...
- CF100965C题解..
求方程 \[ \begin{array}\\ \sum_{i=1}^n x_i & \equiv & a_1 \pmod{p} \\ \sum_{i=1}^n x_i^2 & ...
随机推荐
- Myeclipse设置关键词自动提示功能以及取消空格和"="的自动补全
一.设置Myeclipse从a~z的自动提示功能 1. "window"→"Preferences"2. 选择"Java",展开," ...
- opencv结构IplImage
转载请注明来源:https://www.cnblogs.com/hookjc/ typedef struct _IplImage{int nSize; /* Ip ...
- ANT中的copy和move标签
Copy标签 该标签用于文件或文件集的拷贝,其属性如下: file 表示源文件. tofile 表示目标文件. todir 表示目标目录. overwrite 表示指定是否覆盖目标文件,默认值是不覆盖 ...
- cell重用
少数几个cell可不重用 NSString *CellIdentifier = [NSString stringWithFormat:@"MyCellID_%d",indexPat ...
- 带你十天轻松搞定 Go 微服务系列(九、链路追踪)
序言 我们通过一个系列文章跟大家详细展示一个 go-zero 微服务示例,整个系列分十篇文章,目录结构如下: 环境搭建 服务拆分 用户服务 产品服务 订单服务 支付服务 RPC 服务 Auth 验证 ...
- java+selenium脚本编写规范
2. 源文件规范 2.1 文件名 源文件以最顶层的类名来命名,大小写敏感,文件扩展名为.java 2.2 文件编码 UTF-8 源文件要求编码格式为UTF-8 2.3 源文件结 ...
- [技术干货-算子使用] mindspore.scipy 入门使用指导
1. MindSpore框架的SciPy模块 SciPy 是基于NumPy实现的科学计算库,主要用于数学.物理学.生物学等科学以及工程学领域.诸如高阶迭代,线性代数求解等都会需要用到SicPy.Sci ...
- 我们一起来学Shell - shell的循环控制
文章目录 Shell 循环之 for 语句 Shell 循环之 while 语句 Shell 循环之 until 语句 Shell 循环控制 break指令 continue 指令 exit 指令 s ...
- 通过Dapr实现一个简单的基于.net的微服务电商系统(十九)——分布式事务之Saga模式
在之前的系列文章中聊过分布式事务的一种实现方案,即通过在集群中暴露actor服务来实现分布式事务的本地原子化.但是actor服务本身有其特殊性,场景上并不通用.所以今天来讲讲分布式事务实现方案之sag ...
- swagger 2.0
1.引入jar包 <dependency> <groupId>io.springfox</groupId> <artifactId>springfox- ...