Hdu 3068 Manacher算法求最长回文串长度

最长回文

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

[Manacher算法理解] (http://www.felix021.com/blog/read.php?2040)

想法：只要理解了Manacher算法就很简单了。

#include <cstdio>
#include <cstring>
#include <vector>
#include <algorithm>
#include <map>
#include <iostream>
using namespace std;
char str[2*110000+10];
int p[10+2*110000];
//因为添加了一倍的特殊字符，所以要乘2.
char s[10+110000];
//原字符串。
int len;
void manacher(int len)
{
    str[0]='$';str[1]='#';
    //str[0]是为了不让数组越界，至于‘$’和‘#’可以由不与原字符串重复的任意特殊字符代替。
    for(int i=0;i<len;i++)
        str[2*(i+1)]=s[i],str[2*(i+1)+1]='#';
    //为了排除奇长度与偶长度字符串的区别。
    str[2*len+2]='\0';
    len=2*len+2;
    int id=0,mx=0;
    p[0]=0;
    for(int i=1;i<len;i++)
    {
        if(i>=p[id]+id)
     //超过了以id为中心的最长回文串的右边界。
            p[i]=1;
        else
            p[i]=min(p[2*id-i],mx-i);
     //如果没有，请看链接中的图文解释（我实在想不出比他更好的了，就不献丑了）。
        while(str[i-p[i]]==str[i+p[i]])
            p[i]++;
        if(p[i]+i>mx)
      //更新最长回文串长度。
            mx=p[i]+i,id=i;
    }
}
int main()
{
    while (~scanf("%s",s)) {
        len = strlen(s);
        manacher(len);
        int ans = 0;
        for (int i = 0; i<2*len+2; i++) {
            ans = max(ans,p[i]);
        }
        printf("%d\n",ans-1);
    //减一很重要，可以自己手动举个例子试试。