Description

转眼就要到Karin的生日了!Yuuna她们想为她准备生日礼物!现在有许多礼物被排列成了一个一维序列,每个礼物都有一个价值。Yuuna对这个序列十分感兴趣。因此,你需要多次回答:在某个区间内出现次数第k1少的价值是多少,可能多个不同的价值出现次数均为第k1少,输出其中第k2小的,保证输入合法。注意内存限制。
例如:对于一个区间而言(当然不一定是有序的):
1,2,3,4,5,5,6,6,7,7,7,7,8,8,8,8,9,9,9,9,10,10,10,10,11,11,11,11,11,11,12,12,12,12,12,12,12,12,12,12
权值 出现次数
1 1 出现次数第1少 第1小
2 1 第2小
3 1 第3小
4 1 第4小
5 2 出现次数第2少 第1小
6 2 第2小
7 4 出现次数第3少 第1小
8 4 第2小
9 4 第3小
10 4 第4小
11 6 出现次数第4少 第1小
12 10 出现次数第5少 第1小
若k1=3,k2=2,代表询问这个区间里出现次数第3少的权值中第2小的,则应该输出8。
若k1=5,k2=1,代表询问这个区间里出现次数第5少的权值中第1小的,则应该输出12。
若k1=1,k2=3,代表询问这个区间里出现次数第1少的权值中第3小的,则应该输出3。

Input

第一行包括一个整数n,代表序列的长度。
第二行包括n个整数a1...an,代表该序列。
第三行包括一个整数m,代表询问的次数。
接下来m行,每行包括4个整数l,r,k1,k2,询问al...ar中出现次数第k1少的权值中第k2小的。

Output

对于每个询问,仅输出一行,包括一个整数,代表你的回答。

很显然可以用莫队算法,关键在于维护状态转移。权值分块可以做到O(1)插入/删除,O(√(n))查询第k大,可以用于维护每种出现次数有几个权值以及相同出现次数的不同权值。由均摊分析可知,每种出现次数的权值数之和为n,于是可以预处理离散化并给每个出现次数的权值分块分配对应的空间。可以做到O(n)预处理,O(m√(n))查询,O(n)空间,但常数较大。

#include<cstdio>
#include<algorithm>
#include<cmath>
int _(){
int x=,c=getchar();
while(c<)c=getchar();
while(c>)x=x*+c-,c=getchar();
return x;
}
int n,m,v[],B,id[],ts[],ans[],t1[],t2[],mt[],mc[];
int*l0[],*r0[],*l1[],*l2[],*l3[],_mem[],*mp=_mem;
struct Q{
int l,r,k1,k2,ID;
void init(int i){
l=_();r=_();k1=_();k2=_();ID=i;
}
bool operator<(const Q&w)const{
return id[l]!=id[w.l]?id[l]<id[w.l]:r!=w.r?(r<w.r)!=(id[l]&):ID<w.ID;
}
}q[];
void insv(int x,int y){
if(!t1[x]++)++t2[x/B];
y=l3[y][x];
++l1[x][y];
++l2[x][y/B];
}
void delv(int x,int y){
if(!--t1[x])--t2[x/B];
y=l3[y][x];
--l1[x][y];
--l2[x][y/B];
}
void ins(int x){
if(ts[x]>)delv(ts[x],x);
++ts[x];
if(ts[x]>)insv(ts[x],x);
}
void del(int x){
if(ts[x]>)delv(ts[x],x);
--ts[x];
if(ts[x]>)insv(ts[x],x);
}
int rnk(int x,int k){
for(int i=,s=;;++i){
if(s+l2[x][i]>=k){
for(int j=i*B;;++j)if((s+=l1[x][j])==k)return l0[x][j];
}
s+=l2[x][i];
}
}
int main(){
n=_();
B=sqrt(n+)+;
for(int i=;i<=n;++i)++mc[++mt[v[i]=_()]],id[i]=(i-)/B;
for(int i=;i<=n;++i){
l0[i]=r0[i]=mp;mp+=mc[i];
l1[i]=mp;mp+=mc[i];
l2[i]=mp;mp+=mc[i];
l3[i]=mp-;mp+=mt[i];
}
for(int i=;i<=n;++i)for(int j=;j<=mt[i];++j)l3[i][j]=r0[j]-l0[j],*r0[j]++=i;
m=_();
for(int i=;i<m;++i)q[i].init(i);
std::sort(q,q+m);
int L=,R=;
for(int i=;i<m;++i){
int l=q[i].l,r=q[i].r;
while(L>l)ins(v[--L]);
while(L<l)del(v[L++]);
while(R<r)ins(v[++R]);
while(R>r)del(v[R--]);
for(int j=,s=;;++j){
if(s+t2[j]>=q[i].k1){
int k=B*j;
while(s<q[i].k1)if(t1[k++])++s;--k;
ans[q[i].ID]=rnk(k,q[i].k2);
break;
}
s+=t2[j];
}
}
for(int i=;i<m;++i)printf("%d\n",ans[i]);
return ;
}

bzoj 3920: Yuuna的礼物的更多相关文章

  1. BZOJ.3920.Yuuna的礼物(莫队 分块套分块 分段离散化)

    题目链接 详细题解:https://www.cnblogs.com/autsky-jadek/p/4376091.html 代码参考自:https://www.cnblogs.com/Sakits/p ...

  2. bzoj3920: Yuuna的礼物(莫队+分块套分块)

    思路挺简单的,但是总感觉好难写...码力还是差劲,最后写出来也挺丑的 这题显然是个莫队题,考虑怎么转移和询问... 根据莫队修改多查询少的特点,一般用修改快查询慢的分块来维护.查第$k_1$小的出现次 ...

  3. 【莫队算法】【权值分块】bzoj3920 Yuuna的礼物

    [算法一] 暴力. 可以通过第0.1号测试点. 预计得分:20分. [算法二] 经典问题:区间众数,数据范围也不是很大,因此我们可以: ①分块,离散化,预处理出: <1>前i块中x出现的次 ...

  4. [BZOJ3920]Yuuna的礼物

    题目大意: 给你一个长度为$n(n\le40000)$的数列$\{a_i\}(1\le a_i\le n)$,给出$m(m\le40000)$次询问,每次给出$l,r,k_1,k_2$询问区间$[l, ...

  5. bzoj AC倒序

    Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...

  6. bzoj 3055礼物运送 floyed + 状压DP

    bzoj 3055: 礼物运送 floyed first 设f[i][S]表示取到了S集合中的所有点(不一定是经过的所有点),最后停在了i的最优值. 初始就f[i][{i}] = dis[1][i] ...

  7. bzoj 4827 礼物

    bzoj 4827 礼物 可以看做将其中一个数列(假定为 \(a\) )都加上 \(c\) , \(c\) 可以为负数.易知这里 \(-m\leq c\leq m\). 记要求的答案为 \(ans\) ...

  8. BZOJ 3438 小M的礼物

    BZOJ 3438 小M的礼物 Description 小M在MC里开辟了两块巨大的耕地A和B(你可以认为容量是无穷),现在,小P有n中作物的种子,每种作物的种子 有1个(就是可以种一棵作物)(用1. ...

  9. bzoj 4827: [Hnoi2017]礼物 [fft]

    4827: [Hnoi2017]礼物 题意:略 以前做的了 化一化式子就是一个卷积和一些常数项 我记着确定调整值还要求一下导... #include <iostream> #include ...

随机推荐

  1. Spring3整合Hibernate4-我们到底能走多远系列(30)

    我们到底能走多远系列(30) 扯淡: 30篇啦!从2012-08-15开始的系列,东平西凑将近一年的时间也就这么几篇.目标的100篇,按这个速度也要再搞两年呢. 发博客果然不是件容易的事,怪不得更多的 ...

  2. mysql样例数据库employees

    Oracle和sqlserver都有基于员工信息的样例数据库,同样mysql上也是有的. 给出一个连接地址https://github.com/datacharmer/test_db. 下载后直接调用 ...

  3. UVa 11246 - K-Multiple Free set

    题意大意: 一个{1..n}的集合,求一个子集合,使得元素个数最多,并且不存在有两个元素x * k = y,求出最多的元素个数是多少. 分析: 先要删除k倍的,删除为{k, 2k, 3k, 4k, 5 ...

  4. JQuery中操作Css样式

    //1.获取和设置样式 $("#tow").attr("class")获取ID为tow的class属性 $("#two").attr(&qu ...

  5. Apache配置站点根目录、用户目录及页面访问属性

    一.配置站点根目录及页面访问属性 DocumentRoot "/www/htdoc" <Directory "/www/htdoc"> Option ...

  6. 笨小猴 2008年NOIP全国联赛提高组

    题目描述 Description 笨小猴的词汇量很小,所以每次做英语选择题的时候都很头疼.但是他找到了一种方法,经试验证明,用这种方法去选择选项的时候选对的几率非常大! 这种方法的具体描述如下:假设m ...

  7. Codeforces Round #122 (Div. 2)

    A. Exams 枚举分数为3.4.5的数量,然后计算出2的数量即可. B. Square 相当于求\(\min{x(n+1)\ \%\ 4n=0}\) 打表发现,对\(n\ \%\ 4\)分类讨论即 ...

  8. PHP 将json的stdClass Object转成数组array

    PHP和JS通讯通常都用json,但是PHP要用json的数据,通过json_decode转出来的数组并不是标准的array,所以需要用这个函数进行转换. function object_array( ...

  9. HDU-4747 Mex(线段树区间更新)

    题目大意:给一个长度为n的整数序列,定义mex(i,j)表示区间[i,j]中没有出现过的最小非负整数,求sigma(mex(i,j)),即序列中所有连续非空子区间的mex之和. 题目分析: answe ...

  10. kuangbin_ShortPath D (POJ 3268)

    本来在想 单源多点很好解决但是多源单点怎么解 然后我发现只要倒过来就可以了 把输入存下来然后 处理完dis1 重新init一次 倒着再输入一次 处理dis2 输出max(dis1[i] + dis2[ ...