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.数组要足够大

题目

最近忙,待补

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值