Renatus的博客

不管到了什么时候,人们最缺少的都是时间

Manacher算法详解

Manacher算法详解


RT,Manacher算法今天第一次写,有不少细节,也容易写错,下面进行一个小小的总结

原理

首先,我们应该已经看到过很多有关回文串的题目了,当我们不会Manacher的时候,我们往往是用回文自动机后缀数组代替的,这时候我们会发现有些题它并不能卡过,因为后缀数组是O(nlogn)的,而且编程复杂。这时候我们就需要Manacher算法,它是一个O(n)的算法,而且编程简单

它的主要思想是先将奇偶回文串情况讨论的部分去掉,我们可以通过在每个字符间,包括开头与末尾,插入特殊字符来达到这一点

然后,我们通过扩展KMP的思想,建立一个最远指针,它表示当前的回文串最远达到的地方,这时我们就可以利用回文的性质来加速,求出一个数组f,它表示以i为中心的回文串的半径。由于最远指针最多被推动n次,所以Manacher的复杂度是O(n)

代码

下面是代码

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define maxn 200005
using namespace std;
int f[maxn],len1,pre[maxn*2],suf[maxn*2],len;
char txt[maxn],s[maxn*2];
void Manacher(){
    for(int i=0;i<len1;i++){
        s[i<<1]='@';
        s[(i<<1)+1]=txt[i];
    }
    s[(len1)<<1]='@';
    len=2*len1+1;
    int mx=0,p=0;
    for(int i=0;i<len;i++){
        if(mx>i) f[i]=min(mx-i,f[2*p-i]);
        else f[i]=0;
        while(i+f[i]+1<len&&i-f[i]-1>=0&&s[i+f[i]+1]==s[i-f[i]-1]) f[i]++;
        if(f[i]+i>mx){
            mx=f[i]+i;
            p=i;
        }
    }
    return;
}
int main(){
    freopen("input.txt","r",stdin);
    freopen("output.txt","w",stdout);
    while(scanf("%s",txt)==1){
        len1=strlen(txt);
        Manacher();
        int k=0,op=0;
        for(int i=0;i<len;i++){
            if(i&1)
                f[i]=((f[i]>>1)<<1)+1;
            else
                f[i]=((f[i]+1)>>1)<<1;
            k=max(k,f[i]);
        }
        printf("%d\n",k);
    }
    return 0;
}

它的作用是每次求出一个字符串的最长回文子串

细节

总结一下

1.起点与终点都要插入特殊字符

2.输出时要考虑奇偶情况,即注意转换回原串回文串时的处理,注意不要先求k的最大值,最后再处理k

3.数组要足够大

题目

最近忙,待补

阅读更多
版权声明:本文为博主原创文章,未经博主允许不得转载。 https://blog.csdn.net/Goseqh/article/details/60140765
文章标签: Manacher 字符串
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

不良信息举报

Manacher算法详解

最多只允许输入30个字

加入CSDN,享受更精准的内容推荐,与500万程序员共同成长!
关闭
关闭