HDOJ-3068 最长回文 (manacher求最长回文串)

最长回文
Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)

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
using namespace std;
const int N = 1e6;
int RL[N];//记录每个位置的最长回文数
char origin[N], modified[N];//原始字符串和插入#后的字符串
int manacher(char *s) {
    int res = 0;
    //插入字符#
    int len = strlen(s);
    int l = 0; modified[l++] = '@'; modified[l++] = '#';
    for (int i = 0; i < len; i++) {
        modified[l++] = s[i];
        modified[l++] = '#';
    }
    //初始化能达到的右边边界R和此边界围绕的轴心pos
    int pos = 0; int R = 0;
    for (int i = 0; i < l; i++) {
        //当前点再右边边界左边
        if (i < R) RL[i] = min(RL[2 * pos - i], R - i);
        //如果在边界右边只能暴力搜索
        //先初始化为1
        else RL[i] = 1;
        //从i开始左右对比，找到以i为中心的最长回文串，如果碰到两端字符不相同或者到边界就停止
        while (modified[i - RL[i]] == modified[i + RL[i]]) RL[i]++;
        //更新pos和R
        if (i + RL[i] > R) {
            R = i + RL[i];
            pos = i;
        }
        res = max(res, RL[i] - 1);
    }
    return res;
}
int main() {
    while (scanf("%s", origin) != EOF) {
        printf("%d\n", manacher(origin));
    }
}