回文串算法Manacher

最新推荐文章于 2024-12-07 14:19:22 发布

cyl纤云弄巧

最新推荐文章于 2024-12-07 14:19:22 发布

阅读量507

点赞数

CC 4.0 BY-SA版权

分类专栏： poj 字符串文章标签： Manacher

本文链接：https://blog.csdn.net/Since_natural_ran/article/details/78906189

poj 同时被 2 个专栏收录

171 篇文章

订阅专栏

字符串

2 篇文章

订阅专栏

本文介绍了一种高效的回文串查找算法——Manacher算法，并详细解释了其核心思想及实现过程。通过巧妙地处理字符串，该算法能有效地找到字符串中所有最长回文子串。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

回文串算法Manacher

首先用一个非常巧妙的方式，将所有可能的奇数/偶数长度的回文子串都转换成了奇数长度：在每个字符的两边都插入一个特殊的符号。比如 abba 变成 #a#b#b#a#， aba变成 #a#b#a#。为了进一步减少编码的复杂度，可以在字符串的开始加入另一个特殊字符，这样就不用特殊处理越界问题，比如$#a#b#a#。

下面以字符串12212321为例，经过上一步，变成了 S[] = “$#1#2#2#1#2#3#2#1#”;

然后用一个数组 P[i] 来记录以字符S[i]为中心的最长回文子串向左/右扩张的长度（包括S[i]），比如S和P的对应关系：

S # 1 # 2 # 2 # 1 # 2 # 3 # 2 # 1 #
P 1 2 1 2 5 2 1 4 1 2 1 6 1 2 1 2 1
(p.s. 可以看出，P[i]-1正好是原字符串中回文串的总长度）

下面计算P[i]，该算法增加两个辅助变量id和mx，其中id表示最大回文子串中心的位置，mx则为id+P[id]，也就是最大回文子串的边界。

这个算法的关键点就在这里了：

Mp[i] = mx > i ? min(Mp[2*id-i],mx-i):1;

其实就是对称点和边界的问题。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int maxn = 1000010;

char Ma[maxn*2];
int Mp[maxn*2];
char s[maxn];

void Manacher(char s[],int len)
{
    int l = 0;
    Ma[l++] = '$';
    Ma[l++] = '#';
    for(int i = 0;i < len; i++) {
        Ma[l++] = s[i];
        Ma[l++] = '#';
    }
    Ma[l] = 0;
    int mx = 0,id = 0;
    for(int i = 0;i < l; i++) {
        Mp[i] = mx > i ? min(Mp[2*id-i],mx-i):1;
        while(Ma[i+Mp[i]] == Ma[i-Mp[i]]) Mp[i]++;
        if(i + Mp[i] > mx) {
            mx = i + Mp[i];
            id = i;
        }
    }
}
int main()
{
    freopen("in.txt","r",stdin);
    int ncase = 1;
    while(scanf("%s",s) == 1) {
        if(strcmp(s,"END") == 0) break;
        int len = strlen(s);
        Manacher(s,len);
        int ans = 0;
        for(int i = 0;i < 2*len+2; i++)
            ans = max(ans,Mp[i]-1);
        printf("Case %d: %d\n",ncase++,ans);
    }

    return 0;
}