bzoj1878: [SDOI2009]HH的项链

Description

HH有一串由各种漂亮的贝壳组成的项链。HH相信不同的贝壳会带来好运，所以每次散步 完后，他都会随意取出一段贝壳，思考它们所表达的含义。HH不断地收集新的贝壳，因此他的项链变得越来越长。有一天，他突然提出了一个问题：某一段贝壳中，包含了多少种不同的贝壳？这个问题很难回答。。。因为项链实在是太长了。于是，他只好求助睿智的你，来解决这个问题。

Input

第一行：一个整数N，表示项链的长度。

第二行：N个整数，表示依次表示项链中贝壳的编号（编号为0到1000000之间的整数）。

第三行：一个整数M，表示HH询问的个数。

接下来M行：每行两个整数，L和R（1 ≤ L ≤ R ≤ N），表示询问的区间。

N ≤ 50000，M ≤ 200000

Output

M行，每行一个整数，依次表示询问对应的答案。

Sample Input

6

1 2 3 4 3 5

3

1 2

3 5

2 6

Sample Output

2

2

4

学习莫队算法之后手写的第一道题，，感觉莫队其实不是很难，，，思路非常清晰暴力，只是时间复杂度计算的很玄学。

这题大致题意就是在给你M个询问，让你计算L~R之间有多少个不相同的数，

我记得某位大佬说过，只要是可以离线的区间查询问题，用莫队算法都会产生奇效，

那么这道题，有多组区间查询，还不能用线段树ST表什么乱七八糟的东西，，那我们就用莫队分块来搞一搞

大体思路就是对查询进行分块，以左端点所在的块为关键字排序，左端点在同一块上的查询就按照右端点进行排序，

那么我们固定一个当前区间，枚举排好序的查询中每组查询的左右端点，在移动的过程中便更新完了答案

看起来非常暴力是不是，，，莫队算法最核心的部分就是最开始对查询进行的排序，

我们把数组分为k块；根据我们的排序可知，我们的左端点最多移动(n*q)/k，右端点最多移动 n*k次，

所以说当我们把块的大小设为sqrt（n）时能保证时间复杂度稳定为n*sqrt（n），

其他就没什么好说的了，，看代码吧

#include<cstdio>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
struct jgt{int l,r,k,ans;}qst[1000005];
int len,n,q,num[1000005],cot[1000005],bel[1000005];
int ll,rr,ansl;
bool cmp1(jgt a,jgt b){return bel[a.l]<bel[b.l]||(bel[a.l]==bel[b.l]&&a.r<b.r);}//以左端点所在的块为第一关键字，右端点的位置为第二关键字排序，
bool cmp2(jgt a,jgt b){return a.k<b.k;}
void del(int a)//从当前区间删除元素
{
cot[a]--;
if(cot[a]==0)ansl--;
}
void add(int a)//向当前区间添加元素
{
if(cot[a]==0)ansl++;
cot[a]++;
}
int main()
{
scanf("%d",&n);
len=ceil(sqrt(1.0*n));//确定区间大小

for(int i=1;i<=n;i++)
{
scanf("%d",&num[i]);
bel[i]=(i-1)/len+1;//确定所属的块
}
scanf("%d",&q);
for(int i=1;i<=q;i++)//输入查询
{
qst[i].k=i;
scanf("%d%d",&qst[i].l,&qst[i].r);
}
sort(qst+1,qst+q+1,cmp1);//排序

ll=qst[1].l;rr=ll;add(num[ll]);//确定当前区间

for(int i=1;i<=q;i++)
{
while(ll>qst[i].l) {ll--;add(num[ll]);}//移动左端点
while(ll<qst[i].l) {del(num[ll]);ll++;}

while(rr>qst[i].r) {del(num[rr]);rr--;}//移动右端点
while(rr<qst[i].r) {rr++;add(num[rr]);}

qst[i].ans=ansl;//储存答案
}

sort(qst+1,qst+q+1,cmp2);

for(int i=1;i<=q;i++)   printf("%d\n",qst[i].ans);

return 0;
}