[BZOJ3920]Yuuna的礼物

题目大意：
　　给你一个长度为$n(n\le40000)$的数列$\{a_i\}(1\le a_i\le n)$，给出$m(m\le40000)$次询问，每次给出$l,r,k_1,k_2$询问区间$[l,r]$中出现次数第$k_1$小的数中第$k_2$小的数是多少？

思路：
　　运用莫队算法离线处理所有询问，分块维护数列中每个数的出现次数。考虑如何维护出现次数的出现次数以及出现次数相同的数。同样采用分块，先预处理出对于某一种出现次数，所有可能的数，再将其离散化，对于离散化后的数分块维护。由于题目空间限制只有24M，因此需要手写内存池或者使用vector，可以证明离散化以后的数组空间复杂度是$O(n)$的。这样时间复杂度$O(m\sqrt n)$，空间复杂度是$O(n+m)$。

 #include<cmath>
 #include<cstdio>
 #include<cctype>
 #include<vector>
 #include<algorithm>
 inline int getint() {
     register char ch;
     while(!isdigit(ch=getchar()));
     register int x=ch^'';
     while(isdigit(ch=getchar())) x=(((x<<)+x)<<)+(ch^'');
     return x;
 }
 const int N=,M=;
 int n,a[N],b[N],block,ans[M],cnt[N];
 struct Query {
     int l,r,k1,k2,id;
     bool operator < (const Query &another) const {
         return l/block<another.l/block||(l/block==another.l/block&&r<another.r);
     }
 };
 Query q[M];
 std::vector<int> v[N],rank[N],cnt1,sum1,cnt2[N],sum2[N];
 inline void modify(const int &x,const int &p,const int &d) {
     cnt2[x][p]+=d;
     sum2[x][p/block]+=d;
     if(d==&&!cnt1[x]++) sum1[x/block]++;
     if(d==-&&!--cnt1[x]) sum1[x/block]--;
 }
 inline void modify(const int &x,const int &d) {
     if(cnt[x]) modify(cnt[x],rank[x][cnt[x]],-);
     if(cnt[x]+=d) modify(cnt[x],rank[x][cnt[x]],);
 }
 inline int solve(const std::vector<int> &c,const std::vector<int> &s,const int &k) {
     register int x=,cnt=;
     while((cnt+s[x])<k) cnt+=s[x++];
     for(register int i=x*block;i<(x+)*block;i++) {
         if((cnt+=!!c[i])>=k) return i;
     }
 }
 inline int query(const int &k1,const int &k2) {
     const int x=solve(cnt1,sum1,k1);
     return v[x][solve(cnt2[x],sum2[x],k2)];
 }
 int main() {
     block=sqrt(n=getint());
     for(register int i=;i<=n;i++) {
         cnt[a[i]=b[i]=getint()]++;
     }
     std::sort(&b[],&b[n]+);
     for(register int i=;i<=n;i++) {
         if(b[i]==b[i-]) continue;
         rank[b[i]].resize(cnt[b[i]]+);
         for(register int j=;j<=cnt[b[i]];j++) {
             rank[b[i]][j]=v[j].size();
             v[j].push_back(b[i]);
         }
         cnt[b[i]]=;
     }
     for(register int i=;i<=n;i++) {
         cnt2[i].resize(v[i].size()+);
         sum2[i].resize(cnt2[i].size()/block+);
     }
     const int m=getint();
     for(register int i=;i<m;i++) {
         const int l=getint(),r=getint(),k1=getint(),k2=getint();
         q[i]=(Query){l,r,k1,k2,i};
     }
     cnt1.resize(n+);
     sum1.resize(n/block+);
     std::sort(&q[],&q[m]);
     for(register int i=,l=,r=;i<m;i++) {
         while(r<q[i].r) modify(a[++r],);
         while(l>q[i].l) modify(a[--l],);
         while(r>q[i].r) modify(a[r--],-);
         while(l<q[i].l) modify(a[l++],-);
         ans[q[i].id]=query(q[i].k1,q[i].k2);
     }
     for(register int i=;i<m;i++) {
         printf("%d\n",ans[i]);
     }
     return ;
 }