求字符串的最长回文子串（manacher算法O(n)时间）

题目

描述

小Hi和小Ho是一对好朋友，出生在信息化社会的他们对编程产生了莫大的兴趣，他们约定好互相帮助，在编程的学习道路上一同前进。

这一天，他们遇到了一连串的字符串，于是小Hi就向小Ho提出了那个经典的问题：“小Ho，你能不能分别在这些字符串中找到它们每一个的最长回文子串呢？”

小Ho奇怪的问道：“什么叫做最长回文子串呢？”

小Hi回答道：“一个字符串中连续的一段就是这个字符串的子串，而回文串指的是12421这种从前往后读和从后往前读一模一样的字符串，所以最长回文子串的意思就是这个字符串中最长的身为回文串的子串啦~”

小Ho道：“原来如此！那么我该怎么得到这些字符串呢？我又应该怎么告诉你我所计算出的最长回文子串呢？

小Hi笑着说道：“这个很容易啦，你只需要写一个程序，先从标准输入读取一个整数N（N<=30)，代表我给你的字符串的个数，然后接下来的就是我要给你的那N个字符串（字符串长度<=10^6)啦。而你要告诉我你的答案的话，只要将你计算出的最长回文子串的长度按照我给你的顺序依次输出到标准输出就可以了！你看这就是一个例子。”

样例输入

3

abababa

aaaabaa

acacdas

样例输出

7

5

3

题解

manacher算法详解 <<

把任意字符串的长度变成奇数串的方法：abababa => #a#b#a#b#a#b#a#

#include <iostream>
#include <cstring>
using namespace std;
#define ms(a,b) memset(a,b,sizeof(a))

const int MAXN=1e6+100;
char str[MAXN*2];
int p[MAXN*2];
void Change(char *s,int len)
{//转换字符串函数
int i;
for(i=len;i>=0;i--){
s[i*2+1]='#';
s[i*2+2]=s[i];
}
s[0]='$';
}
int Manacher(char* s,int len)
{//利用Manacher算法求最长回文子串,最大p[i];
p[0]=1;
int mx=0,id=0;
for(int i=1;s[i]!='\0';i++){
p[i]=mx>i?min(p[2*id-i],mx-i):1;
while(s[i + p[i]] == s[i - p[i]]) p[i]++;
if (i + p[i] > mx) {
mx = i + p[i];
id = i;
}
}
int ans=0;
for(int i=0;i<len;i++){
ans=max(ans,p[i]-1);
}
return ans;
}
int main()
{
int T;
scanf("%d",&T);
while(T--){
ms(p,0);
scanf("%s",str);
int len=strlen(str);
Change(str,len);
len=strlen(str);
printf("%d\n",Manacher(str,len));
}
return 0;
}