[AHOI2013]作业 & Gty的二逼妹子序列 莫队

～～～题面～～～

题解：

　　题目要求统计一个区间内数值在[a, b]内的数的个数和种数，而这个是可以用树状数组统计出来的，所以可以考虑莫队。

　　考虑区间[l, r]转移到[l, r + 1],那么对于维护个数的树状数组就直接加即可。

　　对于维护种数的树状数组，我们额外维护一个数组num，表示数a在区间内出现了多少次，如果是新出现的，那么就加入树状数组。
　　如果要删除一个数并且这个数在区间内只出现了一次，那么就删除这个数。注意不论什么情况都要实时维护num数组。

　　然后莫队即可。

　　Gty的二逼妹子序列是洛谷P4867和作业的某一问是一模一样的，把数组开大点就可以过。

 #include<bits/stdc++.h>
 using namespace std;
 #define R register int
 #define AC 100100

 int n, m, cnt, block;
 int ans1[AC], ans2[AC], s[AC], num[AC], tot[AC];

 struct node{
     int l, r, a, b, id;
 }q[AC];

 inline int lowbit(int x)
 {
     return x & (-x);
 }

 struct kkk{
     int a[AC];

     void add(int x, int y)
     {
         for(R i = x; i <= cnt; i += lowbit(i)) a[i] += y;
     }

     int find(int x)
     {
         int rnt = ;
         for(R i = x; i; i -= lowbit(i)) rnt += a[i];
         return rnt;
     }
 }c1, c2;

 inline int read()
 {
     int x = ;char c = getchar();
     while(c > '' || c < '') c = getchar();
     while(c >= '' && c <= '') x = x *  + c - '', c = getchar();
     return x;
 }

 int half(int x)//查询离散化后的值
 {
     int l = , r = cnt, mid;
     while(l < r)
     {
         mid = (l + r) >> ;
         if(num[mid] == x) return mid;
         else if(num[mid] > x) r = mid;
         else l = mid + ;
     }
     return l;
 }

 int half1(int x)//查询最小的大于等于a的值
 {
     int l = , r = cnt, mid;
     if(num[cnt] < x) return cnt + ;
     while(l < r)
     {
         mid = (l + r) >> ;
         if(num[mid] >= x) r = mid;
         else l = mid + ;
     }
     return l;
 }

 int half2(int x)//查询最大的小于等于b的值
 {
     int l = , r = cnt, mid;
     if(num[] > x) return ;
     while(l < r)
     {
         mid = (l + r + ) >> ;//强制偏右
         if(num[mid] > x) r = mid - ;
         else l = mid;
     }
     return l;
 }

 inline bool cmp(node a, node b)
 {
     if(a.l / block != b.l / block) return a.l < b.l;
     else return a.r < b.r;//分块排序
 }

 void pre()
 {
     n = read(), m = read(), block = sqrt(n);
     for(R i = ; i <= n; i ++) s[i] = num[i] = read();
     sort(num + , num + n + );
     for(R i = ; i <= n; i ++)
         if(num[i] != num[i + ]) num[++cnt] = num[i];
     for(R i = ; i <= n; i ++) s[i] = half(s[i]);//在这里离散化，这样后面就不用调用了
     for(R i = ; i <= m; i ++)
     {
         q[i].l = read(), q[i].r = read(), q[i].id = i;
         q[i].a = half1(read()), q[i].b = half2(read());
     }
     sort(q + , q + m + , cmp);
 }

 void add(int x)
 {
     c1.add(x, );
     if(!tot[x]) c2.add(x, );
     ++ tot[x];
 }

 void del(int x)
 {
     c1.add(x, -);
     -- tot[x];
     if(!tot[x]) c2.add(x, -);
 }

 void work()//这里每个点的贡献与区间无关，相对独立，所以不用考虑顺序问题
 {
     int l, r;
     for(R i = q[].l; i <= q[].r; i ++) add(s[i]);
     ans1[q[].id] = c1.find(q[].b) - c1.find(q[].a - );
     ans2[q[].id] = c2.find(q[].b) - c2.find(q[].a - );
     l = q[].l, r = q[].r;
     for(R i = ; i <= m; i ++)
     {
         int ll = q[i].l, rr = q[i].r;
         while(ll < l) -- l, add(s[l]);
         while(ll > l) del(s[l]), ++ l;
         while(rr > r) ++ r, add(s[r]);
         while(rr < r) del(s[r]), -- r;
         ans1[q[i].id] = c1.find(q[i].b) - c1.find(q[i].a - );
         ans2[q[i].id] = c2.find(q[i].b) - c2.find(q[i].a - );
     }
     for(R i = ; i <= m; i ++) printf("%d %d\n", ans1[i], ans2[i]);
 }

 int main()
 {
 //    freopen("in.in", "r", stdin);
     pre();
     work();
 //    fclose(stdin);
     return ;
 }