BZOJ1878 [SDOI2009] HH的项链
题目链接: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 }