[bzoj3809]Gty的二逼妹子序列/[bzoj3236][Ahoi2013]作业

[bzoj3809]Gty的二逼妹子序列/[bzoj3236][Ahoi2013]作业

bzoj   bzoj

题目大意：一个序列，m个询问在$[l,r]$区间的$[x,y]$范围内的数的个数/种类。

思路：

因为看不出来怎么做，所以考虑莫队。

太懒想写个树状数组但是$nsqrt{n}log_{n}$太虚

所以将数分块，修改$O(1)$，查询$O(sqrt{n})$

完结

#include<algorithm>
#include<cstdio>
#include<cmath>
using namespace std;
const int N=114514,SN=366;
template<typename tp>inline void read(tp &kk){
    tp ret=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){ret=ret*10+ch-'0';ch=getchar();}
    kk=ret*f;
}
int n,m,a[N],bl[N],bs;
struct ques
{
    int l,r,x,y,id;
    bool operator < (const ques &a)const{return bl[l]==bl[a.l]?r<a.r:bl[l]<bl[a.l];}
    void init(int i){id=i;read(l),read(r),read(x),read(y);}
}q[N];

struct ans
{
    int a1,a2;
    ans(){a1=a2=0;}
}prt[N];

struct ShangYang
{
    int l[SN],r[SN],buk[SN],bukk[SN],b[N],po[N],sz;
    void start()
    {
        sz=ceil(sqrt(100000));
        int bb=0;
        for(int i=1;i<=n;i+=sz)
        {
            bb++;
            l[bb]=i;
            r[bb]=min(i+sz-1,n);
            for(int j=l[bb];j<=r[bb];j++) b[j]=bb;
        }
    }
    void add(int x)
    {
        if(!po[x]) bukk[b[x]]++;
        po[x]++;
        buk[b[x]]++;
    }
    void mus(int x)
    {
        po[x]--;
        buk[b[x]]--;
        if(!po[x]) bukk[b[x]]--;
    }
    ans query(int x,int y)
    {
        ans ret;
        int bx=b[x],by=b[y];
        if(bx==by)
        {
            for(int i=x;i<=y;i++)
                ret.a1+=po[i],ret.a2+=(!!po[i]);
        }else
        {
            for(int i=bx+1;i<by;i++)
                ret.a1+=buk[i],ret.a2+=bukk[i];
            for(int i=x;i<=r[bx];i++)
                ret.a1+=po[i],ret.a2+=(!!po[i]);
            for(int i=l[by];i<=y;i++)
                ret.a1+=po[i],ret.a2+=(!!po[i]);
        }
        return ret;
    }
}sy;

void icu()
{
    int l=1,r=0;
    for(int i=1;i<=m;i++)
    {
        while(r<q[i].r) sy.add(a[++r]);
        while(l>q[i].l) sy.add(a[--l]);
        while(r>q[i].r) sy.mus(a[r--]);
        while(l<q[i].l) sy.mus(a[l++]);
        prt[q[i].id]=sy.query(q[i].x,q[i].y);
    }
}

int main()
{
    read(n),read(m);
    bs=ceil(sqrt(n));
    for(int i=1;i<=n;i++) read(a[i]),bl[i]=(i-1)/bs+1;
    for(int i=1;i<=m;i++) q[i].init(i);
    // puts("114514");
    sy.start();
    sort(q+1,q+1+m);
    icu();
    // puts("1919810");
    for(int i=1;i<=m;i++) printf("%d %d
",prt[i].a1,prt[i].a2);
    return 0;
}