[SDOI2009]HH的项链解题报告

Posted five20

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了[SDOI2009]HH的项链解题报告相关的知识,希望对你有一定的参考价值。

原题目:洛谷P1972

题目描述

HH 有一串由各种漂亮的贝壳组成的项链。HH 相信不同的贝壳会带来好运,所以每次散步完后,他都会随意取出一段贝壳,思考它们所表达的含义。HH 不断地收集新的贝壳,因此,他的项链变得越来越长。有一天,他突然提出了一个问题:某一段贝壳中,包含了多少种不同的贝壳?这个问题很难回答……因为项链实在是太长了。于是,他只好求助睿智的你,来解决这个问题。

输入输出格式

输入格式:

第一行:一个整数N,表示项链的长度。

第二行:N 个整数,表示依次表示项链中贝壳的编号(编号为0 到1000000 之间的整数)。

第三行:一个整数M,表示HH 询问的个数。

接下来M 行:每行两个整数,L 和R(1 ≤ L ≤ R ≤ N),表示询问的区间。

输出格式:

M 行,每行一个整数,依次表示询问对应的答案。

输入输出样例

输入样例#1:
6
1 2 3 4 3 5
3
1 2
3 5
2 6
输出样例#1:
2
2
4

说明

数据范围:

对于20%的数据,N ≤ 100,M ≤ 1000;
对于40%的数据,N ≤ 3000,M ≤ 200000;
对于100%的数据,N ≤ 50000,M ≤ 200000。

Solution:

1、暴力大法20分——30分~~,思路就是随你乱搞,随便能想出复杂度O(mn2),又能想到计数,于是简化到了O(mn),时间复杂度还是很高,40分都水不了:

 1 #include<cstdio>
 2 int a,b,c,d,e,f,i,m,x[100000]={0},w[2000000]={0};
 3 int main()
 4 {
 5 scanf("%d",&a);
 6 c=0;
 7 for(i=1;i<=a;i++)
 8 {
 9     scanf("%d",&b);
10     if(w[b]==0)
11     {
12         c++;
13         w[b]=c;
14         x[i]=c;
15     }
16     else x[i]=w[b];
17 }
18 scanf("%d",&b);
19 for(i=1;i<=b;i++)
20 {
21     scanf("%d%d",&d,&e);
22     f=0;
23     for(m=1;m<=c;m++)
24     {
25         w[m]=0;
26     }
27     for(m=d;m<=e;m++)
28     {
29         if(w[x[m]]==0)
30         {
31             w[x[m]]=1;
32             f++;
33         }
34     }
35     printf("%d\n",f);
36 }
37 }

那好,我们来讲一下正解(不会告诉你其实我问了两位大佬,学了两种方法),但实际上,这题做法很多。

1、莫队算法(从ka哥处知道的神奇算法,据说是莫涛神犇在比赛时想出来的,很NB)。莫队算法就是解决区间查询问题的离线算法利器。(其实很简单,但是貌似涉及到更改就有些麻烦~~)很容易能想到mn的算法,而莫队就是用奇技淫巧将时间复杂度变为了O(M√n),这里我不过多赘述,不懂莫队的去这里。

下面贴一下莫队的代码(其实就是个暴力,注意每一次要先让R指针右移,再移动L指针,不然可能会出现L>R的情况,然后WA)

 1 #include<cstdio>
 2 #include<cmath>
 3 #include<algorithm>
 4 using namespace std;
 5 const int maxn=50005,maxq=200005;
 6 int h[maxn],n,q,a[maxn],ans[maxq];
 7 struct data{
 8     int l,r,id;
 9     bool operator<(const data&b)const{
10         if (h[l]==h[b.l]) return r<b.r;
11         return l<b.l;
12     }
13 }que[maxq];
14 inline int red(){
15     int tot=0,f=1;char ch=getchar();
16     while (ch<0||9<ch) {if (ch==-) f=-f;ch=getchar();}
17     while (0<=ch&&ch<=9) tot=tot*10+ch-48,ch=getchar();
18     return tot*f;
19 }
20 void blocker(){
21     int k=sqrt(n);
22     for (int i=1;i<=n;i++) h[i]=(i-1)/k+1;
23 }
24 int L=0,R=0,now=0,hsh[1000005];
25 void move(int x,int d){
26     if (d==1){
27         if (hsh[a[x]]==0) now++;
28         hsh[a[x]]++;
29     }else{
30         hsh[a[x]]--;
31         if (hsh[a[x]]==0) now--;
32     }
33 }
34 int main(){
35     n=red();
36     for (int i=1;i<=n;i++) a[i]=red();
37     blocker();
38     q=red();
39     for (int i=1;i<=q;i++) que[i].l=red(),que[i].r=red(),que[i].id=i;
40     sort(que+1,que+1+q);
41     for (int i=1;i<=q;i++){
42         while (L<que[i].l) move(L++,-1);
43         while (L>que[i].l) move(--L,1);
44         while (R<que[i].r) move(++R,1);
45         while (R>que[i].r) move(R--,-1);
46         ans[que[i].id]=now;
47     }
48     for (int i=1;i<=q;i++) printf("%d\n",ans[i]);
49     return 0;
50 }

2、树状数组+前缀和(老余告诉我的思路)。 因为操作中只有询问没有更新,所以可以使用离线算法,对所有问按右区间升序排序。以该数字第一次在区间中出现的点代表所有的点。如果是第一次出现,那么该数字 在之前从未出现或上一次出现不再区间内。记录每个位置i的数字的前一个相同数字出现的位置hash[i],没有前一个相同的hash[i]为0。然后从前到后扫描询问,每次将上一个同值点的值加1,然后求当前区间的左界的前缀和就是答案了。将当前位置下个位置的值减1,这样做可以保证任意一个数字在任意一段区间中最多出现一次。

 

 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 int a[50005],s[50005],has[50005],last[50005],n,m;
 4 struct o{
 5 int x,y,num;
 6 }q[200005];
 7 inline int gi()
 8 {
 9     int a=0;char x=getchar();bool f=0;
10     while((x>9||x<0)&&x!=-)x=getchar();
11     if(x==-)x=getchar(),f=1;
12     while(x>=0&&x<=9)a=a*10+x-0,x=getchar();
13     return f?-a:a;
14 }
15 bool cmp(o a,o b)
16 {return a.y<b.y;}
17 inline void add(int k,int a)
18 {
19     while(k<=n)
20     {
21         s[k]+=a;
22         k+=k&-k;
23     }
24 }
25 int ans(int x)
26 {
27     int sum=0;
28     while(x)
29     {
30         sum+=s[x];
31         x-=x&-x;
32     }
33     return sum;
34 }
35 int main()
36 {    
37     n=gi();
38     for(int i=1;i<=n;i++)
39     {int A=gi();has[i]=last[A]+1;last[A]=i;}
40     m=gi();
41     for(int i=1;i<=m;i++)
42     {
43         q[i].x=gi();q[i].y=gi();q[i].num=i;
44     }
45     sort(q+1,q+m+1,cmp);
46     int now=1;
47     for(int i=1;i<=m;i++)
48     {
49         while(now<=q[i].y)
50         {
51             now++;
52             add(has[now-1],1);
53             add(now,-1);
54         }
55         a[q[i].num]=ans(q[i].x);
56     }
57     for(int i=1;i<=m;i++)
58     printf("%d\n",a[i]);
59     return 0;
60 }

3、主席树(可持久化线段树)。关于主席树戳这里造一个lst[i]表示位置i的这个数上一次出现的位置(如果没有就是0) 那么对于每次询问L~R范围里的数字种数 ,其实就是询问对于L<=i<=R,满足lst[i]<L的个数。那么就很好搞了,对lst造N棵值域线段树询问就直接找到L-1这个位置看它前面有几个 。时间复杂度O(Mlogn)

 1 #include<cstdio>
 2 #include<algorithm>
 3 using namespace std;
 4 const int N=2e5+10;
 5 const int M=N*10;
 6 struct President_Tree{
 7     int L,R,sum;
 8 }T[M<<1];
 9 int last[N],top[M],root[N];
10 int n,m,T_cnt=1;
11 int read(){
12     int x=0,f=1;char ch=getchar();
13     while (ch<0 || ch>9){if (ch==-)f=-1;ch=getchar();}
14     while (0<=ch && ch<=9){x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}
15     return x*f;
16 }
17 void insert(int &now,int x,int l=0,int r=n){
18     T[T_cnt++]=T[now];now=T_cnt-1;
19     T[now].sum++;
20     if (l==r)return;
21     int mid=(l+r)>>1;
22     if (x<=mid)insert(T[now].L,x,l,mid);
23         else insert(T[now].R,x,mid+1,r);
24 }
25 int query(int i,int j,int qr,int l=0,int r=n){
26     if (r<=qr)return T[j].sum-T[i].sum;
27     int mid=(l+r)>>1,t=query(T[i].L,T[j].L,qr,l,mid);
28     if (mid<qr)t+=query(T[i].R,T[j].R,qr,mid+1,r);
29     return t;
30 }
31 int main(){
32     n=read();
33     for (int i=1;i<=n;i++){
34         int x=read();
35         last[i]=top[x];
36         top[x]=i;
37     }
38     root[0]=0;
39     for (int i=1;i<=n;i++){
40         root[i]=root[i-1];
41         insert(root[i],last[i]);
42     }
43     m=read();
44     for (int i=1;i<=m;i++){
45         int l=read(),r=read();
46         printf("%d\n",query(root[l-1],root[r],l-1));
47     }
48     return 0;
49 }

 

O(Qlogn)

 

以上是关于[SDOI2009]HH的项链解题报告的主要内容,如果未能解决你的问题,请参考以下文章

[SDOI2009]HH的项链

[SDOI2009]HH的项链

P1972 [SDOI2009]HH的项链

[BZOJ1878][SDOI2009]HH的项链

P1972 [SDOI2009]HH的项链

[BZOJ1878][SDOI2009]HH的项链