这题可用拓展KMP分治法来做http://blog.sina.cn/dpool/blog/s/blog_677a3eb30100knj8.html 复杂度O(nlogn)这种方法好复杂而且代码很长,不易理解。
相比之下Manacher就简单多了,算法本身也很简单 这里个易懂的资料http://wenku.baidu.com/view/3031d2d3360cba1aa811da42.html
复杂度O(n)
kuangbin的模版(有一处修改)
/*
* 求最长回文子串
*/
const int MAXN=110010;
char Ma[MAXN*2];
int Mp[MAXN*2];
void Manacher(char s[],int len)
{
int l=0;
Ma[l++]='$';
Ma[l++]='#';
for(int i=0;i<len;i++){
Ma[l++] = s[i];
Ma[l++] = '#';
}
Ma[l]=0;
int mx=0,id=0;//mx记录在求i之前的回文串中,延伸至最右端的位置的下一格.
for(int i=1;i<l;i++){
//i=0的位置永远用不到,因为不可能有mx之内的点关于id的对称点是i=0,所以i从1开始,
//从0开始的时候会导致后面i-Mp[i]=-1,binshen的模板在这有点小问题
Mp[i] = mx>i? min(Mp[2*id-i],mx-i):1;
while(Ma[i+Mp[i] ]==Ma[i-Mp[i] ]) Mp[i]++;
if(i+Mp[i]>mx){
mx=i+Mp[i];
id=i;
}
}
}//249MS 2184K
#include<cstdio>
#include<iostream>
#include<cstring>
#define M 111111
using namespace std;
char s[M];
char Ma[M<<1];
int Mp[M<<1],ans;
void Manacher(char s[],int len)
{
int l=0;
Ma[l++]='$';
Ma[l++]='#';
for(int i=0;i<len;i++){
Ma[l++] = s[i];
Ma[l++] = '#';
}
Ma[l]=0;
int mx=0,id=0;
for(int i=1;i<l;i++){
//i=0的位置永远用不到,因为不可能有mx之内的点关于id的对称点是i=0,所以i从1开始,
//从0开始的时候会导致后面i-Mp[i]=-1,binshen的模板在这有点小问题
Mp[i] = mx>i? min(Mp[2*id-i],mx-i):1;
while(Ma[i+Mp[i] ]==Ma[i-Mp[i] ]) Mp[i]++;
if(i+Mp[i]>mx){
mx=i+Mp[i];
id=i;
}
if(Mp[i]>ans) ans=Mp[i];
}
}
int main()
{
while(~scanf("%s",s)){
ans=-(1<<30);
int len=strlen(s);
Manacher(s,len);
printf("%d\n",ans-1);
}
return 0;
}
原文:http://blog.csdn.net/kalilili/article/details/43890181