【BZOJ 3809】Gty的二逼妹子序列 莫队+分块

比较裸，没什么好说的。只是，其实注意一下常数还是可以卡过去的，而且即使是莫队加分块依然可以卡T，所以，就当练练分块把，既然出题人都想这么考了。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#define maxn 100021
#define Q(a) (a>0 ? 1 : 0)
using namespace std;
void read(int& x){
x=0;char c=getchar();
for(;c>'9'||c<'0';c=getchar());
for(;c>='0'&&c<='9';c=getchar())x=x*10+c-'0';
}
int n,m,block,bl[maxn],cnt[maxn],s[maxn],ans[1000021],val[maxn],rs[maxn],ls[maxn];
struct que{
int l,r,id,a,b;
bool operator<(const que& y)const{
return bl[l]==bl[y.l] ? r<y.r : bl[l] < bl[y.l];
}
}q[1000021];
inline void del(int x){cnt[x]--;if(cnt[x]==0)s[bl[x]]--;}
inline void add(int x){cnt[x]++;if(cnt[x]==1)s[bl[x]]++;}
inline int query(int a,int b){
int ans=0;
if(bl[a]==bl[b]){
for(int i=a;i<=b;i++)ans+=Q(cnt[i]);
}else{
for(int i=bl[a]+1;i<=bl[b]-1;i++)ans+=s[i];
for(int i=a;i<=rs[a];i++)ans+=Q(cnt[i]);
for(int i=ls[b];i<=b;i++)ans+=Q(cnt[i]);
}
return ans;
}

void solve(){
int l=1,r=0;
for(int i=1;i<=m;i++){
while(l<q[i].l)del(val[l++]);
while(l>q[i].l)add(val[--l]);
while(r<q[i].r)add(val[++r]);
while(r>q[i].r)del(val[r--]);
ans[q[i].id]=query(q[i].a,q[i].b);
}
}

int main(){
read(n),read(m);block=sqrt(n);
for(int i=1;i<=n;i++){
bl[i]=(i-1)/block+1;
rs[i]=bl[i]*block;ls[i]=(bl[i]-1)*block+1;
}
for(int i=1;i<=n;i++)read(val[i]);
for(int i=1;i<=m;i++){read(q[i].l),read(q[i].r),read(q[i].a),read(q[i].b);q[i].id=i;}
sort(q+1,q+1+m);
solve();
for(int i=1;i<=m;i++)printf("%d\n",ans[i]);
return 0;
}