最长回文(Manacher 马拉车算法模板解析)

给出一个只由小写英文字符a,b,c...y,z组成的字符串S,求S中最长回文串的长度. 
回文就是正反读都是一样的字符串,如aba, abba等

Input

输入有多组case,不超过120组,每组输入为一行小写英文字符a,b,c...y,z组成的字符串S 
两组case之间由空行隔开(该空行不用处理) 
字符串长度len <= 110000

Output

每一行一个整数x,对应一组case,表示该组case的字符串中所包含的最长回文长度. 

Sample Input

aaaa

abab

Sample Output

4
3

题意:求一个字符串的最长回文串。

思路:马拉车算法模板。

马拉车算法:对于任何一个回文串,可能是奇数个(aba),也可能是偶数个(abba)。对于奇数特别好处理,枚举每一个字符为中点,两边找即可。对于偶数个的话就比较麻烦了。为了简化问题,我们把字符串都转化一下,每个字符前后前面都加一个无关的字符.

eg:   aba    -->  #a#b#a# (新串为奇数个)

        abba   --> #a#b#b#a#(新串为奇数个)

其实原串后面还有一个  '\0'   为了简化问题(不特判结尾)我们在最开始再加一个无关字符 $  

eg:    aba    -->  $#a#b#a#0

         abba   --> $#a#b#b#a#0

然后我们就来找回文串了,

我们定义Mp[ i ] 是以 i 为中心的回文串的半径 (Mp[i] = 1 代表这个回文串只有自己)。

id 是当前已知最长回文串的中点,  mx 是当前已知最长回文串的右端点

最长回文(Manacher 马拉车算法模板解析)_第1张图片

对于当前位置 i 它的回文半径怎么算呢?

因为我们是顺序遍历的,所以到达 i 之前的所有 Mp[ i ]都已经得到了。所有我们可以利用 i ,j 关于ID对称的特性来减少“暴力”求解的次数(当然是直求长度了)。

当 mx  > i 时

      ①.  mx - i > p[ j ] :因为上图蓝色部分已经是回文串了,也就是 i 的最小半径应该和 j 的半径相等。(Mp[ i ] >= Mp[ j ])

                  所以最小半径为 p [ j ]

      ②.   mx - i < p[ j ] :黄色部分的长度小于 p [ j ] 那黄色部分再往右就不知道了, 所以最小半径为 mx - i

当 mx <= i 时,就没能利用 对称点 j 了,所以最小半径为 1

接下来就是暴力循环扩增 Mp[ i ]了

最后只需要求Mp的最大值即可,虽然Mp记录的半径,但由于字符串扩增了,所有Mp就是答案。

代码如下:

#include
#include
#include
#define N 110010
using namespace std;
char Ma[N<<1],s[N];
int Mp[N<<1];
void Manacher(int len)
{
    int l=0;
    Ma[l++]='$';
    Ma[l++]='#';
    for(int i=0; ii?min(Mp[2*id-i],mx-i):1;
        while(Ma[i+Mp[i]]==Ma[i-Mp[i]])Mp[i]++;
        if(i+Mp[i]>mx)
        {
            mx=i+Mp[i];
            id=i;
        }
    }
}
int main()
{
    while(~scanf("%s",s))
    {
        int len=strlen(s);
        Manacher(len);
        int ans=0;
        for(int i=0; i<2*len+2; i++)
            ans=max(ans,Mp[i]-1);
        printf("%d\n",ans);
    }
    return 0;
}

 

你可能感兴趣的:(字符串处理题)