Manacher处理字符串

Manacher算法在求回文串的时候是一种很高效的算法，这里讲解一些这个算法的原理以及简单的模板。

1. 原理讲解：

   • 首先对字符串进行预处理，在字符串的中间以及两边都加上一个不会出现的字符，然后这样就可以保证我们求的字符串都是一个长度为奇数的回文串。如下：
     原始字符串s=“absbffbsma",插入一个字符’#'后形成的新字符串sn="#a#b#s#b#f#f#b#s#m#a#". 对与这样的处理，我们可以发现这个字符串有两个性质：
     1.回文串的长度总是一个奇数的
     2.回文串的两端一定是自己所插入的那个字符串。
   • 性质利用：
     第一个性质：首先说明一下如何求新形成的字符串sn的回文串，我们创建一个数组p[inf]，对应的是每一个位置i时，由位置i形成的回文串能向一边延伸的最大长度。另外创键两个整数pos，R，他们所代表的是目前位置前面的回文串中右端点最大的位置R和该回文串中心的位置pos
     求解p[i]，在前面i-1的回文串都已经求得的情况下位置i的值一个是i关于位置pos对称位置j的值p[i],另一种是i的最右端超出了R的右侧，此时暴力求出R右侧的长度。
     
     由此我们可以看出，这个算法的时间复杂度是线性的，因此是一种非常高效的算法。
     第二个性质：第二个性质主要是用来字符串sn的答案转换成原来字符串s的。在sn中，以i为中心的回文串长度为2×p[i]-1 (p[i]最小值1),由于已经保证的两端的字符串一定是我们插入的字符串，所以其中原来的字符串和插入的字符串个数分别为p[i]-1和p[i].

2. 代码实现

下面贴上一份代码：

#include<iostream>
#include<string>
using namespace std;
const int inf = 1e4;
string s;
char sn[inf];
int p[inf];
int main()
{	
	while(cin>>s){
		if(s==".")break;
		for(int i=1;i<=s.length();i++){sn[i*2-1]='$',sn[i*2]=s[i-1];}
		sn[s.length()*2+1]='$';
		int pos=0,R=0;
		int len=s.length()*2+1;
		for(int i=1;i<=len;i++){
			if(i<R)p[i]=min(p[pos*2-i],R-i);else p[i]=1;
			while(i+p[i]<=len&&i-p[i]>=1&&sn[i+p[i]]==sn[i-p[i]])p[i]++;
			if(i+p[i]>R){
				pos=i;R=i+p[i];
			}
		}
		int ans=0;pos=0;
		for(int i=1;i<=len;i++)if(ans<p[i]){ans=p[i];pos=i;}
		//单纯求出最大回文串的长度
		cout<<ans-1<<endl;
		//下面几行是输出原来字符串中的回文串
		// for(int i=pos-p[pos]+1;i<pos+p[pos];i++){
		// 	if(i%2==0)cout<<s[i/2-1];
		// }
		// cout<<endl;
	}
	return 0;
}