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

AC通道：http://www.lydsy.com/JudgeOnline/problem.php?id=3809

【题解】

如果没有区间[a,b]的限制，那么这题就是一道模板题。

查询答案的时候我们需要得到区间[a,b]的答案，很容易想到树状数组。

但是树状数组的时间复杂度是O(n*logn*√n)，不足以卡过此题。

所以我们就用到了分块大法，时间复杂度O(n*√n)/*************
bzoj 3809
by chty
2016.11.29
*************/
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<ctime>
#include<cmath>
#include<algorithm>
using namespace std;
#define FILE "read"
#define MAXN 100010
#define MAXM 1000010
#define up(i,j,n) for(int i=j;i<=n;i++)
#define down(i,j,n) for(int i=j;i>=n;i--)
namespace INIT{
char buf[1<<15],*fs,*ft;
inline char getc(){return (fs==ft&&(ft=(fs=buf)+fread(buf,1,1<<15,stdin),fs==ft))?0:*fs++;}
inline int read(){
int x=0,f=1; char ch=getc();
while(!isdigit(ch)) {if(ch=='-') f=-1; ch=getc();}
while(isdigit(ch)) {x=x*10+ch-'0'; ch=getc();}
return x*f;
}
}using namespace INIT;
struct node{int x,y,a,b,id;}q[MAXM];
int n,m,block,a[MAXN],ans[MAXM],cnt[MAXN],vis[MAXN],blo[MAXN];
bool cmp(node a,node b) {return a.x/block==b.x/block?a.y<b.y:a.x/block<b.x/block;}
void add(int x){
vis[x]^=1;
blo[x/block]+=vis[x]?1:-1;
}
void updata(int x,int v){
if(!cnt[x]) add(x);
cnt[x]+=v;
if(!cnt[x]) add(x);
}
int query(int x,int y){
int l=x/block,r=y/block,sum(0);
if(l==r){
up(i,x,y) sum+=vis[i];
return sum;
}
up(i,l+1,r-1) sum+=blo[i];
up(i,x,(l+1)*block-1) sum+=vis[i];
up(i,r*block,y) sum+=vis[i];
return sum;
}
void solve(){
sort(q+1,q+m+1,cmp);
int l(1),r(0);
up(i,1,m){
while(q[i].y>r) updata(a[++r],1);
while(q[i].y<r) updata(a[r--],-1);
while(q[i].x<l) updata(a[--l],1);
while(q[i].x>l) updata(a[l++],-1);
ans[q[i].id]=query(q[i].a,q[i].b);
}
up(i,1,m) printf("%d\n",ans[i]);
}
void init(){
n=read(); m=read(); block=(int)sqrt(n*1.0);
up(i,1,n) a[i]=read();
up(i,1,m) q[i].x=read(),q[i].y=read(),q[i].a=read(),q[i].b=read(),q[i].id=i;
}
int main(){
freopen(FILE".in","r",stdin);
freopen(FILE".out","w",stdout);
init();
solve();
return 0;
}