Manacher（马拉车）算法详解

博客摘自：Manager算法总结

1.Manacher算法原理与实现

下面介绍Manacher算法的原理与步骤。
首先，Manacher算法提供了一种巧妙地办法，将长度为奇数的回文串和长度为偶数的回文串一起考虑，具体做法是，在原字符串的每个相邻两个字符中间插入一个分隔符，同时在首尾也要添加一个分隔符，分隔符的要求是不在原串中出现，一般情况下可以用#号。下面举一个例子：

（1）Len数组简介与性质

Manacher算法用一个辅助数组Len[i]表示以字符T[i]为中心的最长回文字串的最右字符到T[i]的长度，比如以T[i]为中心的最长回文字串是T[l,r],那么Len[i]=r-i+1。
对于上面的例子，可以得出Len[i]数组为:

Len数组有一个性质，那就是Len[i]-1就是该回文子串在原字符串S中的长度，至于证明，首先在转换得到的字符串T中，所有的回文字串的长度都为奇数，那么对于以T[i]为中心的最长回文字串，其长度就为2*Len[i]-1,经过观察可知，T中所有的回文子串，其中分隔符的数量一定比其他字符的数量多1，也就是有Len[i]个分隔符，剩下Len[i]-1个字符来自原字符串，所以该回文串在原字符串中的长度就为Len[i]-1。
有了这个性质，那么原问题就转化为求所有的Len[i]。下面介绍如何在线性时间复杂度内求出所有的Len。

（2）Len数组的计算（这个才是重点~前面的并没啥子用哟）

首先从左往右依次计算Len[i]，当计算Len[i]时，Lenj已经计算完毕。设P为之前计算中最长回文子串的右端点的最大值，并且设取得这个最大值的位置为po，分两种情况：
第一种情况：i<=P
那么找到i相对于po的对称位置，设为j，那么如果Len[j] < P - i，如下图：

由对称性可知 i 和 j 的周围的字符是一样的，所以Len[i]=Len[j]。
如果Len[j]>=P-i,由对称性，说明以i为中心的回文串可能会延伸到P之外，而大于P的部分我们还没有进行匹配，所以要从P+1位置开始一个一个进行匹配，直到发生失配，从而更新P和对应的po以及Len[i]。

第二种情况: i>P
如果i比P还要大，说明对于中点为i的回文串还一点都没有匹配，这个时候，就只能老老实实地一个一个匹配了，匹配完成后要更新P的位置和对应的po以及Len[i]。

#include<bits/stdc++.h>
using namespace std;
const int MAX_N = 1e6+9;
char str[MAX_N];
char s[MAX_N<<1];
int len[MAX_N<<1];
int val[27];
int presum[MAX_N];

int init(char *str)
{
    int len = strlen(str);
    s[0] = '@';//设置一个特殊字符
    for(int i=1;i<=2*len;i+=2)
    {
        s[i] = '#';
        s[i+1] = str[i/2];
    }
    s[2*len+1] = '#';
    s[2*len+2] = '~';//这里一定要和第一个设置的特殊字符不一样不然会出错
    s[2*len+3] = 0 ;
    return 2*len+1;//返回转换字符串的长度
}
void manacher(char * s,int ll)
{
    int mx = 0 , po = 0;
    for(int i=1;i<=2*ll;i++)
    {
        if(mx > i) len[i] = min(mx-i,len[2*po-i]);
        else len[i] = 1;
        while(s[i-len[i]] == s[i+len[i]] ) len[i] ++ ;
        if(len[i]+i > mx)
        {
            mx = len[i]+i;
            po = i;
        }
    }
    return ;
}