题目链接:点击打开链接
经典的求最长回文子串题目,需要用到manacher算法。
manacher算法的具体介绍:http://blog.csdn.net/ggggiqnypgjg/article/details/6645824
#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=1000005;
char s[maxn];
char str[maxn*2];//用来存储s更新后的字符串
int p[maxn*2]; //p[i]表示以str[i]为中心的最长回文串的半径
int len1,len2;
void init() //初始化操作,更新s到str中,将原字符串两两间隔插入'#'变成奇数串,以便于处理
{
len1=strlen(s);
str[0]='$'; //防止数组越界,字符串首端用'$'标记
str[1]='#';
for(int i=0;i<len1;i++)
{
str[i*2+2]=s[i];
str[i*2+3]='#';
}
str[2*len1+2]='\0';//同样是防止越界的操作,我个人习惯用'\0'
}
void manacher()
{
init();
memset(p,0,sizeof p);
int id=0,mx=0; //id<i,mx=i+p[id]表示在i之前回文串的最大半径距离i的最短距离
len2=strlen(str);
for(int i=1;i<len2;i++)
{
if(mx>i)
p[i]=min(p[2*id-i],mx-i); //2*id-i表示i关于id对称的点,根据对称性可得到此行的关系式
else p[i]=1;
while(str[i+p[i]]==str[i-p[i]]) p[i]++;//从i+p[i]继续向外延伸判断回文串是否能延长
if(i+p[i]>mx)
{
mx=i+p[i];
id=i;
}
}
}
int main()
{
int cnt=1;
while(~scanf("%s",s))
{
if(s[0]=='E') break;
manacher();
int ans=0;
for(int i=0;i<len2;i++)
ans=max(ans,p[i]);
printf("Case %d: %d\n",cnt++,ans-1);
}
return 0;
}