【Manacher】最长回文子串

caioj任意门

hz2016评测传送门

可以的话来一下hz2016评测吧，有的题caioj没有的我也可以给到数据嘛。

#include<map>
#include<queue>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define Maxchar 100000
#define mes(x,y) memset(x,y,sizeof(x));
#define mpy(x,y) memcpy(x,y,sizeof(x))
#define INF 2147483647
using namespace std;
int p[2*Maxchar+1];
char s[2*Maxchar+1],now[2*Maxchar+1];
int Manacher(){
int len=strlen(s+1);
for(int i=1;i<=len;i++)now[2*i-1]='#',now[2*i]=s[i];
len=len*2+1;
now[len]='#';
int pos=0,R=0;//R表示当前访问到的所有回文子串，所能触及的最右一个字符的位置,pos表示R位置的回文串的中心
for(int i=1;i<=len;i++){
int j=2*pos-i;//j为i关于pos的对称点
if(i<=R)p[i]=min(p[j],R-i);
/*
如果i<=R的话，分两种情况
第一种情况p[j]>R-i时，表示j对称点的最长回文子串已经越出R的界限了
这时，因为我们不确定大于R时的情况，所以p[i]暂时等于R-i
第二种情况p[j]<=R-i时，那就可以直接继承p[j]得p[i]=p[j]
综合以上两种情况，我们可以用p[i]=min(p[j],R-i)来归纳
*/
else p[i]=1;//不然就要暴力枚举
while(1<=i-p[i]&&i+p[i]<=len&&now[i-p[i]]==now[i+p[i]])p[i]++;
//向左右枚举，更新p[i]，即是作为继承后的补充，又是作为暴力枚举的操作
if(i+p[i]-1>R){pos=i;R=i+p[i]-1;}//适时更新pos和R
}
int ans=0;
for(int i=1;i<=len;i++)ans=max(ans,p[i]-1);
/*
首先当前的最长回文子串长度为2*p[i]-1
因为我们得到的p数组是在加了#号后的字符串上操作的，所以我们要对答案进行处理
因为#号处于首尾和每个字符之间，所以我们就可以保证所得出的最长回文子串的首尾都为#
这时我们可以得出不带#号的回文串的长度为(2*p[i]-1-1)/2=p[i]-1
所以真正的最长回文子串就是p[i]-1
ans记录最长的回文子串长度
*/
return ans;
}
int main(){
while(scanf("%s",s+1)!=EOF){
printf("%d\n",Manacher());
}
return 0;
}