hihocoder #1419 : 后缀数组四·重复旋律4

http://hihocoder.com/problemset/problem/1419

这个题可以用类似https://www.cnblogs.com/hehe54321/p/8808226.html的方法A掉,分析比题解要简单...(应该吧?)

 1 #include<cstdio>
 2 #include<algorithm>
 3 #include<cstring>
 4 #include<vector>
 5 using namespace std;
 6 typedef long long LL;
 7 typedef unsigned long long ULL;
 8 char s[100010];
 9 LL n,anss;
10 ULL powx[100100],hs[100100];
11 const ULL X=131;
12 ULL ghash(LL x,LL y)//[x,y]的哈希
13 {
14     return hs[x]-hs[y+1]*powx[y-x+1];
15 }
16 LL lcp(LL x,LL y)//后缀x和后缀y的lcp
17 {
18     if(x>n||y>n)    return 0;
19     LL l=0,r=min(n-x+1,n-y+1)+1,mid;
20     while(r-l>1)
21     {
22         mid=l+(r-l)/2;
23         if(ghash(x,x+mid-1)==ghash(y,y+mid-1))    l=mid;
24         else    r=mid;
25     }
26     return l;
27 }
28 LL lcs(LL x,LL y)//前缀x和前缀y的lcs
29 {
30     if(!x||!y)    return 0;
31     LL l=0,r=min(x,y)+1,mid;
32     while(r-l>1)
33     {
34         mid=l+(r-l)/2;
35         if(ghash(x-mid+1,x)==ghash(y-mid+1,y))    l=mid;
36         else    r=mid;
37     }
38     return l;
39 }
40 LL ans;
41 int main()
42 {
43     LL i,j,b,l,r,k;
44     scanf("%s",s+1);n=strlen(s+1);
45     powx[0]=1;
46     for(i=1;i<=100000;i++)    powx[i]=powx[i-1]*X;
47     for(i=n;i>=1;i--)    hs[i]=hs[i+1]*X+s[i];
48     for(i=1;i<=n;i++)
49     {
50         b=n/i;
51         for(j=1,k=0;j<=b;j++)
52         {
53             for(;k<b&&ghash(i*(j-1)+1,i*j)==ghash(i*k+1,i*(k+1));k++);
54             //例如对于dba,bca,bca,bde
55             //当i=3,j=2时,k=3
56             l=min(i-1,lcs(i*(j-1),i*k));r=min(i-1,lcp(i*(j-1)+1,i*k+1));
57             //此时l=1,r=1
58             ans=max(ans,k-j+1);
59             //此时只可能分成两段的倍数,且长度为k-j+1
60             //又例如:当i=3,j=3,k=3时
61             //l=2,r=1
62             //k-j+1=1,因此没有需要统计的
63             if(l+r>=i)    ans=max(ans,k-j+2);
64             //但事实上有cabcab没有统计
65             //可以证明,像这种情况漏统计的串的循环次数最多比k-j+1大1
66             //漏统计的数量是从第i*(j-1)-l+1位开始的到第i*(j-1)+r-i+1位开始的
67             //共出现r+l-i+1次,每一次均循环k-j+2次
68             //显然只有当r+l-i+1大于0时才会有此情况
69             //又例如:例如ac,cc,c,当i=2,j=2,k=2,l=1,r=1时
70             //此时还有cccc被漏统计
71             //从第2位开始,到第2位结束,重复2次
72         }
73     }
74     printf("%lld",ans);
75     return 0;
76 }
原文地址:https://www.cnblogs.com/hehe54321/p/8808281.html