最长回文 HDU - 3068 (马拉车算法模板)

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

马拉车算法是查找一个字符串的最长回文子串的线性算法

回文串有奇数回文串和偶数回文串，为了避免这个差异，构造样例
然后从字符串开头到结尾，先找到当前右边界最大的回文子串，记对称中心为id，右界为mx， 将以i为中心的回文子串的半径记为p[i]，然后再以i为对称中心，在半径为p[i]的基础上向两边遍历，更新以i为中心的回文子串半径，如果对称中心加上回文子串半径（在while循环中是不断增大的）大于回文子串的右界时，更新mx和id。

代码：

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
char a[110010];
char Ma[220020];  //样例构造数组要开到字符串长度的2倍多
int p[220020];
void manacher(char s[],int l1) 
{
    int l=0,i;
   
    Ma[l++]='$';           //样例构造为奇数串
    Ma[l++]='#';
    for(i=0;i<l1;i++)
    {
        Ma[l++]=s[i];
        Ma[l++]='#';
    }
    Ma[l]=0;   // 清空功能


    int mx=0,id=0;   //mx是以id为中心的最长回文子串的右界 id+p[id]=mx
    for(i=1;i<l;i++)
    {
       /*核心代码*/
       if(mx>i)
        p[i]=min(p[2*id-i],mx-i);  //p[i]表示以i为对称中心能得到的最长回文半径
       else
        p[i]=1;
        
       while(Ma[i+p[i]]==Ma[i-p[i]])  //以i为对称中心，在p[i]的基础上向两边遍历，更新以i为中心的回文子串半径 最终求得最长回文子串半径
        p[i]++;
        
        /*向两边遍历，不断更新id和mx的值*/
       if(mx<i+p[i])   
       {
           mx=p[i]+i;
           id=i;
       }
    }
}
int main()
{
    while(~scanf("%s",a))
    {
        int l1=strlen(a);
        int s=0,i;
        manacher(a,l1);
        for(i=1;i<2*l1+2;i++)
            s=max(s,p[i]-1);   /*长度为原字符串中最长回文串的长度-1*/
        printf("%d\n",s);
    }
    return 0;
}