manacher算法（仅用于回文）——算法详解

算法简介：

$manacher$算法是一个能够在线性时间内得到一个字符串中最长回文子串的算法，性能比其他的算法，如字符串哈希还要优秀，且不会被造数据卡住。是处理回文串最有效的工具。

算法分析：

1.首先有一个问题，那就是我们如何处理奇数回文串和偶数回文串这两种不同情况。我们可以采取如下处理方式：

• 假设存在一个字符串$str$:$abbabba$
• 将其转化为一个新的字符串$s$：$$@a@b@b@a@b@b@a@\lambda$
• 原则是首尾两个字符不同，和中间的字符也不同（防止出现访问越界），中间所有字符相同。

2.现在我们来观察一下

• 以我们的s[8]='a'，为中心的回文串向右最大拓展长度$+1$即为$len[8]=8$，而回到原字符串，我们知道以str[3]为中心的最大回文字串长度为$7=len[8]-1$。
• 以我们的s[5]='#'，为中心的回文串向右最大拓展长度$+1$即为$len[5]=5$，而回到原字符串，我们知道以str[1,2]为中心的最大回文字串长度为$4=len[5]-1$。
  由此可以得到结论，$ans=max(len[i]-1)$。

3.下面正式开始$manacher$算法，上核心代码，假设此时我们已经得到前$i-1$个$len$，那么现在来求第$i$个$len$。

len[0]=0;
ll mx=0,ans=0,id=0;//mx,id：前面的所有结果中向右扩展得最多的len和中心
for(ll i=1;i<count1;i++)
{
	if(i<mx)//如果i在最大被拓展的内部，则有两种情况
		len[i]=min(mx-i,len[2*id-i]);
	else//如果不在最大被拓展的内部，则暴力拓展匹配
		len[i]=1;
	while(s[i-len[i]]==s[i+len[i]])//这个循环解释了为什么len表示最右边+1
		len[i]++;
	if(len[i]+i>mx)
	{
		mx=len[i]+i;
		id=i;
		ans=max(ans,len[i]-1);
	}
}

4.我们只分析$i<mx$的情况，因为$i>=mx$时很简单，很容易理解。

• 我们可以找到$i$关于$id$的对称点$i^{\prime }=2\times id-i$，而$len[i^{\prime }]$我们是知道的（即红色部分），根据对称性，理想的回文即为绿色部分。

• 如果$i+len[i^{\prime }]>mx$，那么$len[i]=mx-i$，因为，超过的部分（紫色+绿色和红色+蓝色）根据关于$id$的对称性可知道是以id为中心的回文串没有拓展到的位置因为$mx$就是最大拓展，那就说明不对称，而在$mx$以内是满足对称的，所以$len[i]=mx-i$。
  

• 而当$i+len[i^{\prime }]<=mx$时，显然，这时候以$id$为中心，$mx$为半径的字串是对称的，$i+len[i^{\prime }]$被包含在这个区间内，而$i+len[i]$和$mx$中间部分无法计入答案的原因是，根据对称性，$i^{\prime }$无法匹配到，那么$i$也无法匹配到。所以$len[i]=len[i^{\prime }]$。

5.所以综上所述，在$i<mx$时，$len[i]=min(mx-i,len[2\times id-i])$。

模板程序：

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>

using namespace std;
typedef long long ll;
const ll maxn=11000010;
char fr[maxn],s[maxn];
ll len[maxn];//len[i]表示以i为中心能扩展的最大回文串的右端+1
int main()
{
	scanf("%s",fr);
	ll frlen=strlen(fr),count1=0;
	s[count1++]='*';
	for(ll i=0;i<frlen;i++)
	{
		s[count1++]='#';
		s[count1++]=fr[i];
	}
	s[count1++]='#';
	s[count1++]='!';
	len[0]=0;
	ll mx=0,ans=0,id=0;
	for(ll i=1;i<count1;i++)
	{
		if(i<mx)
			len[i]=min(mx-i,len[2*id-i]);
		else
			len[i]=1;
		while(s[i-len[i]]==s[i+len[i]])
			len[i]++;
		if(len[i]+i>mx)
		{
			mx=len[i]+i;
			id=i;
			ans=max(ans,len[i]-1);
		}
	}
	printf("%lld",ans);
	
	return 0;
}