[bzoj1878][SDOI2009][HH的项链] (莫队算法)

Description

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

Input

第一行：一个整数N，表示项链的长度。 第二行：N个整数，表示依次表示项链中贝壳的编号（编号为0到1000000之间的整数）。 第三行：一个整数M，表示HH询问的个数。 接下来M行：每行两个整数，L和R（1 ≤ L ≤ R ≤ N），表示询问的区间。

Output

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

Sample Input

6
1 2 3 4 3 5
3
1 2 
3 5
2 6

Sample Output

2
2
4

HINT

对于20%的数据，N ≤ 100，M ≤ 1000；
对于40%的数据，N ≤ 3000，M ≤ 200000；
对于100%的数据，N ≤ 50000，M ≤ 200000。

Solution

分块大小开n^(2/3)就跪了，乖乖开n^(1/2)

//Kaiba_Seto 20170116
//orz cjkmao
#include <math.h>
#include <stdio.h>
#include <string.h>
#include <algorithm>
#define MaxN 50010
#define MaxM 200010
#define MaxS 1000010
#define MaxBuf 1<<22
#define RG register
#define inline __inline__ __attribute__((always_inline))
#define Blue() (S == T&&(T=(S=B)+fread(B,1,MaxBuf,stdin),S == T)?0:*S++)

char B[MaxBuf],*S=B,*T=B;

inline void Rin(RG int &x) {
    x=0;RG int c=Blue(),f=1;
    for(; c < 48||c > 57; c=Blue())
        if(c == 45)f=-1;
    for(; c > 47&&c < 58; c=Blue())
        x=(x<<1)+(x<<3)+c-48;
    x*=f; }

int n,m,a[MaxN],cnt[MaxS],block_size,ans[MaxM];

struct Request{
    int l,r,id;
    bool operator < (const Request &other)const {
        if(l/block_size == other.l/block_size)
            return r/block_size < other.r/block_size;
        return l/block_size < other.l/block_size; } }Q[MaxM];

inline void extend(RG int &tmp,RG int pos,RG int dir) {
    if(dir == 1) {
        if(++cnt[a[pos]] == 1)tmp++; }
    else
        if(--cnt[a[pos]] == 0)tmp--; }

inline void block_solve() {
    for(RG int l=1,r=0,ans=0,i=1; i<=m; i++) {
        while(l < Q[i].l)extend(ans,l++,-1);
        while(l > Q[i].l)extend(ans,--l,1);
        while(r > Q[i].r)extend(ans,r--,-1);
        while(r < Q[i].r)extend(ans,++r,1);
        :: ans[Q[i].id]=ans; } }

int main() {

    Rin(n);
    for(RG int i=1; i<=n; i++)
        Rin(a[i]);
    Rin(m);
    for(RG int i=1; i<=m; i++)
        Rin(Q[i].l),Rin(Q[i].r),Q[i].id=i;

    //block_size=floor(pow(n,2.0/3)+1);
    block_size=static_cast<int>(sqrt(n));
    std::sort(Q+1,Q+1+m);
    block_solve();

    for(RG int i=1; i<=m; i++)
        printf("%d\n",ans[i]);
    return 0; }