马拉车算法

标题

求最长回文串长度
给出一个只由小写英文字符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<stdio.h>
#include<string.h>
#include<algorithm>
#define d 110005
using namespace std;
char a[d*2];//加#号之后的字符串
int b[d*2];//回文半径
char s[d];//原字符串
void manacher(char s[],int le)
{
    int l=0;
    a[l++]='$';//防止while 越界
    a[l++]='#';//解决奇偶问题  把所有的回文串都变成奇数。
    for(int i=0;i<le;i++)
    {
        a[l++]=s[i];
        a[l++]='#';//每两个字符间加一个#
    }
    a[l]=0;
    int mx=0;//上一次操作的最长回文子串的最右端
    int id=0;//上一次操作的位置，也就是i的位置；
    for(int i=0;i<l;i++)
    {
        //my    j     id      i   mx    五个点位置
        if(i<mx)
            b[i]=min(b[2*id-i],mx-i);//也就是id 到i  或者 i到mx   二选一
        //j=2*id-i  是i关于id的对称点  两者的b[]相等。
        //如果b[j]<mx-i;说明j的长度在id的范围内；
        //my是mx关于id对称点
        //如果b[j]>mx-i;说明j的长度超出了my
        else
            b[i]=1;//在i之前没有回文串出现，初始为1；
        while(a[i+b[i]]==a[i-b[i]])//以i为中心扩展长度
            b[i]++;
        if(i+b[i]>mx)//更新
        {
            mx=i+b[i];//向右扩大mx；
            id=i;//记录上次的中心。
        }
    }
}
int main()
{
    while(scanf("%s",s)!=EOF)
    {
        int ll=strlen(s);
        manacher(s,ll);
        int ss=0;
        for(int i=0;i<2*ll+2;i++)//便利取最长的那个回文半径。
            ss=max(ss,b[i]-1);//加井号的比原来的  回文半径大一个单位。
        printf("%d\n",ss);
    }
    return 0;
}