BZOJ 3809 Gty的二逼妹子序列

暴力。。。严格的说是mlogn+m√nlogn的。。。就是莫队+BIT。

然后我们可以按照值域分块，这样修改O（1）查询O（√n），这样总复杂度就是2*m*√n。

有点像做了一个时间复杂度的权衡。。。看起来查询变慢了，事实上总复杂度是变快了。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 100050
#define maxm 1000050
using namespace std;
int n,m,len=300,blo[maxn],st[maxn],ret=0,val[maxn],cnt[maxn],a[maxn],ans[maxm],p1,p2;
int a1,a2,a3,a4;
struct query
{
    int l,r,a,b,id;
    query (int l,int r,int a,int b,int id):l(l),r(r),a(a),b(b),id(id) {}
    query () {}
}q[maxm];
int read()
{
    char ch;int data=0;
    while (ch<'0' || ch>'9') ch=getchar();
    while (ch>='0' && ch<='9')
    {
        data=data*10+ch-'0';
        ch=getchar();
    }
    return data;
}
bool cmp(query x,query y)
{
    if (blo[x.l]==blo[y.l]) return x.r<y.r;
    return blo[x.l]<blo[y.l];
}
void build()
{
    ret=1;
    for (int i=1;i<=n;i++)
    {
        blo[i]=ret;
        if (!(i%len)) {st[ret]=i;if (i!=n) ret++;}
    }
    if (n%len) st[ret]=n;
}
void modify(int x,int y)
{
    bool flag=cnt[x];
    cnt[x]+=y;
    if ((!cnt[x]) && (flag)) val[blo[x]]--;
    if ((cnt[x]) && (!flag)) val[blo[x]]++;
}
int ask(int l,int r)
{
    int ret1=0,ret2=0;
    for (int i=1;i<=blo[l-1]-1;i++) ret1+=val[i];
    for (int i=st[blo[l-1]-1]+1;i<=l-1;i++) ret1+=(cnt[i]>0);
    for (int i=1;i<=blo[r]-1;i++) ret2+=val[i];
    for (int i=st[blo[r]-1]+1;i<=r;i++) ret2+=(cnt[i]>0);
    return ret2-ret1;
}
int main()
{
    n=read();m=read();
    for (int i=1;i<=n;i++) a[i]=read();
    for (int i=1;i<=m;i++)
    {
        scanf("%d%d%d%d",&a1,&a2,&a3,&a4);
        q[i]=query(a1,a2,a3,a4,i);
    }
    build();
    sort(q+1,q+m+1,cmp);
    for (int i=1;i<=m;i++)
    {
        for (;p1>q[i].l;p1--) modify(a[p1-1],1);
        for (;p1<q[i].l;p1++) modify(a[p1],-1);
        for (;p2>q[i].r;p2--) modify(a[p2],-1);
        for (;p2<q[i].r;p2++) modify(a[p2+1],1);
        ans[q[i].id]=ask(q[i].a,q[i].b);
    }
    for (int i=1;i<=m;i++) printf("%d\n",ans[i]);
    return 0;
}