wannafly 练习赛10 E 数列查找（莫队+分块分块分块......）

链接：https://www.nowcoder.net/acm/contest/58/E

时间限制：C/C++ 1秒，其他语言2秒
空间限制：C/C++ 32768K，其他语言65536K
64bit IO Format: %lld

题目描述

给你一个长为n的序列a，m次查询区间[l,r]内出现次数第k1小的数中值第k2小的数是多少？
保证输入合法

输入描述:

第一行一个数n
第二行n个数表示序列a
第三行一个数m
之后m行每行四个数表示l r k1 k2

输出描述:

对于每次询问输出一行一个数表示答案

输入例子:

10

3 6 6 8 3 10 1 6 5 6

10

4 7 1 2

5 7 1 1

5 6 1 2

2 6 2 1

8 9 1 1

6 9 1 2

1 2 1 1

1 4 2 1

5 7 1 3

2 6 1 3

输出例子:

-->

示例1

输入

10

3 6 6 8 3 10 1 6 5 6

10

4 7 1 2

5 7 1 1

5 6 1 2

2 6 2 1

8 9 1 1

6 9 1 2

1 2 1 1

1 4 2 1

5 7 1 3

2 6 1 3

输出

说明

3 6 6 8 3 10 1 6 5 6
[4,7]中出现1次的有1,3,8,10，第2小的是3
[5,7]中出现1次的有1,3,10，第1小的是1
[5,6]中出现1次的有3,10，第2小的是10
[2,6]中出现2次的有6，第1小的是6
[8,9]中出现1次的有5,6，第1小的是5
[6,9]中出现1次的有1,5,6,10，第2小的是5
[1,2]中出现1次的有3,6，第1小的是3
[1,4]中出现2次的有6，第1小的是6
[5,7]中出现1次的有1,3,10，第3小的是10
[2,6]中出现1次的有3,8,10，第3小的是10

备注:

对于100%的数据，
有1<=n,a[i],m<=40000
数据保证一定能找到那个数

这一题没有码出来，去看了一下别人提交的代码看懂就没自己写了

/////////////////////////////////////////////////

可以考虑进行高维离散化

比如说，对于一个数x，其在序列中出现了y次

开个vector < int > v[ MAXN ]

在v[1] , v[2] ... v[y]中都push_back( x )

然后对于每个vector，分别进行离散化

这样就保证了空间线性

在高维离散化的基础上进行值域分块，然后跑莫队即可

O( nsqrtm + msqrtn ) = O( msqrtn )

/////////////////////////////////////////////////

题解里是高维离散化，但是我看到的代码我感觉那位大佬写的好像更好一点，

他是用了两次分块，次数分块和值分块

#include<cstdio>

#include<algorithm>

#include<cmath>

using namespace std;

int read(){

    int x=,f=;char ch=getchar();

    while (ch<''||ch>'') {if (ch=='-') f=-;ch=getchar();}

    while (''<=ch&&ch<='') x=(x<<)+(x<<)+ch-'',ch=getchar();

    return x*f;

}

const int N=;

int n,m,have[N],a[N],ap[N],cnt[N][],cntk[],ans[N],pos,bb,blo,L,R;

//have[i]表示是否有次数为i的值出现

//ap[i]表示值为i的数字出现了几次

//cntk[i]表示有多少个不同数值出现了k次，k处在i块里

//cnt[i][j]表示出现了i次的值里，值处在j块有多少

struct node{

    int l,r,blo,id,k1,k2;

}q[N];

bool operator < (const node &A,const node &B){

    return A.blo<B.blo||A.blo==B.blo&&A.r<B.r;

}

void work(int x,int w){

    int val=a[x],cc=ap[val];

    if (!--have[cc]) cntk[(cc-)/bb+]--;

    cnt[cc][(val-)/bb+]--;

    ap[val]+=w;cc=ap[val];

    if (++have[cc]==) cntk[(cc-)/bb+]++;

    cnt[cc][(val-)/bb+]++; //[次数][值所属的块]

}

int qry(int k1,int k2){

    int i=,j;

    while (k1>cntk[i]) k1-=cntk[i],i++;//sqrt(n),次数所属的块，每个块有多少个出现了

    for (j=(i-)*bb+;j<=i*bb;j++)//sqrt(n)

        if (have[j])

        if (!--k1) break;

    pos=j;i=;//k1小的是j次/pos次

    while (k2>cnt[pos][i]) k2-=cnt[pos][i],i++;//sqrt(n)

    for (j=(i-)*bb+;j<=i*bb;j++)//sqrt(n)

        if (ap[j]==pos)

        if (!--k2) break;

    return j;

}

int main(){

    n=read();bb=(int)sqrt(n);

    for (int i=;i<=n;i++) a[i]=read();

    m=read();blo=(int)sqrt(m);

    for (int i=;i<=m;i++)

        q[i].l=read(),q[i].r=read(),q[i].k1=read(),

        q[i].k2=read(),q[i].id=i,q[i].blo=(q[i].l-)/blo+;

    sort(q+,q+m+);

    L=;R=;

    for (int i=;i<=m;i++)

    {

        while (L>q[i].l) work(--L,);

        while (R<q[i].r) work(++R,);

        while (L<q[i].l) work(L++,-);

        while (R>q[i].r) work(R--,-);

        ans[q[i].id]=qry(q[i].k1,q[i].k2);

    }

    for (int i=;i<=m;i++) printf("%d\n",ans[i]);

    return ;

}

引用了别的大佬的代码，加了一点注释

分块自己用的还是很生涩，惨

假如是高维离散化的话，应该得跟f题一样用链表存一下出现的次数，然后每个次数里面的值再分块查询，emm就这样