BZOJ.2724.[Violet 6]蒲公英(静态分块)
区间众数 强制在线
考虑什么样的数会成为众数
如果一个区间S1的众数为x,那么S1与新区间S2的并的众数只会是x或S2中的数
所以我们可以分块先预处理f[i][j]表示第i到第j块的众数
对于零散部分,我们还需要知道它们在区间中的出现次数。这部分至多有2sqrt(n)个
由于没有修改,离散化后对于每个数x开一个vector,把x出现位置push_back进去,查x时二分即可。
像普通分块一样更新即可。之前写了一堆特判 醉了
这个块大小怎么算。。没算出来。
//2520kb 19652ms 这么慢。。
#include <cmath>
#include <cctype>
#include <cstdio>
#include <vector>
#include <cstring>
#include <algorithm>
#define gc() getchar()
const int N=4e4+5,S=205;
int n,m,size,f[S][S],t[S][S],A[N],ref[N],bel[N],tm[N];
std::vector<int> v[N];
inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now;
}
int Find(int x,int r)
{
int l=1,mid;
while(l<r)
if(ref[mid=l+r>>1]>=x) r=mid;
else l=mid+1;
return l;
}
void Discrete()
{
for(int i=1; i<=n; ++i) ref[i]=A[i]=read();
std::sort(ref+1,ref+1+n);
int cnt=1;
for(int i=2; i<=n; ++i)
if(ref[i]!=ref[i-1]) ref[++cnt]=ref[i];
for(int i=1; i<=n; ++i) A[i]=Find(A[i],cnt);
}
void Init()
{
Discrete();
for(int i=1; i<=n; ++i) bel[i]=(i-1)/size+1;
for(int i=1; i<=n; ++i) v[A[i]].push_back(i);
for(int i=1; i<=bel[n]; ++i)
{
memset(tm,0,sizeof tm);
int id=0,mxt=0;
for(int j=i; j<=bel[n]; ++j)
{
int r=std::min(size*j,n);
for(int k=(j-1)*size+1; k<=r; ++k)
if(++tm[A[k]]>mxt) id=A[k],mxt=tm[id];
else if(tm[A[k]]==mxt && A[k]<id) id=A[k];
f[i][j]=id, t[i][j]=mxt;
}
}
}
int Lower(int id,int x)
{
int l=0,r=v[id].size()-1,mid;
while(l<r)
if(v[id][mid=l+r>>1]>=x) r=mid;
else l=mid+1;
return l;
}
int Upper(int id,int x)
{
int l=0,r=v[id].size(),mid;
while(l+1<r)
{
if(v[id][mid=l+r>>1]<=x) l=mid;
else r=mid;
}
return l;
}
inline int Query(int l,int r,int id){
if(!id) return 0;
return Upper(id,r)-Lower(id,l)+1;
}
void Update(int l,int r,int L,int R,int &id,int &mxt)
{
for(int tmp,i=l; i<=r; ++i)
if((tmp=Query(L,R,A[i]))>mxt||(tmp==mxt&&id>A[i]))
id=A[i], mxt=tmp;
}
int Solve(int l,int r)
{
int id=f[bel[l]+1][bel[r]-1],mxt=t[bel[l]+1][bel[r]-1];//mxt=Query(l,r,id);//这样id可能不存在 为0,Query中必须要特判
Update(l,std::min(r,bel[l]*size),l,r,id,mxt);
if(bel[l]!=bel[r]) Update((bel[r]-1)*size+1,r,l,r,id,mxt);
return id;
}
int main()
{
n=read(),m=read(); size=sqrt(n);//200
Init();
int res=0,l,r;
while(m--)
{
l=(read()+res-1)%n+1, r=(read()+res-1)%n+1;
if(l>r) std::swap(l,r);
printf("%d\n",res=ref[Solve(l,r)]);
}
return 0;
}
BZOJ.2724.[Violet 6]蒲公英(静态分块)的更多相关文章
- [BZOJ 2724] [Violet 6] 蒲公英 【分块】
题目链接:BZOJ - 2724 题目分析 这道题和 BZOJ-2821 作诗 那道题几乎是一样的,就是直接分块,每块大小 sqrt(n) ,然后将数字按照数值为第一关键字,位置为第二关键字排序,方便 ...
- BZOJ 2724: [Violet 6]蒲公英( 分块 )
虽然AC了但是时间惨不忍睹...不科学....怎么会那么慢呢... 无修改的区间众数..分块, 预处理出Mode[i][j]表示第i块到第j块的众数, sum[i][j]表示前i块j出现次数(前缀和, ...
- BZOJ 2724: [Violet 6]蒲公英
2724: [Violet 6]蒲公英 Time Limit: 40 Sec Memory Limit: 512 MBSubmit: 1633 Solved: 563[Submit][Status ...
- BZOJ 2724: [Violet 6]蒲公英 [分块 区间众数]
传送门 题面太美不忍不放 分块分块 这种题的一个特点是只有查询,通常需要预处理:加入修改的话需要暴力重构预处理 预处理$f[i][j]$为第i块到第j块的众数,显然$f[i][j]=max{f[i][ ...
- BZOJ 2724 [Violet 6]蒲公英(分块)
题意 在线区间众数 思路 预处理出 f[i][j] 即从第 i 块到第 j 块的答案.对于每个询问,中间的整块直接用预处理出的,两端的 sqrtn 级别的数暴力做,用二分查找它们出现的次数.每次询问的 ...
- 【刷题】BZOJ 2724 [Violet 6]蒲公英
Description Input 修正一下 l = (l_0 + x - 1) mod n + 1, r = (r_0 + x - 1) mod n + 1 Output Sample Input ...
- 【BZOJ 2724】 2724: [Violet 6]蒲公英 (区间众数不带修改版本)
2724: [Violet 6]蒲公英 Time Limit: 40 Sec Memory Limit: 512 MBSubmit: 1908 Solved: 678 Description In ...
- 【BZOJ】2724: [Violet 6]蒲公英
2724: [Violet 6]蒲公英 Time Limit: 40 Sec Memory Limit: 512 MBSubmit: 2900 Solved: 1031[Submit][Statu ...
- BZOJ_2724_[Violet 6]蒲公英_分块
BZOJ_2724_[Violet 6]蒲公英_分块 Description Input 修正一下 l = (l_0 + x - 1) mod n + 1, r = (r_0 + x - 1) mod ...
随机推荐
- 简单透彻理解JSONP原理及使用
首先提一下JSON这个概念,JSON是一种轻量级的数据传输格式,被广泛应用于当前Web应用中.JSON格式数据的编码和解析基本在所有主流语言中都被实现,所以现在大部分前后端分离的架构都以JSON格式进 ...
- 图文并茂的Python教程-numpy.pad
图文并茂的Python教程-numpy.pad np.pad()常用与深度学习中的数据预处理,可以将numpy数组按指定的方法填充成指定的形状. 声明: 需要读者了解一点numpy数组的知识np.pa ...
- day1 diff命令递归比较目录下的文件
一.diff实战 (1)递归比较文件夹下所有的文件及目录的不同 diff --brief -Nr dir1/ dir2/ Reference ...
- [译]bootstrap-select (selectpicker)方法
方法 .selectpicker('val') 您可以通过调用val元素上的方法来设置所选值. $('.selectpicker').selectpicker('val', 'Mustard'); $ ...
- 关卡得分(if 嵌套for)与(for嵌套if)
- Android 截屏与 WebView 长图分享经验总结
最近在做新业务需求的同时,我们在 Android 上遇到了一些之前没有碰到过的问题,截屏分享. WebView 生成长图以及长图在各个分享渠道分享时图片模糊甚至分享失败等问题,在这过程中踩了很多坑,到 ...
- how to detect circles and rectangle?
opencv中对圆检测的函数为:HoughCircles(src_gray,circles,CV_HOUGHT_GRADIENT,1,src_gray.cols/8,200,100,0,0) circ ...
- LeetCode(53):最大子序和
Easy! 题目描述: 给定一个整数数组 nums ,找到一个具有最大和的连续子数组(子数组最少包含一个元素),返回其最大和. 示例: 输入: [-2,1,-3,4,-1,2,1,-5,4], 输出: ...
- 使用caffe模型测试图片(python接口)
1.加载相关模块 1.1 加载numpy import numpy as np 1.2 加载caffe 有两种方法. 方法一(静态导入): 找到当前环境使用的python的site-packages目 ...
- OGNL(对象图导航语言)学习
一.关于OGNL(Object-Graph Navigation Language),一种可以方便地操作对象属性的开源表达式语言. 特点: 1)支持对象方法调用,形式如:objName.method ...