HDOJ3068最长回文 ---- Manacher算法

题目链接 http://acm.hdu.edu.cn/showproblem.php?pid=3068

 

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

题目是一道马拉车的模板题，直接用manacher模板就可以在O（n）时间内找出最长回文子串的长度。

详情参考 https://segmentfault.com/a/1190000008484167

 

 

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#define MAXN 110005
using namespace std;
char old_s[MAXN];// 输入的字符串
char new_s[2*MAXN];// 转换的新字符串
int p[2*MAXN];
// 初始化
void Init() {
    memset(p,0,sizeof(p));
    new_s[0] = '$';// 避免越界
    new_s[1] = '#';
    int len = strlen(old_s);
    int i = 0;
    int j = 2;
    while(i < len) {
        new_s[j++] = old_s[i++];
        new_s[j++] = '#';
    }
    new_s[j++] = '!';// 避免越界
    new_s[j] = 0;// 字符串分隔符
}

// 马拉车算法
int Manacher() {
    int ret = 0;// 返回值
    int mx = 2;
    int id = 1;
    p[1] = 1;
    int len = strlen(new_s);
    for (int i = 2; i < len - 1; i++) {
        if(i < mx) {
            p[i] = min(p[2*id-i],mx-i);
        } else {
            p[i] = 1;
        }
        while( new_s[i+p[i]] == new_s[i-p[i]])
            p[i]++;
        if(i+p[i] > mx) {
            // 更新
            mx = i + p[i];
            id = i;
        }
        ret = max(ret,p[i]-1);
    }
    return ret;
}

int main()
{
    while(scanf("%s",old_s) != EOF) {
        Init();     
        printf("%d\n",Manacher());
    }
    
    return 0;
}