[BZOJ]3809 3809: Gty的二逼妹子序列 莫队+分块

3809: Gty的二逼妹子序列

Time Limit: 80 Sec   Memory Limit: 28 MB
Submit: 2082   Solved: 628
[ Submit][ Status][ Discuss]

Description

Autumn和Bakser又在研究Gty的妹子序列了！但他们遇到了一个难题。

对于一段妹子们，他们想让你帮忙求出这之内美丽度∈[a,b]的妹子的美丽度的种类数。

为了方便，我们规定妹子们的美丽度全都在[1,n]中。

给定一个长度为n(1<=n<=100000)的正整数序列s(1<=si<=n)，对于m(1<=m<=1000000)次询问“l,r,a,b”，每次输出sl...sr中，权值∈[a,b]的权值的种类数。

Input

第一行包括两个整数n,m(1<=n<=100000,1<=m<=1000000),表示数列s中的元素数和询问数。

第二行包括n个整数s1...sn(1<=si<=n)。

接下来m行,每行包括4个整数l,r,a,b(1<=l<=r<=n,1<=a<=b<=n),意义见题目描述。

保证涉及的所有数在C++的int内。

保证输入合法。

Output

对每个询问，单独输出一行，表示sl...sr中权值∈[a,b]的权值的种类数。

Sample Input

10 10
4 4 5 1 4 1 5 1 2 1
5 9 1 2
3 4 7 9
4 4 2 5
2 3 4 7
5 10 4 4
3 9 1 1
1 4 5 9
8 9 3 3
2 2 1 6
8 9 1 4

Sample Output

2
0
0
2
1
1
1
0
1
2

HINT

样例的部分解释：

5 9 1 2

子序列为4 1 5 1 2

在[1,2]里的权值有1,1,2，有2种，因此答案为2。

3 4 7 9

子序列为5 1

在[7,9]里的权值有5，有1种，因此答案为1。

4 4 2 5

子序列为1

没有权值在[2,5]中的，因此答案为0。

2 3 4 7

子序列为4 5

权值在[4,7]中的有4,5，因此答案为2。

建议使用输入/输出优化。

Source

[ Submit][ Status][ Discuss] 

HOME Back

      树状数组听说要T...要不然就是非常简单的莫队用值域树状数组维护区间和查询即可. 因为莫队是m根号n的时间， 每移动一步就要修改， 那么树状数组修改时log， 查询也是log， 那么时间复杂度就是m根号n * log值域.

      肯定要降移动的复杂度. 那么比log还小的就是o(1)的了. 也就是说我们要o(1)的修改， 那么什么线段树平衡树等log级别的数据结构就都可以滚了， 要用一种比较骨骼惊奇的数据结构--分块.分块往往会往根号的复杂度想， 但是这里莫队每一步移动是单点修改， 那么修改就是o(1)的. 查询还是根号的 区间中间整块直接查询， 两端暴力算即可. 移动o(1)， 查询根号， 复杂度降为 O(m根号n). 

#include<stdio.h>
#include<cmath>
#include<algorithm>
using namespace std;
const int maxn = 100005;
int bloc, n, m;
int c[maxn], cnt[maxn], cntt[maxn], blo[maxn], ans[maxn * 10];
inline const int read(){
    register int x = 0;
    register char ch = getchar();
    while(ch < '0' || ch > '9') ch = getchar();
    while(ch >= '0' && ch <= '9') x = (x << 3) + (x << 1) + ch - '0', ch = getchar();
    return x;
}
struct point{
    int l, r, a, b, id;
}q[maxn * 10];
inline bool cmp(point x, point y){
    return (blo[x.l] == blo[y.l]) ? x.r < y.r : blo[x.l] < blo[y.l];
}
inline void ins(int x){
    cnt[x]++;
    if(cnt[x] == 1) cntt[blo[x]]++;
}
inline void del(int x){
    cnt[x]--;
    if(!cnt[x]) cntt[blo[x]]--;
}
inline int query(int l, int r){
    int ed = 0;
    if(blo[r] - blo[l] < 2){
        for(int i = l; i <= r; ++i)
            if(cnt[i]) ++ed;
        return ed;
    }
    for(int i = blo[l] + 1; i < blo[r]; ++i)
        ed += cntt[i];
    for(int i = l; blo[i] == blo[l]; ++i)
        if(cnt[i]) ++ed;
    for(int i = r; blo[i] == blo[r]; --i)
        if(cnt[i]) ++ed;
    return ed;
}
inline void Captain_Mo(){
    for(register int l = 1, r = 0, i = 1; i <= m; ++i){
        while(r < q[i].r) ins(c[++r]);
        while(r > q[i].r) del(c[r--]);
        while(l < q[i].l) del(c[l++]);
        while(l > q[i].l) ins(c[--l]);
        ans[q[i].id] = query(q[i].a, q[i].b);
    }
}
int main(){
    n = read(), m = read(), bloc = (int)sqrt(n);
    for(register int i = 1; i <= n; ++i) c[i] = read();
    for(register int i = 1; i <= m; ++i) q[i].l = read(), q[i].r = read(), q[i].a = read(), q[i].b = read(), q[i].id = i;
    for(register int i = 1; i <= n; ++i) blo[i] = (i - 1) / bloc + 1;
    sort(q + 1, q + m + 1, cmp);
    Captain_Mo();
    for(register int i = 1; i <= m; ++i) printf("%d\n", ans[i]);
}