[BZOJ3809]Gty的二逼妹子序列

题目大意：
　　给定一个长度为$n(nleq10^5)$的数列，$m(mleq10^6)$组询问，每次求出$A_{l..r}$中，权值$in[a,b]$的权值的种类数。

思路：
　　莫队处理询问，分块维护每个区间内的权值种类数。时间复杂度$O(msqrt n)$。如果用树状数组是$O(msqrt nlog n)$的，会TLE。

#include<cmath>
#include<cstdio>
#include<cctype>
#include<algorithm>
inline int getint() {
    register char ch;
    while(!isdigit(ch=getchar()));
    register int x=ch^'0';
    while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
    return x;
}
const int N=100001,M=1000000;
int v[N],ans[M],cnt[N],sum[N],n,block,bel[N],begin[N],end[N];
struct Query {
    int l,r,a,b,id;
    bool operator < (const Query &another) const {
        return l/block<another.l/block||(l/block==another.l/block&&r/block<another.r/block);
    }
};
Query q[M];
inline void ins(const int &x) {
    if(!cnt[v[x]]++) sum[bel[v[x]]]++;
}
inline void del(const int &x) {
    if(!--cnt[v[x]]) sum[bel[v[x]]]--;
}
inline int query(const int &l,const int &r) {
    int ret=0;
    if(bel[l]==bel[r]) {
        for(register int i=l;i<=r;i++) ret+=(bool)cnt[i];
        return ret;
    }
    for(register int i=l;bel[i]==bel[l];i++) ret+=(bool)cnt[i];
    for(register int i=r;bel[i]==bel[r];i--) ret+=(bool)cnt[i];
    for(register int i=bel[l]+1;i<bel[r];i++) ret+=sum[i];
    return ret;
}
int main() {
    block=sqrt(n=getint());
    const int m=getint();
    for(register int i=1;i<=n;i++) {
        v[i]=getint();
        bel[i]=i/block;
        if(!begin[bel[i]]) begin[bel[i]]=i;
        end[bel[i]]=i;
    }
    for(register int i=0;i<m;i++) {
        const int l=getint(),r=getint(),a=getint(),b=getint();
        q[i]=(Query){l,r,a,b,i};
    }
    std::sort(&q[0],&q[m]);
    for(register int i=0,l=1,r=0;i<m;i++) {
        while(r<q[i].r) ins(++r);
        while(l>q[i].l) ins(--l);
        while(r>q[i].r) del(r--);
        while(l<q[i].l) del(l++);
        ans[q[i].id]=query(q[i].a,q[i].b);
    }
    for(register int i=0;i<m;i++) printf("%d
",ans[i]);
    return 0;
}