BZOJ3809: Gty的二逼妹子序列

【传送门：BZOJ3809】

---

简要题意：

　　给出一个n个数的序列，序列中的数范围为1到n，给出m个操作，每个操作输入l,r,a,b，输出l到r中权值为a到b的数的种类

---

题解：

　　肯定不能在线求，离线就用莫队

　　一开始想法是用树状数组维护权值总类，结果T了

　　旁边的lxj大佬D了我，说：这道题卡了树状数组，因为它修改带log

　　认真地听了一波讲解，就直接把权值也给分块了

　　因为权值和数的个数的区间是一样的，所以可以一并分块

　　然后修改O(1)，而求值用n^0.5，直接卡了过去

　　伏地膜lxj大佬

---

参考代码：

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
struct node
{
    int l,r,a,b,id,d;
}q[1100000];
int s[110000];
int bl[1100],br[1100],belong[110000];
bool cmp1(node n1,node n2)
{
    if(belong[n1.l]<belong[n2.l]) return true;
    if(belong[n1.l]>belong[n2.l]) return false;
    if(belong[n1.l]==belong[n2.l])
    {
        if(n1.r<n2.r) return true;
        if(n1.r>n2.r) return false;
    }
    return false;
}
bool cmp2(node n1,node n2)
{
    return n1.id<n2.id;
}
int d[1100],sum[110000];
void add(int x)
{
    sum[x]++;
    if(sum[x]==1) d[belong[x]]++;
}
void del(int x)
{
    sum[x]--;
    if(sum[x]==0) d[belong[x]]--;
}
int solve(int x,int y)
{
    int bx=belong[x],by=belong[y];
    int ans=0;
    if(bx==by)
    {
        for(int i=x;i<=y;i++)
        {
            if(sum[i]>0) ans++;
        }
        return ans;
    }
    for(int i=bx+1;i<=by-1;i++) ans+=d[i];
    for(int i=x;i<=br[bx];i++) if(sum[i]>0) ans++;
    for(int i=bl[by];i<=y;i++) if(sum[i]>0) ans++;
    return ans;
}
int main()
{
    int n,m;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&s[i]);
    int block=int(sqrt(n));
    for(int i=1;i<=n;i++)
    {
        int t=(i-1)/block+1;
        belong[i]=t;
        if(bl[t]==0) br[t-1]=i-1,bl[t]=i;
    }
    br[belong[n]]=n;
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d%d%d",&q[i].l,&q[i].r,&q[i].a,&q[i].b);
        q[i].id=i;
    }
    sort(q+1,q+m+1,cmp1);
    memset(sum,0,sizeof(sum));
    memset(d,0,sizeof(d));
    int l=1,r=0;
    for(int i=1;i<=m;i++)
    {
        while(l>q[i].l){l--;add(s[l]);}
        while(l<q[i].l){del(s[l]);l++;}
        while(r>q[i].r){del(s[r]);r--;}
        while(r<q[i].r){r++;add(s[r]);}
        q[i].d=solve(q[i].a,q[i].b);
    }
    sort(q+1,q+m+1,cmp2);
    for(int i=1;i<=m;i++) printf("%d\n",q[i].d);
    return 0;
}