bzoj1878: [SDOI2009]HH的项链

传送门
大力莫队算法一点不虚。
我们先将序列分成sqrt（N）块
然后将询问按照左端点坐在快编号为第一关键字，右端点为第二关键字进行排序。
然后我们只需要在两次询问之间进行修改就可以了
可以证明时间复杂度是O（Nsqrt（N））的。

#include<cstring>  
#include<cmath>   
#include<cstdio>  
#include<iostream>  
#include<cstdlib>   
#include<algorithm>
using namespace std;
struct data{
    int l,r,ans,id;
}e[200005];
int b[1000005],a[100005],c[100005],d[100005];
int n,ma,l,m;
inline int read(){
    int k=0;
    char ch=getchar();
    for (;ch<'0'||ch>'9';ch=getchar());
    for (;ch>='0'&&ch<='9';ch=getchar()) k=k*10+ch-48;
    return k;
}
inline bool cmp1(data a,data b){
    if (a.l==b.l) return a.r<b.r;
    return a.l<b.l;
}
inline bool cmp2(data a,data b){return a.id<b.id;}
inline void change(int x,int y){
    for (;x<=n;x+=x&(-x)) d[x]+=y;
}
inline int ask(int x){
    int s=0;
    for (;x>=1;x-=x&(-x)) s+=d[x];
    return s;
}
int main(){
    n=read();
    for (int i=1;i<=n;i++){
        a[i]=read();
        ma=max(ma,a[i]);
    }
    for (int i=n;i>=1;i--){
        c[i]=b[a[i]];
        b[a[i]]=i;
    }
    m=read();
    for (int i=1;i<=m;i++){
        e[i].l=read();
        e[i].r=read();
        e[i].id=i;
    }
    sort(e+1,e+m+1,cmp1);
    for (int i=1;i<=ma;i++)
        if (b[i]) change(b[i],1);
    l=1;
    for (int i=1;i<=m;i++){
        while (l<e[i].l){
            if (c[l]) change(c[l],1);
            l++;
        }
        e[i].ans=ask(e[i].r)-ask(e[i].l-1);
    }
    sort(e+1,e+m+1,cmp2);
    for (int i=1;i<=m;i++) printf("%d\n",e[i].ans);
}