manacher是一种优秀的,可以在O(n)时间复杂度内求最长回文子串个数的算法,又叫做“马拉车”
先膜一下大佬Orz,我是看ZigZagK学长的博客看懂的,神奇的传送门 。
预处理
大家都知道求回文串时可以枚举中心对称点,然后向左右两边暴力拓展,这时候,中心点有时在字符上,有时又在相邻字符中间,这比较烦人。其实解决方法很简单,在字符串中,每两个字符中间都插入一个无关紧要的字符,比如#
,%
之类的这样子所有的回文子串长度就一定是奇数的了。(是不是特别简单,读入时预处理一下就可以了)注意头和尾也要塞一个字符,打头儿的那个新加进来的字符我喜欢放到位置0上。
讨论情况
假设我们已经求好了以1~i-1为中心的回文子串长,现在要求以i为中心的回文子串长,可以怎么求呢?
先维护一个pos,表示以1~i-1为中心的回文子串的右端点最远的点是哪个,而以pos为中心的回文子串的右端点为R,左端点为L。
画个图来理解一下。
上图是一个最容易处理的形式,因为以j为中心的回文串是被以pos为中心的回文串包含的,且i和j以pos为中心对称,所以以i为中心的回文串与以j为中心的回文串一定以pos为中心对称,所以以i为中心最大的回文串就很好求了。
还有一种情况,是以j为中心的回文串的左端点在小于等于L,这时,i的右端点只能从R开始,向右暴力推进了。
如果i在R右边呢?i最多为R+1,暴力推进呗。
一个不识时务的分割线(雾
主体思想就讲完了,没了,真的没了,虽然前面很多时暴力,但真的没了。所以怎么证明复杂度为O(n)?因为R只会往右推进,而且每个点R只会经过一次。
模板代码
lougu3805【模板】manacher算法 水过 (逃
#include<cstdio>
#include<algorithm>
using namespace std;
int len=-1,R,pos,f[22000005],ans=1;
char s[22000005];
int main(){
freopen("manacher.in","r",stdin);
freopen("manacher.out","w",stdout);
char ch=getchar();
while(ch>='a'&&ch<='z'){//读入字符串
s[++len]='#';s[++len]=ch;
ch=getchar();
}s[++len]='#';
for(int i=1;i<=len;i++){
if(i>R){
pos=i;R=i;
for(int j=1;i+j<=len&&i-j>=0;j++) if(s[i+j]==s[i-j]) R=i+j;else break;//暴力推进
f[i]=R;
ans=max(ans,f[i]-i);
continue;
}
int L=pos*2-R,j=pos*2-i;
if(j*2-f[j]<=L){//一定是<=而不是<,要注意了
pos=i;
for(int k=R-i+1;i+k<=len&&i-k>=0;k++) if(s[i+k]==s[i-k]) R=i+k;else break;
f[i]=R;
}else f[i]=i+f[j]-j;
ans=max(ans,f[i]-i);
}
printf("%d\n",ans);
return 0;
}