洛谷 P4168 [Violet]蒲公英 解题报告
P4168 [Violet]蒲公英
题目背景
亲爱的哥哥:
你在那个城市里面过得好吗?
我在家里面最近很开心呢。昨天晚上奶奶给我讲了那个叫「绝望」的大坏蛋的故事的说!它把人们的房子和田地搞坏,还有好多小朋友也被它杀掉了。我觉得把那么可怕的怪物召唤出来的那个坏蛋也很坏呢。不过奶奶说他是很难受的时候才做出这样的事的……
最近村子里长出了一大片一大片的蒲公英。一刮风,这些蒲公英就能飘到好远的地方了呢。我觉得要是它们能飘到那个城市里面,让哥哥看看就好了呢!
哥哥你要快点回来哦!
爱你的妹妹 Violet
Azure 读完这封信之后微笑了一下。
“蒲公英吗……”
题目描述
在乡下的小路旁种着许多蒲公英,而我们的问题正是与这些蒲公英有关。
为了简化起见,我们把所有的蒲公英看成一个长度为n的序列 \((a_1,a_2..a_n)\),其中 \(a_i\)为一个正整数,表示第\(i\)棵蒲公英的种类编号。
而每次询问一个区间\([l,r]\),你需要回答区间里出现次数最多的是哪种蒲公英,如果有若干种蒲公英出现次数相同,则输出种类编号最小的那个。
注意,你的算法必须是在线的
输入输出格式
输入格式:
第一行两个整数 \(n\),\(m\) ,表示有\(n\)株蒲公英,\(m\)次询问。
接下来一行\(n\)个空格分隔的整数 \(a_i\),表示蒲公英的种类
再接下来\(m\)行每行两个整数 \(l_0,r_0\),我们令上次询问的结果为 \(x\)(如果这是第一次询问, 则 \(x=0\))。
令 \(l=(l_0+x-1)\bmod n + 1,r=(r_0+x-1) \bmod n + 1\),如果 \(l>r\),则交换 \(l,r\) 。
最终的询问区间为\([l,r]\)。
输出格式:
输出m 行。每行一个整数,表示每次询问的结果。
说明
对于 20% 的数据,保证 \(1\le n,m \le 3000\)。
对于 100% 的数据,保证 \(1\le n \le 40000,1\le m \le 50000,1\le a_i \le 10^9\)。
分块有时候还是蛮考思维哒?
这是一道经典的在线区间求众数的问题
设分成大小为\(S\)块
当询问区间\([l.r]\)时,我们把区间拆成\([l,L)\),\([L,R]\),\((R,r]\)三个区间
其中\(L,R\)为块的边界
答案只可能是 \([L,R]\)的众数以及在区间\([l,L)\)和区间\((R,r]\)出现的数字
我们对任意两个块所包含的大区间维护它的众数和每个数字的出现个数
后者需要一个长为\(n\)的数组存储
这样查询的时候,我们只需要枚举边角的每个数字出现个数就行啦,单次复杂度\(O(S)\)
预处理的话,我们需要枚举块的左右边界以及数字集,复杂度\(O(N*T^2)\)
则总复杂度为\(O(MS+N*T^2)\)
考虑块大小取多少时达到平衡
当\(MS=N*T^2\)时平衡
\(N ≈ M\)且\(S*T=N\)
所以\(S^2=N^3\)
\(S=N^{\frac{2}{3}}\)
总复杂度为\(O(N^{\frac{5}{3}})\)
如果用vector存+二分找似乎更快
Code:
#include <cstdio>
#include <cmath>
#include <map>
using namespace std;
const int N=4e4+10;
const int T=40;
int n,m,r,a[N],b[N],seg[T][T][N],num[T][T],L[T],R[T],pos[N],buct[N];
std::map <int,int > ma;
int query(int l,int r)
{
int mx=0,ans,ll=pos[l]+1,rr=pos[r]-1;
if(ll>rr)
{
for(int i=l;i<=r;i++)
{
++buct[a[i]];
if(buct[a[i]]>mx||(buct[a[i]]==mx&&a[i]<ans))
{
mx=buct[a[i]];
ans=a[i];
}
}
for(int i=l;i<=r;i++)
buct[a[i]]=0;
}
else
{
mx=seg[ll][rr][0],ans=num[ll][rr];
for(int i=l;i<L[ll];i++)
{
++buct[a[i]];
if(buct[a[i]]+seg[ll][rr][a[i]]>mx||(buct[a[i]]+seg[ll][rr][a[i]]==mx&&a[i]<ans))
{
mx=buct[a[i]]+seg[ll][rr][a[i]];
ans=a[i];
}
}
for(int i=R[rr]+1;i<=r;i++)
{
++buct[a[i]];
if(buct[a[i]]+seg[ll][rr][a[i]]>mx||(buct[a[i]]+seg[ll][rr][a[i]]==mx&&a[i]<ans))
{
mx=buct[a[i]]+seg[ll][rr][a[i]];
ans=a[i];
}
}
for(int i=l;i<L[ll];i++)
buct[a[i]]=0;
for(int i=R[rr]+1;i<=r;i++)
buct[a[i]]=0;
}
return ans;
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
scanf("%d",a+i),ma[a[i]]=1;
for(map <int,int>::iterator it=ma.begin();it!=ma.end();it++)
it->second=++r;
for(int i=1;i<=n;i++)
b[ma[a[i]]]=a[i],a[i]=ma[a[i]];
int s=pow(n,0.6666667),t;
for(t=1;t*s<=n;t++)
L[t]=(t-1)*s+1,R[t]=t*s;
if(R[t-1]<n) L[t]=R[t-1]+1,R[t]=n;
else --t;
for(int i=1;i<=t;i++)
for(int j=L[i];j<=R[i];j++)
pos[j]=i;
for(int i=1;i<=t;i++)
for(int j=i;j<=t;j++)
{
for(int k=L[i];k<=R[j];k++)
{
++seg[i][j][a[k]];
if(seg[i][j][a[k]]>seg[i][j][0]||(seg[i][j][a[k]]==seg[i][j][0]&&a[k]<num[i][j]))
{
num[i][j]=a[k];
seg[i][j][0]=seg[i][j][a[k]];
}
}
}
int lastans=0;
for(int l,r,i=1;i<=m;i++)
{
scanf("%d%d",&l,&r);
l=(l+lastans-1)%n+1,r=(r+lastans-1)%n+1;
if(l>r) swap(l,r);
printf("%d\n",lastans=b[query(l,r)]);
}
return 0;
}
2018.8.27
洛谷 P4168 [Violet]蒲公英 解题报告的更多相关文章
- 洛谷 P4168 [Violet] 蒲公英
历尽千辛万苦终于AC了这道题目... 我们考虑1个区间\([l,r]\), 被其完整包含的块的区间为\([L,R]\) 那么众数的来源? 1.\([l,L)\)或\((R,r]\)中出现的数字 2.\ ...
- [洛谷P4168][Violet]蒲公英
题目大意:有$n(n\leqslant4\times10^4)$个数,$m(m\leqslant5\times10^4)$个询问,每次问区间$[l,r]$内的众数,若相同输出最小的,强制在线. 题解: ...
- 洛谷 P1783 海滩防御 解题报告
P1783 海滩防御 题目描述 WLP同学最近迷上了一款网络联机对战游戏(终于知道为毛JOHNKRAM每天刷洛谷效率那么低了),但是他却为了这个游戏很苦恼,因为他在海边的造船厂和仓库总是被敌方派人偷袭 ...
- 洛谷 P4597 序列sequence 解题报告
P4597 序列sequence 题目背景 原题\(\tt{cf13c}\)数据加强版 题目描述 给定一个序列,每次操作可以把某个数\(+1\)或\(-1\).要求把序列变成非降数列.而且要求修改后的 ...
- 洛谷1087 FBI树 解题报告
洛谷1087 FBI树 本题地址:http://www.luogu.org/problem/show?pid=1087 题目描述 我们可以把由“0”和“1”组成的字符串分为三类:全“0”串称为B串,全 ...
- 「分块系列」「洛谷P4168 [Violet]」蒲公英 解题报告
蒲公英 Description 我们把所有的蒲公英看成一个长度为\(n\)的序列(\(a_1,a_2,...a_n\)),其中\(a_i\)为一个正整数,表示第i棵蒲公英的种类的编号. 每次询问一个区 ...
- 洛谷 P3349 [ZJOI2016]小星星 解题报告
P3349 [ZJOI2016]小星星 题目描述 小\(Y\)是一个心灵手巧的女孩子,她喜欢手工制作一些小饰品.她有\(n\)颗小星星,用\(m\)条彩色的细线串了起来,每条细线连着两颗小星星. 有一 ...
- 洛谷 P3177 树上染色 解题报告
P3177 [HAOI2015]树上染色 题目描述 有一棵点数为\(N\)的树,树边有边权.给你一个在\(0\) ~ \(N\)之内的正整数\(K\),你要在这棵树中选择\(K\)个点,将其染成黑色, ...
- 洛谷 P4705 玩游戏 解题报告
P4705 玩游戏 题意:给长为\(n\)的\(\{a_i\}\)和长为\(m\)的\(\{b_i\}\),设 \[ f(x)=\sum_{k\ge 0}\sum_{i=1}^n\sum_{j=1}^ ...
随机推荐
- 【tp5.1】composer安装PHPExcel以及导入\导出Excel
一.安装PHPExcel 1.下载:PHPExcel https://github.com/PHPOffice/PHPExcel 2.解压后:Classes文件夹改名为PHPExcel 3.把文件夹 ...
- Typora -- 书写即美学
#Typora -- 书写即美学 ##基本快捷键--需要在所见即所想界面进行输入 加粗 Ctrl + B 加粗 斜体 Ctrl + I 斜体 下划线 Ctrl + U 下划线 删除线 Ctrl + S ...
- linux下磁盘分区、格式化、挂载
(1)fdisk /dev/sdb进行分区 (2)选择n表示添加一个分区,选择d表示删除一个分区.可通过m获取帮助信息 (3)选择p表示主分区,然后输入分区大小 (4)分区完成后,可通过fdisk - ...
- HDU 5530:Pipes Selection
题意: 给定长度为\(L\),元素总和为\(S\)的非负整数序列\(A\),对于每一个\(1 \leq i \leq S\),求出:所有满足\(\sum_{j=l}^rA_j=i\)的二元组\((l, ...
- 1、Java多线程基础:进程和线程之由来
Java多线程基础:进程和线程之由来 在前面,已经介绍了Java的基础知识,现在我们来讨论一点稍微难一点的问题:Java并发编程.当然,Java并发编程涉及到很多方面的内容,不是一朝一夕就能够融会贯通 ...
- Hystrix入门指南
Introduction 1.Where does the name come from? hystrix对应的中文名字是“豪猪”,豪猪周身长满了刺,能保护自己不受天敌的伤害,代表了一种防御机制,这与 ...
- 完整的vue+vuex+api-router+database请求流程
- 虚拟接VMnet1 和VMnet8的区别
vmnet1是host-only,也就是说,选择用vmnet1的话就相当于VMware给你提供了一个虚拟交换机,仅将虚拟机和真实系统连上了,虚拟机可以与真实系统相互共享文件,但是虚拟机无法访问外部互联 ...
- XML与Object的范型转换
前段时间做object转换xml想了很多,所有打算整理下 做成以下的通用方法. public static bool ObjectToXml<T>(string filePath, T t ...
- 【题解搬运】PAT_A1016 Phone Bills
从我原来的博客上搬运.原先blog作废. 题目 A long-distance telephone company charges its customers by the following rul ...