Manacher算法+注释

Manacher算法是用来求一个字符串中最长回文串的算法。

考虑暴力求最长回文串的做法：

暴力枚举字符串中的所有字串判断是否回文，然后求最大值。

时间复杂度O(n^3)，考虑优化。

我们从枚举所有字串改成枚举所有回文串的对称轴，向左右扩展直到不相等，得到最长回文串。

优化到O(n^2)，还是不够优秀。

于是我们引出Manacher算法。

先向字符串s中插入特殊字符得到字符串str，这样我们就不用讨论字符串长度是奇是偶了。

用一个辅助数组p表示每个点可以扩展出去的最长回文长度

从str[1]扫到str[strlen(str)]，再设置两个变量mr表示已触及的最右边的字符，mid表示包含mr的回文串的对称轴位置。

当i属于(mid,mr)时，显然i关于mid的对称点是(mid<<1)-i（中点坐标公式简单推一下），由于回文串对称串的全等性，我们令p[i]=p[(mid<<1)-i]，然后接着尝试扩展：str[i+p[i]]==str[i-p[i]]（前后是否对称），p[i]++

若i>mid，我们就设置mid=i,mr=当前扩展到的最右字符。

给出代码结束本篇博客

#include<bits/stdc++.h>
using namespace std;
inline int read(){
    int data=,w=;char ch=;
    while(ch!='-' && (ch<''||ch>''))ch=getchar();
    if(ch=='-')w=-,ch=getchar();
    while(ch>='' && ch<='')data=data*+ch-'',ch=getchar();
    return data*w;
}
const int maxn=5e7+;
int n,p[maxn],ans;
char s[maxn],str[maxn];
void init(){
    str[]=str[]='$';//ccf喜欢这个
    for(int i=;i<n;i++)
        str[(i<<)+]=s[i],str[(i<<)+]='$';
    n=(n<<)+;
    str[n]=;
}
void Manacher(){
    int mr=,mid;
    for(int i=;i<n;i++){
        if(i<mr)
            p[i]=min(p[(mid<<)-i],p[mid]+mid-i);
        else p[i]=;
        for(;str[i+p[i]]==str[i-p[i]];p[i]++);
        if(i+p[i]>mr)
            mr=p[i]+i,mid=i;
    }
}
int main(){
    scanf("%s",s);
    n=strlen(s);
    init();Manacher();
    ans=;
    for(int i=;i<n;i++)
        ans=max(ans,p[i]);
    printf("%d\n",ans-);
    return ;
}

下一篇更新一些数论知识