[bzoj2743][HEOI2012]采花(树状数组+离线)

2743: [HEOI2012]采花

Time Limit: 15 Sec  Memory Limit: 128 MB
Submit: 1832  Solved: 954
[Submit][Status][Discuss]

Description

萧芸斓是Z国的公主，平时的一大爱好是采花。

今天天气晴朗，阳光明媚，公主清晨便去了皇宫中新建的花园采花。花园足够大，容纳了n朵花，花有c种颜色（用整数1-c表示），且花是排成一排的，以便于公主采花。公主每次采花后会统计采到的花的颜色数，颜色数越多她会越高兴！同时，她有一癖好，她不允许最后自己采到的花中，某一颜色的花只有一朵。为此，公主每采一朵花，要么此前已采到此颜色的花，要么有相当正确的直觉告诉她，她必能再次采到此颜色的花。由于时间关系，公主只能走过花园连续的一段进行采花，便让女仆福涵洁安排行程。福涵洁综合各种因素拟定了m个行程，然后一一向你询问公主能采到多少朵花（她知道你是编程高手，定能快速给出答案！），最后会选择令公主最高兴的行程（为了拿到更多奖金！）。

Input

 第一行四个空格隔开的整数n、c以及m。接下来一行n个空格隔开的整数，每个数在[1, c]间，第i个数表示第i朵花的颜色。接下来m行每行两个空格隔开的整数l和r（l ≤ r），表示女仆安排的行程为公主经过第l到第r朵花进行采花。

Output

 

共m行，每行一个整数，第i个数表示公主在女仆的第i个行程中能采到的花的颜色数。

Sample Input

5 3 5
1 2 2 3 1
1 5
1 2
2 2
2 3
3 5

Sample Output

2
0 0 1 0
【样例说明】
询问[1, 5]：公主采颜色为1和2的花，由于颜色3的花只有一朵，公主不采；询问[1, 2]：颜色1和颜色2的花均只有一朵，公主不采；
询问[2, 2]：颜色2的花只有一朵，公主不采；
询问[2, 3]：由于颜色2的花有两朵，公主采颜色2的花；
询问[3, 5]：颜色1、2、3的花各一朵，公主不采。

HINT

【数据范围】

对于100%的数据，1 ≤ n ≤    10^6，c ≤ n，m ≤10^6。

Source

这道题乍一看很像bzoj1878

但其实要更难

(请先参阅bzoj1878  http://www.cnblogs.com/Pumbit-Legion/p/5874133.html)

上一道题其实是弱化版本

“求区间内一个及以上的颜色数量”

而这道题是“两个及以上”

如何变形是一个问题

对于每个要加入的元素，可以判断它下一个元素是否存在，存在再加一

同理，要出去的话，先把因为这个元素而加入的元素删掉，再看看删掉的能不能带来新的元素

因此，这个问题还能推广到3次，4次等

 #include<stdio.h>
 #include<stdlib.h>
 #include<string.h>
 #include<algorithm>
 #define LL long long
 using namespace std;
 typedef struct{
     int l,r;
     int pos;
 }qry;
 qry qr[];
 int a[];
 int now[],next[];
 int bit[],n;
 int ans[];
 int cmp1(const qry a,const qry b){
     if(a.l<b.l)return ;
     else if(a.l==b.l)return a.r<b.r;
     else return ;
 }
 int lb(int x){
     return x&(-x);
 }
 int c1(int x){
     while(x<=n){
         bit[x]++;
         x+=lb(x);
     }
     return ;
 }
 int c2(int x){
     while(x<=n){
         bit[x]--;
         x+=lb(x);
     }
     return ;
 }
 int q(int x){
     int ans=;
     while(x){
         ans+=bit[x];
         x-=lb(x);
     }
     return ans;
 }
 int main(){
     int m,c;
     scanf("%d %d %d",&n,&c,&m);
     for(int i=;i<=n;i++){
         scanf("%d",&a[i]);
     }
     for(int i=n;i>;i--){
         next[i]=now[a[i]];
         now[a[i]]=i;
     }
     for(int i=;i<=c;i++)if(next[now[i]])c1(next[now[i]]);
     for(int i=;i<=m;i++){
         scanf("%d %d",&qr[i].l,&qr[i].r);
         qr[i].pos=i;
     }
     sort(qr+,qr+m+,cmp1);
     int l=;
     for(int i=;i<=m;i++){
         //?????
         while(l<qr[i].l){
             if(next[l])c2(next[l]);
             if(next[next[l]])c1(next[next[l]]);
             l++;
         }
         ans[qr[i].pos]=q(qr[i].r)-q(qr[i].l-);
     }
     for(int i=;i<=m;i++)printf("%d\n",ans[i]);
     return ;
 }

//ps：后来尹神代码不科学的比我快了将近两秒，于是继续改...

 #include<stdio.h>
 #include<stdlib.h>
 #include<string.h>
 #include<algorithm>
 #define LL long long
 using namespace std;
 typedef struct{
     int l,r;
     int pos;
 }qry;
 qry qr[];
 int a[];
 int now[],next[];
 int bit[];
 int n;
 int ans[];
 inline int read(){
     char ch=getchar();
     int sum=;
     while(ch<''||ch>'')ch=getchar();
     while(ch<=''&&ch>='')sum=sum*+ch-'',ch=getchar();
     return sum;
 }
 int cmp1(const qry &a,const qry &b){
     return a.l<b.l;
 }
 inline int lb(int x){
     return x&(-x);
 }
 inline int c1(int x){
     while(x<=n){
         ++bit[x];
         x+=lb(x);
     }
     return ;
 }
 inline int c2(int x){
     while(x<=n){
         --bit[x];
         x+=lb(x);
     }
     return ;
 }
 inline int q(int x){
     int ans=;
     while(x){
         ans+=bit[x];
         x-=lb(x);
     }
     return ans;
 }
 int main(){
     register int m,c,i,n1;
     n=read(),c=read(),m=read();
      n1=n;
     for(i=;i<=n1;++i){
         a[i]=read();
     }
     for(i=n1;i>;--i){
         next[i]=now[a[i]];
         now[a[i]]=i;
     }
     for(i=;i<=c;++i)if(next[now[i]])c1(next[now[i]]);
     for(i=;i<=m;++i){

         qr[i].l=read(),qr[i].r=read();
         qr[i].pos=i;
     }
     sort(qr+,qr+m+,cmp1);
     register int l=;
     for(i=;i<=m;++i){
         //?????
         while(l<qr[i].l){
             if(next[l])c2(next[l]);
             if(next[next[l]])c1(next[next[l]]);
             ++l;
         }
         ans[qr[i].pos]=q(qr[i].r)-q(qr[i].l-);
     }
     for(i=;i<=m;++i)printf("%d\n",ans[i]);
     return ;
 }