bzoj 3809: Gty的二逼妹子序列

3809: Gty的二逼妹子序列

Time Limit: 80 Sec  Memory Limit: 28 MB

Description

Autumn和Bakser又在研究Gty的妹子序列了！但他们遇到了一个难题。

 

对于一段妹子们，他们想让你帮忙求出这之内美丽度∈[a,b]的妹子的美丽度的种类数。

 

为了方便，我们规定妹子们的美丽度全都在[1,n]中。

 

给定一个长度为n(1<=n<=100000)的正整数序列s(1<=si<=n)，对于m(1<=m<=1000000)次询问“l,r,a,b”，每次输出sl...sr中，权值∈[a,b]的权值的种类数。

Input

第一行包括两个整数n,m(1<=n<=100000,1<=m<=1000000),表示数列s中的元素数和询问数。

 

第二行包括n个整数s1...sn(1<=si<=n)。

 

接下来m行,每行包括4个整数l,r,a,b(1<=l<=r<=n,1<=a<=b<=n),意义见题目描述。

 

保证涉及的所有数在C++的int内。

 

保证输入合法。

Output

对每个询问，单独输出一行，表示sl...sr中权值∈[a,b]的权值的种类数。

Sample Input

10 10
4 4 5 1 4 1 5 1 2 1
5 9 1 2
3 4 7 9
4 4 2 5
2 3 4 7
5 10 4 4
3 9 1 1
1 4 5 9
8 9 3 3
2 2 1 6
8 9 1 4

Sample Output

2
0
0
2
1
1
1
0
1
2

HINT

样例的部分解释：

 

5 9 1 2

子序列为4 1 5 1 2

在[1,2]里的权值有1,1,2，有2种，因此答案为2。

 

3 4 7 9

子序列为5 1

在[7,9]里的权值有5，有1种，因此答案为1。

 

4 4 2 5

子序列为1

没有权值在[2,5]中的，因此答案为0。

 

2 3 4 7

子序列为4 5

权值在[4,7]中的有4,5，因此答案为2。

 

建议使用输入/输出优化。

 Solution:

  这里对于数字进行分块，然后我们将询问根据左端点排序后，就可以用莫队在很可观的的时间复杂度内算出答案了

 本蒟蒻打出的第一个莫队（+分块），how amazing.

 算法的魅力阿 以及就是不写读入优化，你咬我呀

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
const int maxn=100020;
const int maxm=1000020;
struct haha
{
    int l,r,a,b,pos;
}q[maxm];
int a[maxn],ans[maxm],bel[maxn],c[maxn],sum[maxn],n,m;
bool cmp(const haha lu,const haha d) {return bel[lu.l] == bel[d.l] ? lu.r < d.r : lu.l < d.l;}
void inti()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&a[i]);
    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].pos=i;
    int block;
    block=sqrt(n+0.5);
    for(int i=1;i<=n;i++) bel[i]=(i-1)/block+1;
}
int query(int a,int b)
{
    int l=bel[a],r=bel[b],tot=0;
    if(l==r)
    {
        for(int i=a;i<=b;i++)
        if(c[i]>0) tot++;
        return tot;
    }
    for(int i=a;bel[i]==bel[a];i++)
     if(c[i]>0) tot++;
    for(int i=b;bel[i]==bel[b];i--)
     if(c[i]>0) tot++;
    for(int i=bel[a]+1;i<bel[b];i++)
    tot+=sum[i];
    return tot;
}
void updata(int x,int flag)
{
    c[x]+=flag;
    if(flag==-1&&c[x]==0) sum[bel[x]]--;
    if(flag==1&&c[x]==1) sum[bel[x]]++;
}
void work()
{
    sort(q+1,q+m+1,cmp);
    int l,r,ll,rr; l=1;r=0;
    for(int i=1;i<=m;i++)
    {
        ll=q[i].l;rr=q[i].r;
        while(l<ll) updata(a[l++],-1);
        while(l>ll) updata(a[--l],1);
        while(r<rr) updata(a[++r],1);
        while(r>rr) updata(a[r--],-1);
        ans[q[i].pos]=query(q[i].a,q[i].b);
    }
}
int main()
{
    inti();
    work();
    for(int i=1;i<=m;i++)
    printf("%d
",ans[i]);
    return 0;
}