Manacher 学习笔记

引入

求解最长回文子串问题：给出一个字符串 $s$ ，求出 $s$​ 最长的回文子串长度
我们的暴力做法就是依次以每个字符或者两个字符之间的空隙为中心，当两侧对应字符相同时向两侧进行扩展，扩展到极限展时，我们就得到了该中心对应的回文子串最大长度，与 $ans$ 取 $\max$ 即可
核心代码如下（默认字符串下标从 $1$ 开始）：

#define rep(i,n) for(int i = 1;i <= n;i++)
inline int solve(char * s){
	int n = strlen(s + 1),ans = 0;
	rep(i,n){
		int l = i,r = i;
		while(l > 1 &&r < n && s[l-1] == s[r+1]) ++l,--r;
		ans = max(ans,r - l + 1);
		if(i < n && s[i] == s[i+1]){
			l = i,r = i + 1;
			while(l > 1 &&r < n && s[l-1] == s[r+1]) ++l,--r;
			ans = max(ans,r - l + 1);
		}
	}
	return ans;
}

但是这个暴力做法的复杂度显然是 $O(n^2)$ 的，太慢了，有没有更优一点的做法呢？

Manacher

首先，我们发现我们的暴力程序在判断回文子串的时候需要分长度奇偶，这个分讨有点麻烦，我们不妨在原串的每两个字符之间（包括字符首尾）加入分隔符（通常是不在字符串中出现的 $\#$）构造出新串 $t$ 。这样一来，我们就只用讨论长度为奇数的回文子串求解方法了。
我们使用一个数组 $p$ 来记录以每个字符为中心的最长回文半径（即用 $p_i$ 记录以 $t_i$ 为中心的最长回文子串的半径， $p_i$ 最小为 $1$ ）。接下来证明 $p_i-1$ 就是以 $t_i$ 为中心的最长回文子串在原串中的长度

1. 显然 $len=p_i\times 2-1$ 是 $t$ 中以 $t_i$ 为中心的最长回文子串长度。
2. 以 $t_i$ 为中心的最长回文子串一定是以 $\#$ 开头和结尾的。所以 $len$ 减去最前或者最后的 $\#$ 字符之后，就是原串中对应回文子串长度的两倍。即原串长度为 $(len-1)/2=p_i-1$ ，得证。

所以接下来就是如何求 $p$ 数组了。假设我们已经求出了 $p_1,p_2...p_{i-1}$ ，现在要求 $p_i$
我们用一个变量 $maxid$ 记录在求解 $p_i$ 之前，我们找到的回文串中，延伸到最右端的回文串的右端点位置，再用一个变量 $id$ 记录这个回文串的中心的下标。这两个变量的初值都是零。
接下来有如下代码，实现原理会稍后解释

#define rep(i,n) for(int i = 1;i <= n;i++)
void manacher(){
	……
	rep(i,m){
		p[i] = 1;
		if(maxid > i) p[i] = min(p[(id<<1)-i],maxid - i);
		while(t[i-p[i]] == t[i+p[i]]) p[i]++;
		if(i + p[i] > maxid) id = i,maxid = i + p[i];
		……
	}
	……
}

$j=id\times 2$ 是 $i$ 关于 $id$ 的对称点，根据回文串的对称性，$p_j$ 代表的回文串可以对称到 $i$ 这边，但是如果 $p_j$ 代表的回文串对称过来以后超过 $maxid$ 的话，超出的部分就不能对称过来了，所以 $p_i$ 的下限就取 $p_{id\times 2-i}$ 和 $maxid-i$ 的较小者。随后我们再暴力扩展 $p_i$ ，并更新 $id$ 和 $maxid$

以下是我的 Manacher 模板 代码

#include <bits/stdc++.h>
using namespace std;
#define rep(i,n) for(int i = 1;i <= n;i++)
#define itn int
#define ll long long
const int N = 2.3e7;
int p[N],n,m,id,maxid,ans;
char s[N],t[N];
void manacher(){
	n = strlen(s + 1),m = n << 1 | 1;
	t[0] = '*',t[m+1] = '^',t[1] = '#';//注意边界情况
	rep(i,n) t[i<<1] = s[i],t[i<<1|1] = '#';
	rep(i,m){
		p[i] = 1;
		if(maxid > i) p[i] = min(p[(id<<1)-i],maxid - i);
		while(t[i-p[i]] == t[i+p[i]]) p[i]++;
		if(i + p[i] > maxid) id = i,maxid = i + p[i];
		ans = max(ans,p[i]);
	}
	printf("%d",ans - 1);
}
int main(){
	scanf("%s",s + 1);
	manacher();
	return 0;
}

复杂度

我们的字符串 $t$ 只被遍历了一次，复杂度为 $O(n)$ （因为 $m=n\times 2+1$ ）
现在考虑我们暴力扩展 $p_i$ 的次数
需要暴力扩展 $p_i$ ，就意味着在暴力扩展前以 $t_i$ 为回文中心的回文串右端点就已经到达了 $maxid$ 。注意到我们每暴力扩展一次就会使得 $maxid$ 增加 $1$ ，而 $maxid$ 最多为 $m$ ，因此在整个过程中最多只会暴力扩展 $O(n)$ （ $O(m)$ ） 次