码迷,mamicode.com
首页 > 其他好文 > 详细

hihocoder 1419 重复旋律4

时间:2018-03-06 17:03:48      阅读:208      评论:0      收藏:0      [点我收藏+]

标签:include   ring   weight   iostream   coder   code   output   最大   ref   

描述

小Hi平时的一大兴趣爱好就是演奏钢琴。我们知道一个音乐旋律被表示为长度为 N 的数构成的数列。小Hi在练习过很多曲子以后发现很多作品中的旋律有重复的部分。

我们把一段旋律称为(k,l)-重复的,如果它满足由一个长度为l的字符串重复了k次组成。 如旋律abaabaabaaba是(4,3)重复的,因为它由aba重复4次组成。

小Hi想知道一部作品中k最大的(k,l)-重复旋律。

解题方法提示

输入

一行一个仅包含小写字母的字符串。字符串长度不超过 100000。

输出

一行一个整数,表示答案k。

Sample Input
babbabaabaabaabab
Sample Output
4
假如循环节长度为L
如果求出了后缀i和i+L的lcp
那么重复次数k=lcp(i,i+L)/L+1
那么O(n^2)的暴力有了,枚举L和i,计算k更新最大值
但是枚举i可以只枚举L的倍数,算出k
对于i和i+L,最多只会相差1
如果最优串的开始位置恰好在L的倍数上,那我们找到的最大的k就是正确答案
如果不在L的倍数上,那么最优的开始位置肯定在(i-L,i)上
如果lcp模L有余数,说明前面还需要补上L-lcp%L位,如果可行答案可以加1
于是可以判断开始位置为i-L+lcp%L与i+lcp%L的lcp
如果构成重复,那么算出来的
k=lcp(i-L+lcp%L,i+lcp%L)/L+1=[lcp(i,i+L)/L+1]+1,比开始位置为i要好
于是复杂度:
$O(\frac{n}{1}+\frac{n}{2}+\frac{n}{3}+...\frac{n}{n})$
$=O(nlogn)$
 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<algorithm>
 5 #include<cmath>
 6 using namespace std;
 7 int n,m,c[200001],x[200001],y[200001],SA[200001],s[200001],rk[200001],h[200001];
 8 int Min[200001][21],Log[200001],ans;
 9 char ch[200001];
10 void radix_sort()
11 {int i;
12   for (i=0;i<=m;i++)
13     c[i]=0;
14   for (i=1;i<=n;i++)
15     c[x[y[i]]]++;
16   for (i=2;i<=m;i++)
17     c[i]+=c[i-1];
18   for (i=n;i>=1;i--)
19     SA[c[x[y[i]]]--]=y[i];
20 }
21 void build_SA()
22 {int i,j,k,p;
23   for (i=1;i<=n;i++)
24     x[i]=s[i],y[i]=i;
25   m=100000;
26   radix_sort();
27   for (k=1;k<=n;k<<=1)
28     {
29       p=0;
30       for (i=n-k+1;i<=n;i++)
31     y[++p]=i;
32       for (i=1;i<=n;i++)
33     if (SA[i]>k) y[++p]=SA[i]-k;
34       radix_sort();
35       p=1;swap(x,y);
36       x[SA[1]]=1;
37       for (i=2;i<=n;i++)
38     x[SA[i]]=((y[SA[i]]==y[SA[i-1]])&&(y[SA[i]+k]==y[SA[i-1]+k]))?p:++p;
39       if (p>=n) break;
40       m=p;
41     }
42   for (i=1;i<=n;i++)
43     rk[SA[i]]=i;
44   int L=0;
45   for (i=1;i<=n;i++)
46     {
47       if (L>0) L--;
48       j=SA[rk[i]-1];
49       while (i+L<=n&&j+L<=n&&(s[j+L]==s[i+L])) L++;
50       h[rk[i]]=L;
51     }
52 }
53 int rmq(int x,int y)
54 {
55   int L=Log[y-x+1];
56   return min(Min[x][L],Min[y-(1<<L)+1][L]);
57 }
58 int lcp(int x,int y)
59 {
60   if (rk[x]>rk[y]) swap(x,y);
61   return rmq(rk[x]+1,rk[y]);
62 }
63 int main()
64 {int i,j,L;
65   scanf("%s",ch);
66   n=strlen(ch);
67   for (i=1;i<=n;i++)
68     {
69       s[i]=ch[i-1]-a+1;
70     }
71   build_SA();
72   memset(Min,127/2,sizeof(Min));
73   for (i=1;i<=n;i++)
74     Min[i][0]=h[i];
75   Log[1]=0;
76   for (i=2;i<=n;i++)
77     Log[i]=Log[i/2]+1;
78   for (j=1;(1<<j)<=n;j++)
79     {
80       for (i=1;i<=n-(1<<j)+1;i++)
81     Min[i][j]=min(Min[i][j-1],Min[i+(1<<j-1)][j-1]);
82     }
83   for (L=1;L<=n;L++)
84     {
85       for (i=1;i+L<=n;i+=L)
86     {
87       int p=lcp(i,i+L);
88       ans=max(ans,p/L+1);
89       if (i-L+p%L>=1)
90         {
91           ans=max(ans,lcp(i-L+p%L,i+p%L)/L+1);
92         }
93     }
94     }
95   printf("%d\n",ans);
96 }

 

hihocoder 1419 重复旋律4

标签:include   ring   weight   iostream   coder   code   output   最大   ref   

原文地址:https://www.cnblogs.com/Y-E-T-I/p/8514661.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!