BZOJ1878 [SDOI2009] HH的项链

Posted 2020-06-17

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=1878

Description

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

Input

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

Output

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

第一眼感觉是莫队，看了PoPoQQQ也写了莫队……

然而网上的题解铺天盖地的树状数组？

卧槽世界上竟然有如此巧妙的离线做法，我居然没想到

 1 #include <iostream>
 2 #include <cstdio>
 3 #include <algorithm>
 4 #include <cstring>
 5 #define rep(i,l,r) for(int i=l; i<=r; i++)
 6 #define clr(x,y) memset(x,y,sizeof(x))
 7 using namespace std;
 8 const int maxn = 50010;
 9 inline int read(){
10     int ans = 0, f = 1;
11     char c = getchar();
12     for(; !isdigit(c); c = getchar())
13     if (c == ‘-‘) f = -1;
14     for(; isdigit(c); c = getchar())
15     ans = ans * 10 + c - ‘0‘;
16     return ans * f;
17 }
18 struct Query{
19     int l,r,id,ans;
20 }q[200010];
21 int n,m,a[maxn],t[maxn],next[maxn],p[1000010],mx=-1;
22 inline int lowbit(int x){
23     return x & (-x);
24 }
25 inline bool cmp1(Query x,Query y){
26     if (x.l == y.l) return x.r < y.r; return x.l < y.l;
27 }
28 inline bool cmp2(Query x,Query y){
29     return x.id < y.id;
30 }
31 inline void update(int x,int v){
32     for (int i = x; i <= n; i += lowbit(i)) t[i] += v;
33 }
34 inline int query(int x){
35     int ret = 0;
36     for(int i = x; i; i -= lowbit(i)) ret += t[i];
37     return ret;
38 }
39 int main(){
40     n = read();
41     rep(i,1,n) a[i] = read(), mx = max(mx,a[i]);
42     for(int i=n; i; i--) next[i] = p[a[i]], p[a[i]] = i;
43     rep(i,1,mx) if (p[i]) update(p[i],1);
44     m = read();
45     rep(i,1,m) q[i].l = read(), q[i].r = read(), q[i].id = i;
46     sort(q+1,q+m+1,cmp1);
47     int l = 1;
48     rep(i,1,m){
49         while (l < q[i].l){
50             if (next[l]) update(next[l],1);
51             l++;
52         }
53         q[i].ans = query(q[i].r) - query(q[i].l - 1);
54     }
55     sort(q+1,q+m+1,cmp2);
56     rep(i,1,m) printf("%d\n",q[i].ans);
57     return 0;
58 }

View Code