manacher模板（最长回文子串）

版权声明：版权声明：本文为博主原创文章，未经博主允许不得转载，欢迎添加友链。 https://blog.csdn.net/zzk_233/article/details/82871956

设p[i]为当前节点i最长回文串的右节点与i的距离，i为回文中心，而有些回文串是以空格为回文中心的，比如bb，那么为了使

每个回文串都是以字符为中心，把每个字符的前后都加上一个#（操作上就是后面加上#，最后在最前面加#）。

设maxp为最长的长度，p0为maxp的回文中心。那么更新p[i]的时候有两中情况。如果i<maxp时，分两种

一种是可以直接到达maxp的位置也就是p0+p[p0]-i的长度，还有一种找到i关于p0对称的位置2*p0-i，

因为那个位置的长度和i关于p0对称，所以p[2*p0-i]也是p[i]的解，但是这种情况在大于maxp时有不确定性，

所以取最小值。但为了保证答案的准确性和考虑在i>maxp时的情况，暴力验证左右是否满足回文。并时时更新maxp和p0。

#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
#define N 51000005
using namespace std;
int p[N],l,maxp,p0,ans;
char ss[N*2],a[N];
void manacher()
{
	for(int i=1;i<l;i++)
	{
		if(i<maxp)
		{
			p[i]=min(p[2*p0-i],p[p0]+p0-i);
		}else
		{
			p[i]=1;
		}
		while(ss[i-p[i]]==ss[i+p[i]])
		{
			p[i]++;
		}
		if(i+p[i]>maxp)
		{
			maxp=i+p[i];p0=i;
		}
	}
}
int main()
{
	scanf("%s",a+1);
	l=strlen(a+1);
	for(int i=1;i<=l;i++)
	{
		ss[i*2]=a[i];
		ss[i*2+1]='#';
	}
	l=l*2+2; 
	ss[0]='^',ss[1]='#',ss[l]='0';
	manacher();
	for(int i=1;i<=l;i++)
	{
		ans=max(ans,p[i]-1);
	}
	printf("%d\n",ans);
	return 0;
}