Manacher Algorithm

Manacher算法介绍

---

• Manacher算法是查找一个字符串的最长回文子串的线性算法。
• 在介绍算法之前，首先介绍一下什么是回文串，所谓回文串，简单来说就是正着读和反着读都是一样的字符串，比如abba，noon等等，一个字符串的最长回文子串即为这个字符串的子串中，是回文串的最长的那个。
• 计算字符串的最长回文字串最简单的算法就是枚举该字符串的每一个子串，并且判断这个子串是否为回文串，这个算法的时间复杂度为O(n^3)的，显然无法令人满意，稍微优化的一个算法是枚举回文串的中点，这里要分为两种情况，一种是回文串长度是奇数的情况，另一种是回文串长度是偶数的情况，枚举中点再判断是否是回文串，这样能把算法的时间复杂度降为O(n^2)，但是当n比较大的时候仍然无法令人满意，Manacher算法可以在线性时间复杂度内求出一个字符串的最长回文字串，达到了理论上的下界。

 

Manacher算法原理与实现

---

• Manacher算法提供了一种巧妙地办法，将长度为奇数的回文串和长度为偶数的回文串一起考虑，具体做法是在原字符串的每个相邻两个字符中间插入一个分隔符，同时在首尾也要添加一个分隔符，分隔符的要求是不在原串中出现。

• Manacher算法用一个辅助数组len[i]表示以字符s[i]为中心的最长回文字串的半径（包括中心字符）

• len数组有一个性质，那就是p[i]-1就是该回文子串在原字符串S中的长度，至于证明，首先在转换得到的字符串T中，所有的回文字串的长度都为奇数，那么对于以T[i]为中心的最长回文字串，其长度就为2*len[i]-1,经过观察可知，T中所有的回文子串，其中分隔符的数量一定比其他字符的数量多1，也就是有p[i]个分隔符，剩下p[i]-1个字符来自原字符串，所以该回文串在原字符串中的长度就为Len[i]-1。
• 有了这个性质，那么原问题就转化为求所有的Len[i]。
• 如下图

• 设置两个变量，maxn 和 id 。maxn 代表以 id 为中心的最长回文的右边界，也就是maxn=id+len[id]
• 假设我们现在求len[i]，也就是以 i 为中心的最长回文半径，如果 i<maxn，如上图，那么：

if(i<maxn)p[i]=min(p[2*id-i],maxn-i);

• 2*id-1为 i 关于 id 的对称点 j ,因此我们可以利用len[j]来加快查找

 

Code

---

#include<bits/stdc++.h>
#define rep(i,a,b) for(int i=(a);i<=(b);i++)
using namespace std;
const int N=1e8;
char a[N],s[N];
int maxn,len,id,ans,p[N];
int init(){
	len=strlen(a+1);
	int tmp=2;
	rep(i,1,len)s[++tmp]=a[i],s[++tmp]='#';
	s[1]='@',s[2]='#',s[++tmp]='&';
	return tmp;
}
void Manacher(){
	len=init();maxn=0;
	rep(i,1,len){
		if(i<maxn)p[i]=min(p[2*id-i],maxn-i);
		else p[i]=1;
		while(s[i+p[i]]==s[i-p[i]])p[i]++;
		if(maxn<i+p[i])id=i,maxn=i+p[i];
		ans=max(ans,p[i]-1);
	}
	printf("%d",ans);
}
int main()
{
	scanf("%s",a+1);
	Manacher();
	return 0;
}

练习题目

---

• Manacher模板题目

 

鸣谢

---

• https://subetter.com/articles/manacher-algorithm.html
• https://www.cnblogs.com/z360/p/6375514.html