manacher算法总结

manacher是求最长回文子串的算法，时间复杂度仅为O(n)。在了解其原理之前先考虑求最长回文子串的朴素做法：中心扩展法。中心扩展法很简单，就是枚举每个点，暴力地向两边匹配，不能匹配时就停止，复杂度O(n^2)。而manacher就是对中心扩展法的优化版本，其核心思想在于充分利用已找到的回文串来减少中心扩展匹配次数。

粗略介绍一下大体思路：

考虑字符串abababa，从1开始标号，为举例方便只考虑奇数长度回文子串。开始时用中心扩展法暴力地跑完第3个字符'a'时，可以得到以第1个，第2个以及第3个字符为中心的最长回文串，分别是"a", "aba", "ababa"。接下来要对第4个字符'b'中心扩展，此时我们已知"ababa"是回文串且第4个字符在这个回文串中，由于回文的特性，我们可以找到第4个字符关于中心字符对称的字符，也就是第2个字符'b'，而且第2个字符为中心的回文串长度已知，那么第4个字符为中心的回文串长度至少应该等于第2个字符为中心的回文串长度（镜像对称），具体多长只需要在此基础上继续中心扩展就好了。时间复杂度O(n)，因为for循环中的while最多总共执行n次。

最后放一个模板：

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cctype>
using namespace std;

int p[22000010];//p[i]记录s[i]能单侧延伸多长 
char s[22000010];

int manacher(char *s)
{
	int r = 0;//最右边界
	int c = 0;//当前中心 
	int len = strlen(s)-1, ans = 0;
	for(int i = 1; i <= len; i++)
	{
		if(i < r) 
			p[i] = min(r-i, p[2*c-i]);//点i关于c的对称点2*c-i 
		else
			p[i] = 1;//单个字符延伸长度为1 
		while(s[i+p[i]] == s[i-p[i]])
			p[i]++;
		if(i+p[i] > r)//更新中心点及最右边界 
		{
			r = i+p[i];
			c = i;
		} 
		ans = max(ans, p[i]);
	}
	return ans-1;//当前位置的p[i]-1就是该位置能构成的最长回文串 
} 

signed main()
{
	//预处理s，处理为%#a#b#c#这样形式 
	char ch;
	int cnt = 0;
	while(isalpha(ch = getchar()))
	{
		s[++cnt] = '#';
		s[++cnt] = ch;
	}
	s[0] = '%';
	s[++cnt] = '#';
	s[++cnt] = '\0';
	cout << manacher(s) << endl;
	
    return 0;
}