（ 字符串专题 ）【 Manacher算法 】

（ 字符串专题 ）【 Manacher算法 】

推荐视频：https://www.bilibili.com/video/av41883533?p=2

1975年，Manacher发明了Manacher算法（中文名：马拉车算法），是一个可以在O(n)的复杂度中返回字符串s中最长回文子串长度的算法，十分巧妙。

让我们举个栗子：

1.字符串：abbababa        最长回文子串：5（abbababa）

2.字符串：abcbbabbc      最长回文子串：7（abcbbabbc）

3.字符串：abccbaba        最长回文子串：6（abccbaba）

算法过程分析

回文分为奇回文（ababa）和偶回文（abba），这里比较难以处理，我们使用一个小技巧（很重要）。我们将字符串首尾和每个字符间插入一个字符（注意：这个自符在串中并未出现）例如：'#'.

这样原串中的偶回文（abba）与奇回文（adcacda），变成了（#a#d#d#a#）与（#a#d#c#a#c#d#a#）两个奇回文。

定义数组Mp[]，用Mp[i]表示以i为中心的最长回文半径，求出的所有Mp中的最大值减1就是原字符串的最长回文子串

 

现在的问题就是如何 O(n) 的求Mp了。

对于上图的定义：

1、 i 是我当前需要求的Mp[ i ] 的位置，这时候前 i-1 的mp我们都已经知道了。

2、 id是前 i-1 中Mp最大的位置， mx 是以id为中心的回文的右边界。

3、 j 点是 i 点关于id 点的对称点。因为此时mx的对称点到mx是关于id的回文串。所以 i 的回文串至少等于 j 的回文串长度。之后我们在现有长度上在扩展。

 

模板代码：

#include <bits/stdc++.h>

using namespace std;

char s[51000100];
char Ma[51000100];
int Mp[51000100];

int Manacher()
{
    int len = strlen(s);
    int l = 0;
    Ma[l++] = '$';
    Ma[l++] = '#';
    for ( int i=0; i<len; i++ ) {
        Ma[l++] = s[i];
        Ma[l++] = '#';
    }
    Ma[l] = 0;
    int mx=0,id=0;
    for ( int i=0; i<l; i++ ) {
        Mp[i] = mx>i?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 ans = 0;
    for ( int i=0; i<l; i++ ) ans=max(ans,Mp[i]);
    return ans-1;
}

int main()
{
    scanf("%s",s);
    cout << Manacher() << endl;

    return 0;
}

交题连接：https://www.luogu.com.cn/problem/P3805