首页 > 代码库 > [hdu3068]最长回文(Manacher算法)

[hdu3068]最长回文(Manacher算法)

http://acm.hdu.edu.cn/showproblem.php?pid=3068

题目大意:求最长回文串的长度。

解题关键:Manacher算法

引用一个较好的解释

p[i] = mx > i ? min(p[2 * id - i], mx - i) : 1;

 

可以这么说,这行要是理解了,那么马拉车算法基本上就没啥问题了,那么这一行代码拆开来看就是

如果mx > i, 则 p[i] = min(p[2 * id - i], mx - i)

否则, p[i] = 1

当 mx - i > P[j] 的时候,以S[j]为中心的回文子串包含在以S[id]为中心的回文子串中,由于 i 和 j 对称,以S[i]为中心的回文子串必然包含在以S[id]为中心的回文子串中,所以必有 P[i] = P[j],见下图。

 

技术分享

 


当 P[j] >= mx - i 的时候,以S[j]为中心的回文子串不一定完全包含于以S[id]为中心的回文子串中,但是基于对称性可知,下图中两个绿框所包围的部分是相同的,也就是说以S[i]为中心的回文子串,其向右至少会扩张到mx的位置,也就是说 P[i] >= mx - i。至于mx之后的部分是否对称,就只能老老实实去匹配了。



技术分享


对于 mx <= i 的情况,无法对 P[i]做更多的假设,只能P[i] = 1,然后再去匹配了。

 

 1 #include<bits/stdc++.h> 2 using namespace std; 3 typedef long long ll; 4 char s[110010*2]; 5 int p[110010*2]; 6 int main(){ 7     while(scanf("%s",s)!=EOF){ 8         int len=strlen(s); 9         for(int i=len;i>=0;i--){10             s[2*i+1]=#;11             s[2*i+2]=s[i];12         }13         s[0]=*;14         int id=0,maxlen=0;15         for(int i=2;i<2*len+1;i++){16             if(p[id]+id>i) p[i]=min(p[2*id-i],p[id]+id-i);17             else p[i]=1;18             while(s[i-p[i]]==s[i+p[i]])++p[i];19             if(p[i]+i>p[id]+id) id=i;20             if(p[i]>maxlen) maxlen=p[i];21         }22         printf("%d\n",maxlen-1);23     }24 }

 

[hdu3068]最长回文(Manacher算法)