『Manacher算法』学习总结

$\mathrm{M}\mathrm{a}\mathrm{n}\mathrm{a}\mathrm{c}\mathrm{h}\mathrm{e}\mathrm{r}$ 算法的作用&原理

求一个序列的最大回文子串，我们需要需要用到Manacher算法。其作用在于能够在O(n)的时间内求出最长回文子串，同时也能够求出回文子串的个数。且时间效率高，代码十分简洁。

我们知道，回文子串分为奇数回文子串和偶数回文子串。由于两种情况的处理较为麻烦，我们可以考虑在期间插入字符在简化问题。

例如例子：$abababc$中，如果我们队每相邻的两个数之间插入字符#,我们可以得到字符串：$$\#a\#b\#a\#b\#a\#b\#c\#$$

这样，我们就将原问题转化为了奇数的回文子串。对于在第$i$个位置的字符$a_i$，我们可以求出最长的回文半径$p_i$.那么我们有如下结论：

• 结论1：字符串的最长回文子串为$\max (p_i-1)$.
• 结论2：字符串的回文子串个数为：$$\sum \lfloor \frac{p_i}{2}\rfloor$$

那么我们只要完成一遍namacher算法，就能够得到相应的$p$数组解决相应的问题。现在我们我们考虑如何求解$p$数组。

---

$\mathrm{M}\mathrm{a}\mathrm{n}\mathrm{a}\mathrm{c}\mathrm{h}\mathrm{e}\mathrm{r}$ 算法的实现

对于Manacher算法，我们需要维护一个扩展到的最右边的回文位置$mx$和回文中心$id$。

即如果当前位置为$i$，有位置$i+p[i]$为当前的最右位置，那么$mx=\max (i+p[i])$,$id$是对应的$i$。那么这个算法有什么用呢？如图所示：

我们根据之前对称中心$id$和当前位置$i$，得到$id$左边与i对应的位置$j$。$j$的坐标是：$$j=id-(i-id)=2\times id-i$$
我们会发现：

• 当$i+p[j]\le mx$时，i的回文半径至少是$p[j]$.根据这个大回文的性质可以得出。
• 当$i+p[j]>mx$时，i的回文半径至少是$mx-i$.其实也很显然。

那么在初始化时，我们就有：$$p[i]=\min (p[id\times 2-i],mx-i)$$

然后我们进行暴力扩展即可。

---

$\mathrm{M}\mathrm{a}\mathrm{n}\mathrm{a}\mathrm{c}\mathrm{h}\mathrm{e}\mathrm{r}$ 算法的复杂度分析

当$p[i]=p[j]$时，一次的求解复杂度为$O(1)$.

当$p[i]!=p[j]$时，一定有$p[i]+i>mx$则更新了$mx$，而$mx$的更新复杂度不超过$n$，所以时间复杂度不超过$O(n)$.

故manacher算法的总时间复杂度为$O(n)$.

---

$\mathrm{M}\mathrm{a}\mathrm{n}\mathrm{a}\mathrm{c}\mathrm{h}\mathrm{e}\mathrm{r}$ 算法代码实现

为了防止超出了字符串的范围后继续拓展，我们应该在字符串的末尾加一个不一样的字符即可，当时开头也可以加一个不一样的字符，但需要保证开头结尾不同。

代码如下：

#include <cstdio>
#include <string.h>
#include <iostream>

using namespace std;
const int N = 23000000;

int p[N];
char s[N], a[N];

void Malacher(int n)
{
	int mx = 0, id = 0;
	for (int i=1;i<=n;++i)
	{
		p[i] = i < mx ? min(p[id*2-i],mx-i) : 1;
		while (a[i-p[i]] == a[i+p[i]]) p[i] ++;
		if (i + p[i] > mx) mx = i + p[i], id = i;
	}
	return;
}

int main(void)
{
	cin>>s+1;
	int n = strlen(s+1);
	a[1] = '#';
	for (int i=1;i<=n;++i) 
	{
		a[i*2] = s[i];
		a[i*2+1] = '#';	
	}	
	a[n*2+2] = '*';
	Malacher(2*n+1);
	int ans = 0;
	for (int i=1;i<=n*2+1;++i)
	    ans = max(ans,p[i]-1);
	cout<<ans<<endl;
	return 0;
}