Manacher算法 O(n) 求最长回文子串
2014-03-19 23:35
375 查看
转自:http://bbs.dlut.edu.cn/bbstcon.php?board=Competition&gid=23474
其实原文说得是比较清楚的,只是英文的,我这里写一份中文的吧。
首先:大家都知道什么叫回文串吧,这个算法要解决的就是一个字符串中最长的回文子串有多长。这个算法可以在O(n)的时间复杂度内既线性时间复杂度的情况下,求出以每个字符为中心的最长回文有多长,
这个算法有一个很巧妙的地方,它把奇数的回文串和偶数的回文串统一起来考虑了。这一点一直是在做回文串问题中时比较烦的地方。这个算法还有一个很好的地方就是充分利用了字符匹配的特殊性,避免了大量不必要的重复匹配。
算法大致过程是这样。先在每两个相邻字符中间插入一个分隔符,当然这个分隔符要在原串中没有出现过。一般可以用‘#’分隔。这样就非常巧妙的将奇数长度回文串与偶数长度回文串统一起来考虑了(见下面的一个例子,回文串长度全为奇数了),然后用一个辅助数组P记录以每个字符为中心的最长回文串的信息。P[id]记录的是以字符str[id]为中心的最长回文串,当以str[id]为第一个字符,这个最长回文串向右延伸了P[id]个字符(包含str[id]在内的P[id]个字符)。
原串: w aa bwsw f d
新串: # w # a # a # b # w # s # w # f # d #
辅助数组P: 1 2 1 2 3 2 1 2 1 2 1 4 1 2 1 2 1 2 1
这里有一个很好的性质,P[id]-1就是该回文子串在原串中的长度(不包括‘#’)。如果这里不是特别清楚,可以自己拿出纸来画一画,自己体会体会。当然这里可能每个人写法不尽相同,不过我想大致思路应该是一样的吧。
好,我们继续。现在的关键问题就在于怎么在O(n)时间复杂度内求出P数组了。只要把这个P数组求出来,最长回文子串就可以直接扫一遍得出来了。
由于这个算法是线性从前往后扫的。那么当我们准备求P[i]的时候,i以前的P[j]我们是已经得到了的。我们用mx记在i之前的回文串中,延伸至最右端的位置。同时用id这个变量记下取得这个最优mx时的id值。(注:为了防止字符比较的时候越界,我在这个加了‘#’的字符串之前还加了另一个特殊字符‘$’,故我的新串下标是从1开始的)
好,到这里,我们可以先贴一份代码了。
View Code
其实原文说得是比较清楚的,只是英文的,我这里写一份中文的吧。
首先:大家都知道什么叫回文串吧,这个算法要解决的就是一个字符串中最长的回文子串有多长。这个算法可以在O(n)的时间复杂度内既线性时间复杂度的情况下,求出以每个字符为中心的最长回文有多长,
这个算法有一个很巧妙的地方,它把奇数的回文串和偶数的回文串统一起来考虑了。这一点一直是在做回文串问题中时比较烦的地方。这个算法还有一个很好的地方就是充分利用了字符匹配的特殊性,避免了大量不必要的重复匹配。
算法大致过程是这样。先在每两个相邻字符中间插入一个分隔符,当然这个分隔符要在原串中没有出现过。一般可以用‘#’分隔。这样就非常巧妙的将奇数长度回文串与偶数长度回文串统一起来考虑了(见下面的一个例子,回文串长度全为奇数了),然后用一个辅助数组P记录以每个字符为中心的最长回文串的信息。P[id]记录的是以字符str[id]为中心的最长回文串,当以str[id]为第一个字符,这个最长回文串向右延伸了P[id]个字符(包含str[id]在内的P[id]个字符)。
原串: w aa bwsw f d
新串: # w # a # a # b # w # s # w # f # d #
辅助数组P: 1 2 1 2 3 2 1 2 1 2 1 4 1 2 1 2 1 2 1
这里有一个很好的性质,P[id]-1就是该回文子串在原串中的长度(不包括‘#’)。如果这里不是特别清楚,可以自己拿出纸来画一画,自己体会体会。当然这里可能每个人写法不尽相同,不过我想大致思路应该是一样的吧。
好,我们继续。现在的关键问题就在于怎么在O(n)时间复杂度内求出P数组了。只要把这个P数组求出来,最长回文子串就可以直接扫一遍得出来了。
由于这个算法是线性从前往后扫的。那么当我们准备求P[i]的时候,i以前的P[j]我们是已经得到了的。我们用mx记在i之前的回文串中,延伸至最右端的位置。同时用id这个变量记下取得这个最优mx时的id值。(注:为了防止字符比较的时候越界,我在这个加了‘#’的字符串之前还加了另一个特殊字符‘$’,故我的新串下标是从1开始的)
好,到这里,我们可以先贴一份代码了。
#include<cstdio> #include<cstring> const int M = 110010*2; char str[M];//start from index 1 int p[M]; char s[M]; int n; void checkmax(int &ans,int b){ if(b>ans) ans=b; } inline int min(int a,int b){ return a<b?a:b; } void kp(){ int i; int mx = 0; int id; for(i=1; i<n; i++){ if( mx > i ) p[i] = min( p[2*id-i], p[id]+id-i ); else p[i] = 1; for(; str[i+p[i]] == str[i-p[i]]; p[i]++) ; if( p[i] + i > mx ) { mx = p[i] + i; id = i; } } } void pre() { int i,j,k; n = strlen(s); str[0] = '$'; str[1] = '#'; for(i=0;i<n;i++) { str[i*2 + 2] = s[i]; str[i*2 + 3] = '#'; } n = n*2 + 2; str = 0; } void pt() { int i; int ans = 0; for(i=0;i<n;i++) checkmax(ans, p[i]); printf("%d\n", ans-1); } int main() { int T,_=0; while( scanf("%s", s) !=EOF ) { pre(); kp(); pt(); } return 0; }
View Code
相关文章推荐
- 最长回文(manacher算法,最大回文子串)
- [51Nod](1089)最长回文子串 V2 ---- Manacher算法
- hdu【最长回文子串——Manacher算法】POJ 3974/HDU 3068
- Manacher算法求解最长回文子串
- 51nod1089(最长回文子串之manacher算法)
- Manacher算法——找字符串最长的回文子串
- Manacher算法------求最长回文子串(Java)
- 最长回文子串,Manacher算法
- 最长回文子串 manacher算法
- [hdu 3068] Manacher算法O(n)最长回文子串
- Manacher算法: O(n)时间求字符串的最长回文子串
- hihoCoder 1032 最长回文子串 (manacher算法)
- hdu5371 最长回文子串变形(Manacher算法)
- KT学算法(三)——最长回文子串与Manacher算法
- 求最长回文子串的Manacher算法
- HiHo 1032 最长回文子串 (Manacher算法求解)
- 练习:最长回文子串(Manacher算法)
- Poj3974 最长回文子串 Manacher算法
- 最长回文子串的Manacher算法
- 51nod_1089 最长回文子串 V2(Manacher算法)