最长回文子串 ( manacher算法 ) HDU3068

manacher算法

【最长回文子串】

给定一个字符串，求最长的回文子串。回文的意思即字符串关于中心对称。

【引入】

解决这个问题，一般思维是枚举中心，向两边扩展。还要分奇偶，偶数长度的子串关于中缝对称。这种解决方法的时间复杂度为O(n^2)，对于较长的字符串还是不能接受。manacher算法提供了时间复杂度O(n)的解决方案。

【manacher】

在【引入】中提到的枚举中心的方法，加一些优化，即可达到manacher的目的。
1、为了方便起见，我们先处理一下奇偶长度的问题。

我们在给定字符串的所有空隙中都插入一个没出现过的字符，不妨用‘#’

例如，“ab”变为“#a#b#”，“aba”变为“#a#b#a#”。

是不是巧妙的都转化为奇数长度了呢？

设转变后的字符串为s

2、用数组p[i]记录以第i个字符为中心时，最多可以向两边扩展多长。只有关于中心对称才扩展。

s[i]	#	a	#	b	#	a	#
p[i]	1	2	1	4	1	2	1

假设当前在为s[i]计算p[i]值，用变量id记录i之前的右边界最大的下标。即以s[id]对称的回文串的右边界是最大的。看起来跟暴力还是一样。而暴力会重复访问一些子串，我们可以避免这样的事情发生。加一个辅助变量id

for(int i=0;i<strlen(s);i++){

第一步，用已找到的回文串信息，给当前位置i一个可能的最大p[i]值。

考虑当前i位置，在不在已找到的最长回文串内部，即 i < id+p[id]。

若 i < id+p[id]，则p[i]=p[id*2-i]（i关于id对称的位置），否则，令p[i]=1;

第二步，根据已得到的p[i]值，继续向两边扩展，以得到p[i]的正确的最大值。

第三步，更新id值

}

3、对于计算完成的p数组，找到其最大值用ans记录，ans-1即为最长回文子串的长度，可自行验证。

【例题HDU3068】http://acm.hdu.edu.cn/showproblem.php?pid=3068

【代码】

#include<bits/stdc++.h>
using namespace std;
const int MAX=221000;
char str[MAX];//原字符串
char s[MAX*2];//插空串
int p[MAX*2];//上述数组下标均从0开始
int manacher()
{
    int n=strlen(str)*2+1;//插空串长度
    int id=0; //最长回文串的中心
    s[0]='#';
    for(int i=0,j=1;str[i];i++)
    {
        s[j++]=str[i];
        s[j++]='#';
    }
    int ans=0;
    p[0]=1;//初始
    for(int i=1;i<n;i++)//对p[i]计算
    {
        if(i < id+p[id])//i包含在最大回文串中
            p[i]=min(p[id*2-i],id+p[id]-i);//对称点的p值,否则i到p[id]右边界的距离
        else p[i]=1;
        while(i-p[i]>=0&&s[i-p[i]]==s[i+p[i]])p[i]++;//往外扩
        if(id+p[id]<i+p[i])id=i;//比较右边界
        if(ans<p[i])ans=p[i];
    }
    return ans-1;
}
int main()
{
    while(~scanf("%s",str))
        printf("%d\n",manacher());
}