最长回文 HDU - 3068（马拉车算法）

Problem Description

给出一个只由小写英文字符$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

先看代码：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
#define maxx 101000
char s[maxx],s1[maxx*2];
int mp[maxx*2];
void Manacher(int l)
{
    int k=0;
    s1[k++]='$';s1[k++]='#';
    for(int i=0;i<l;i++)
    {
        s1[k++]=s[i];
        s1[k++]='#';
    }
    s1[k]=0;
    int mx=0,id=0;
    for(int i=0;i<k;i++)
    {
        mp[i]=mx>i?min(mp[2*id-i],mx-i):1;
        while(s1[i+mp[i]]==s1[i-mp[i]])
            mp[i]++;
         if(i+mp[i]>mx)
         {
             mx=i+mp[i];
             id=i;
         }
    }
}
int main()
{
    while(~scanf("%s",s))
    {
        int l=strlen(s);
       int ans=0;
        Manacher(l);
        for(int i=0;i<2*l+2;i++)
            ans=max(ans,mp[i]-1);
        cout<<ans<<endl;
    }
    return 0;
}

思路：

代码虽不长，主要需要理解的就是:
mp[i]=mx>i?min(mp[2*id-i],mx-i):1;

一般想法都是从1开始慢慢向两边移动来试探，而马拉车主要就是优化了每次试探P[i]的时候不一定需要从1开始，P[i]代表当前下标i为中心的字符串的回文串半径。
下图中 j 点代表是 i 关于 id 的对称点 ,j=2*id-i;
　　　　　mx的对称点 　　　 j 　　　　　　d 　　　　　　　　i 　　　　　　mx

id表示的就是最长回文串的中心，从图观察i与j关于Id对称，i是从2开始枚举过来已经经过了j的位置，那么j位置的最长回文串就可以确定了，如图所示；如果回文串完全被id的回文串所包围，那么根据对称原理i点的回文串的长度最少就是j点回文串的长度。即如果回文串的子串也是回文串，那么这个子串关于主串中心对称而得的子串也是一个回文串。接下来要确定的就是通过j点所能确定的i点回文串的长度最多是多少。首先应该明确，如果i点跑到mx（id点回文串所确定的范围边界）外面去了，那么j点无论如何缩减范围都不可能是id回文串的子串，就不满足上面加粗的结论了。就一定只能从1开始慢慢试探。这就是当mx < i的时候，MP[i] = 1的原因了。

还有两种情况

一种就是上图中，j所确定的回文串完全被包含，即整个串都是其子串。那么i的可确定回文串范围就是j的回文串范围，MP[i]就变成了MP[j]。

还有一种情况就是j的回文串已经超出了mx的范围
　　　mx的对称点　　　　　 j　　　　　　id　　　　　　　　 i 　　　　　　mx

对于绿线以外的区域完全未知，所以必须将MP[j]减去红线外的范围才是i的可确定范围。或者理解为只有两端都去掉外面的部分之后，剩下的才是id回文串的子串，才可以对称过去成为i的回文串。然后再在已确定的范围基础上向两边扩展。

实践是检验真理的唯一标准