HDU 3068 最长回文(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
  
  

 

Source

2009 Multi-University Training Contest 16 - Host by NIT

PS：

manacher算法

数组p[i]表示以i为中心的(包含i这个字符)回文串半径长

代码如下：

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define maxn 200017
int N;
int p[maxn];
char str[maxn];//原字符串
char b[maxn];//转换后的字符串
//转换原始串
void init()
{
    int i;
    for(i = 0; str[i]; i ++)
    {
        b[2 * i + 1] = '#', b[2 * i + 2] = str[i];
    }
    N = 2 * i + 1;
    b[0] = '$';
    b[N] = b[N + 1] = '#';
}
int manacher()
{
    int id;
    int maxx = 0;//maxx即为当前计算回文串最右边字符的最大值
    int ans = 0;
    for(int i = 1; i <= N; i ++)
    {
        if(maxx > i)
            p[i] = min(maxx-i, p[2*id-i]);//在Len[j]和mx-i中取个小
        else
            p[i] = 1;//如果i>=mx，要从头开始匹配
        while(b[i + p[i]] == b[i - p[i]])
            ++ p[i];
        if(i + p[i] > maxx)//若新计算的回文串右端点位置大于maxx，要更新id和maxx的值
        {
            maxx = i + p[i];
            id = i;
        }
        if((p[i]-1)>ans)
        {
            ans = p[i]-1;//返回p[i]中的最大值-1即为原串的最长回文子串额长度
        }
    }
    return ans;
}
int main()
{
    while(~scanf("%s",str))
    {
        init();//转换原始串
        int ans = manacher();
        printf("%d\n",ans);
    }
    return 0;
}