HDU3068 最长回文【manacher算法】

最长回文

Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 32452 Accepted Submission(s): 11908

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

Source
2009 Multi-University Training Contest 16 - Host by NIT

问题链接：HDU3068 最长回文
问题简述：（略）
问题分析：
    该题是计算字符串中最长回文的模板题，套manacher算法的函数计算就可以了。
    需要注意的是使用的存储空间要合适。字符串长度为n时，需要2n+3的存储空间。
程序说明：（略）
参考链接：（略）
题记：（略）

AC的C++语言程序如下：

/* HDU3068 最长回文 */

#include <iostream>
#include <stdio.h>
#include <string.h>

using namespace std;

#define N 110000
char s[N + N + 3];
int p[N + N + 3];

int manacher(char s[], int p[])
{
    int len = strlen(s), id = 0, maxlen = 0;
    for(int i = len; i >= 0; i--) {     /* 插入'#' */
        s[i + i + 2] = s[i];
        s[i + i + 1] = '#';
    }
    s[0] = '*';
    for(int i = 2; i < 2 * len + 1; i++) {
        if(p[id] + id > i)
            p[i] = min(p[2 * id - i], p[id] + id - i);
        else
            p[i] = 1;
        while(s[i - p[i]] == s[i + p[i]])
            p[i]++;
        if(id + p[id] < i + p[i])
            id = i;
        maxlen = max(maxlen, p[i]);
    }
    return maxlen - 1;
}

int main(void)
{
    while(scanf("%s", s) != EOF)
        printf("%d\n", manacher(s, p));

    return 0;
}