Manacher‘s Algorithm 马拉车算法——查找一个字符串的最长回文子串的线性方法

Manacher’s Algorithm 马拉车算法

这个博主写的非常清楚了
给定字符串 s ,将处理后的字符串放入数组 t.
处理前：ababa 处理后：@#a#b#a#b#a#

始终注意各个元素的含义，数组 P[i] 用来存储以字符 t[i] 为中心的回文子串的半径，若p[i] = 1，则该回文子串就是t[i]本身。

辅助变量mx和id，其中id为能延伸到最右端的位置的那个回文子串的中心点位置，mx是回文串能延伸到的最右端的位置

最长子串的长度是半径减1，在原字符串的起始位置是中间位置减去半径再除以2。

重点解释一下这部分代码：

int mx = 0, id = 0, resLen = 0, resCenter = 0;
    for (int i = 1; i < t.size(); ++i) {
        p[i] = mx > i ? min(p[2 * id - i], mx - i) : 1;
        while (t[i + p[i]] == t[i - p[i]]) ++p[i];
        if (mx < i + p[i]) {
            mx = i + p[i];
            id = i;
        }
        if (resLen < p[i]) {
            resLen = p[i];
            resCenter = i;
        }
    }

mx 是当下找到的回文串中，右边界在最右的位置的右边（mx不属于回文串内部）。id就是这个回文串的中心位置。

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

这是一段对于p[i]进行初始化的过程。
对于当下讨论的第i个字符（以t[i]为中心的回文串）：

• 如果mx>i，也就是说t[i]这个字符是以t[id]为中心的回文串内部的一个字符，那么我们能找到与之对称的另一个字符t[j]，且这个t[j]的索引为 id-( i - id ) = 2 * id -i.

若p[j] < mx -i ，也就是说在这个半径内，i 和 j 的情况是一样的，则p[i]至少为p[j]，其右的字符需要继续遍历讨论。

若p[j] > mx - i，半径延伸到mx之外了，由于还不知道mx右边字符的情况，则p[i]至少为 mx - i ， mx右边的字符要继续遍历讨论。

• 如果 mx < i ，说明 t[i]不在t[id]为中心的回文串中，那么无法确定p[i]的情况，无法初始化，只能对其继续遍历讨论。
• 因此下面通过一个while循环来确定化p[i]（老老实实的一一匹配）。

while (t[i + p[i]] == t[i - p[i]]) ++p[i];

至此p[i]就完全确定了。

接下来，更新id、mx，以及最长回文串中心resCenter及其半径resLen。

//更新最右边界和中心点
if (mx < i + p[i]) {
            mx = i + p[i];
            id = i;
        }
        //更新最长回文串
if (resLen < p[i]) {
            resLen = p[i];
            resCenter = i;
        }

整体代码如下：

#include <vector>
#include <iostream>
#include <string>

using namespace std;

string Manacher(string s) {
    // Insert '#'
    string t = "$#";
    for (int i = 0; i < s.size(); ++i) {
        t += s[i];
        t += "#";
    }
    // Process t
    vector<int> p(t.size(), 0);
    int mx = 0, id = 0, resLen = 0, resCenter = 0;
    for (int i = 1; i < t.size(); ++i) {
        p[i] = mx > i ? min(p[2 * id - i], mx - i) : 1;
        while (t[i + p[i]] == t[i - p[i]]) ++p[i];
        if (mx < i + p[i]) {
            mx = i + p[i];
            id = i;
        }
        if (resLen < p[i]) {
            resLen = p[i];
            resCenter = i;
        }
    }
    return s.substr((resCenter - resLen) / 2, resLen - 1);
}

int main() {
    string s1 = "12212";
    cout << Manacher(s1) << endl;
    string s2 = "122122";
    cout << Manacher(s2) << endl;
    string s = "waabwswfd";
    cout << Manacher(s) << endl;
}