Manacher算法-求字符串中最长回文串

一、算法原理

Manacher算法在对求字符串中最长回文串问题中，具有O(n)时间和空间复杂度。算法的精妙之处在于巧妙的利用了回文串的对偶性质。

第一步：对字符间添加间隔符，（例如，字符串aababcdab，通过添加间隔符转化为 #a#a#b#a#b#c#d#a#b#）从而避免了字符串的奇偶性问题，为了避免算法中考虑边界问题，可以在字符串首尾加入奇异符号（例如，$#a#a#b#a#b#c#d#a#b#&，因为C/C++字符串中有结束符\0，所以可可以省略尾端加入的奇异符）；

第二步：更新模式数组p[]，这里也是Manacher算法的核心所在，我当初看这块时纠结了好长时间。首先通过图片来说明一下：
假设当前处理位置在i点，且 i < mx（mx为关于坐标k对称子串的右边界为），即P[k] = mx - k。此时忽略左边界mx’又 i 关于 k 对称的坐标为 j ，这时就可以分成两种情况：

1. P[j] > mx - i + 1 时，P[i] = mx - i + 1；因为 mx 和 mx’ ，i 和 j 为关于k对称，所以 S[ mx’, ~ ,j ] = S[ i, ~ ,mx ]。如果 P[i] > mx - i + 1，也即 S[mx+1] = S[2*i - mx - 1]，而已知：S[mx’ -1] = S[2*j - mx’ + 1 ]，j + i= 2*k，mx’ + mx= 2*k ，所以有S[mx’-1] = S[mx+1]，这时以k为中心对称子串右边界变为了mx+1，矛盾。
   

2. P[j] == mx - i + 1 时，P[i] >= P[j]；证明方法和上面一样。

3. P[j] <”mx - i + 1 时，P[i] = P[j]；证明方法和上面一样。
   对于其他情况就需要重新匹配了。
   

第三步：找到（也可以在更新P[]时记录最大数）模式数组P[]中最大数减一，即为字符串中最大回文串长度。

二、程序

#include <iostream>#include <string>using namespace std;int LongestPalindromicSubstring(const string& S, int& maxPos){    int *p = new int[S.size() + 1];    memset(p, 0, sizeof(p));    int mx = 0, k = 0;    int maxLength = 0;    bool reCompute = false;    for(int i = 1; i < S.size(); i++)    {        if( mx > i ) {            p[i] = (p[2*k - i] < (mx - i) ? p[2*k - i] : (mx - i + 1));            reCompute = p[2*k - i] == (mx - i + 1) ? true : false;        }        else {            p[i] = 1;            reCompute = true;        }        if( reCompute )        {            for(; S[ i+p[i] ] == S[i - p[i] ]; ++p[i]);         }        if( p[i] > maxLength)        {            maxPos = i;            maxLength = p[i];        }        if(i + p[i] > mx)        {            mx = i + p[i] - 1;            k = i;        }        cout << S[i] << "  ";        reCompute = false;    }    cout << endl;    for(int i = 1; i < S.size(); i++)    {        cout << p[i] << "  ";    }    cout << endl;    maxPos = maxPos/2 - 1 - --maxLength/2;    delete  p;    return maxLength;}int main(){    string S = "aabababab";    string S_Process;    S_Process += "$#";    for(int i = 0; i < S.size(); i++)    {        S_Process += S[i];        S_Process += "#";    }    cout << S_Process << endl;    int pos = 0;    int length = LongestPalindromicSubstring(S_Process, pos);    cout << "max Longest Palindromic Substring is : " << S.substr(pos, length) << endl;    cout << "max length : " << length << endl;     return 0;}

结果：