8.31 Manacher(马拉车算法)理解+模板例题代码

---

一、概念

Manacher算法是查找一个字符串的最长回文子串的线性算法。

计算字符串的最长回文字串：
枚举该字符串的每一个子串，并且判断这个子串是否为回文串 O(n^3) - - -优化：枚举回文串的中点，该算法有两种情况，一种是回文串长度是奇数的情况，另一种是回文串长度是偶数的情况， O(n^2) - - -再优化：Manacher算法 O(n)，将长度为奇数的回文串和长度为偶数的回文串一起考虑

原理

在原字符串的每个相邻两个字符中间插入一个分隔符，同时在首尾也要添加一个分隔符，分隔符的要求是不在原串中出现，一般情况下可以用#号。在转换得到的字符串T中，所有的回文字串的长度都为奇数

Manacher算法有一个辅助数组Len[ ]，Len[i]表示以字符T[i]为中心的最长回文字串的最右字符到T[i]的长度，比如以T[i]为中心的最长回文字串是T[l,r],那么Len[i]=r-i+1

注意，Len[i]-1就是该回文子串在原字符串S中的长度（证明略）

在manacher算法中Len数组的计算：首先从左往右依次计算Len[i]，当计算Len[i]时，Lenj已经计算完毕。设mx为之前计算中最长回文子串的右端点的最大值，并且设取得这个最大值的位置为po，分两种情况：第一种情况：i<=mx ;第二种情况: i>P
(看了其他博客都有图解，说实话我看图 越看越不懂，看代码理解了 qwq)

二、【模板】manacher 算法

链接: P3805 【模板】manacher 算法

/**/
#include<bits/stdc++.h>
using namespace std;
const int maxn=11000005;
char s[maxn],t[maxn<<1];
int Len[maxn<<1];
//转换原始串
void init(char *st){
	int len=strlen(st);
	t[0]='@';//字符串开头增加一个特殊字符，防止越界
	for(int i=1;i<=2*len;i+=2){//转换后字符串是从1开始
		t[i]='#';
		t[i+1]=st[i/2];
	}
	t[2*len+1] ='#';
	t[2*len+2] ='$';//字符串结尾加一个字符，防止越界
	//return 2*len+1;
}
//Manacher算法
int manacher(char *st,int len){
	int mx=0,ans=0,po;//mx即为当前计算回文串最右边字符的最大值,po为取得这个最大值的位置中心点 
	for(int i=1;i<len;i++){
		if(mx>i) Len[i]=min(mx-i,Len[2*po-i]);//在Len[j]和mx-i中取个小;i j关于po对称 
		else Len[i]=1; //如果i>=mx，要从头开始匹配
		while(st[i-Len[i]]==st[i+Len[i]]) Len[i]++;
		if(i+Len[i]>mx){//若新计算的回文串右端点位置大于mx，要更新po和mx的值
			po=i;
			mx=i+Len[i];
		}
		ans=max(ans,Len[i]);//返回Len[i]中的最大值-1即为原串的最长回文子串额长度
	}
	return ans-1;
}
int main(){
	cin>>s;
	init(s);
	cout<<manacher(t,2*strlen(s)+1);
}