[bzoj2743][HEOI2012]采花(樹狀數組+離線) - Pumbit-Legion－IT工程師數位筆記本

文章出處

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次等

 1 #include<stdio.h>
 2 #include<stdlib.h>
 3 #include<string.h>
 4 #include<algorithm>
 5 #define LL long long
 6 using namespace std;
 7 typedef struct{
 8     int l,r;
 9     int pos;
10 }qry;
11 qry qr[1000010];
12 int a[1000050];
13 int now[1001000],next[1000050];
14 int bit[1000050],n;
15 int ans[1000050];
16 int cmp1(const qry a,const qry b){
17     if(a.l<b.l)return 1;
18     else if(a.l==b.l)return a.r<b.r;
19     else return 0;
20 }
21 int lb(int x){
22     return x&(-x);
23 }
24 int c1(int x){
25     while(x<=n){
26         bit[x]++;
27         x+=lb(x);
28     }
29     return 0;
30 }
31 int c2(int x){
32     while(x<=n){
33         bit[x]--;
34         x+=lb(x);
35     }
36     return 0;
37 }
38 int q(int x){
39     int ans=0;
40     while(x){
41         ans+=bit[x];
42         x-=lb(x);
43     }
44     return ans;
45 }
46 int main(){
47     int m,c;
48     scanf("%d %d %d",&n,&c,&m);
49     for(int i=1;i<=n;i++){
50         scanf("%d",&a[i]);
51     }
52     for(int i=n;i>0;i--){
53         next[i]=now[a[i]];
54         now[a[i]]=i;
55     }
56     for(int i=1;i<=c;i++)if(next[now[i]])c1(next[now[i]]);
57     for(int i=1;i<=m;i++){
58         scanf("%d %d",&qr[i].l,&qr[i].r);
59         qr[i].pos=i;
60     }
61     sort(qr+1,qr+m+1,cmp1);
62     int l=1;
63     for(int i=1;i<=m;i++){
64         //?????
65         while(l<qr[i].l){
66             if(next[l])c2(next[l]);
67             if(next[next[l]])c1(next[next[l]]);
68             l++;
69         }
70         ans[qr[i].pos]=q(qr[i].r)-q(qr[i].l-1);
71     }
72     for(int i=1;i<=m;i++)printf("%d\n",ans[i]);
73     return 0;
74 }

View Code

//ps：后來尹神代碼不科學的比我快了將近兩秒，于是繼續改...

 1 #include<stdio.h>
 2 #include<stdlib.h>
 3 #include<string.h>
 4 #include<algorithm>
 5 #define LL long long
 6 using namespace std;
 7 typedef struct{
 8     int l,r;
 9     int pos;
10 }qry;
11 qry qr[1000010];
12 int a[1000010];
13 int now[1000010],next[1000010];
14 int bit[1000010];
15 int n;
16 int ans[1000010];
17 inline int read(){
18     char ch=getchar();
19     int sum=0;
20     while(ch<'0'||ch>'9')ch=getchar();
21     while(ch<='9'&&ch>='0')sum=sum*10+ch-'0',ch=getchar();
22     return sum;
23 }
24 int cmp1(const qry &a,const qry &b){
25     return a.l<b.l;
26 }
27 inline int lb(int x){
28     return x&(-x);
29 }
30 inline int c1(int x){
31     while(x<=n){
32         ++bit[x];
33         x+=lb(x);
34     }
35     return 0;
36 }
37 inline int c2(int x){
38     while(x<=n){
39         --bit[x];
40         x+=lb(x);
41     }
42     return 0;
43 }
44 inline int q(int x){
45     int ans=0;
46     while(x){
47         ans+=bit[x];
48         x-=lb(x);
49     }
50     return ans;
51 }
52 int main(){
53     register int m,c,i,n1;
54     n=read(),c=read(),m=read();
55      n1=n;
56     for(i=1;i<=n1;++i){
57         a[i]=read();
58     }
59     for(i=n1;i>0;--i){
60         next[i]=now[a[i]];
61         now[a[i]]=i;
62     }
63     for(i=1;i<=c;++i)if(next[now[i]])c1(next[now[i]]);
64     for(i=1;i<=m;++i){
65 
66         qr[i].l=read(),qr[i].r=read();
67         qr[i].pos=i;
68     }
69     sort(qr+1,qr+m+1,cmp1);
70     register int l=1;
71     for(i=1;i<=m;++i){
72         //?????
73         while(l<qr[i].l){
74             if(next[l])c2(next[l]);
75             if(next[next[l]])c1(next[next[l]]);
76             ++l;
77         }
78         ans[qr[i].pos]=q(qr[i].r)-q(qr[i].l-1);
79     }
80     for(i=1;i<=m;++i)printf("%d\n",ans[i]);
81     return 0;
82 }