[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[1000010];
int a[1000050];
int now[1001000],next[1000050];
int bit[1000050],n;
int ans[1000050];
int cmp1(const qry a,const qry b){
    if(a.l<b.l)return 1;
    else if(a.l==b.l)return a.r<b.r;
    else return 0;
}
int lb(int x){
    return x&(-x);
}
int c1(int x){
    while(x<=n){
        bit[x]++;
        x+=lb(x);
    }
    return 0;
}
int c2(int x){
    while(x<=n){
        bit[x]--;
        x+=lb(x);
    }
    return 0;
}
int q(int x){
    int ans=0;
    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=1;i<=n;i++){
        scanf("%d",&a[i]);
    }
    for(int i=n;i>0;i--){
        next[i]=now[a[i]];
        now[a[i]]=i;
    }
    for(int i=1;i<=c;i++)if(next[now[i]])c1(next[now[i]]);
    for(int i=1;i<=m;i++){
        scanf("%d %d",&qr[i].l,&qr[i].r);
        qr[i].pos=i;
    }
    sort(qr+1,qr+m+1,cmp1);
    int l=1;
    for(int i=1;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-1);
    }
    for(int i=1;i<=m;i++)printf("%d
",ans[i]);
    return 0;
}

//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[1000010];
int a[1000010];
int now[1000010],next[1000010];
int bit[1000010];
int n;
int ans[1000010];
inline int read(){
    char ch=getchar();
    int sum=0;
    while(ch<'0'||ch>'9')ch=getchar();
    while(ch<='9'&&ch>='0')sum=sum*10+ch-'0',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 0;
}
inline int c2(int x){
    while(x<=n){
        --bit[x];
        x+=lb(x);
    }
    return 0;
}
inline int q(int x){
    int ans=0;
    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=1;i<=n1;++i){
        a[i]=read();
    }
    for(i=n1;i>0;--i){
        next[i]=now[a[i]];
        now[a[i]]=i;
    }
    for(i=1;i<=c;++i)if(next[now[i]])c1(next[now[i]]);
    for(i=1;i<=m;++i){

        qr[i].l=read(),qr[i].r=read();
        qr[i].pos=i;
    }
    sort(qr+1,qr+m+1,cmp1);
    register int l=1;
    for(i=1;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-1);
    }
    for(i=1;i<=m;++i)printf("%d
",ans[i]);
    return 0;
}