Manacher算法(马拉车)

什么是Manacher
Manacher算法是用来求一个字符串中的最长回文串的算法。
时间复杂度为:O(2*n) 比较高效。

Manacher算法的原理
如果直接对原有字符串进行操作的话,需要考虑回文串的长度是奇数还是偶数(偶数长度的回文串没有中心字符,不易处理)。暴力匹配的复杂度是:O(n^3)。所以有了马拉车算法。

1. 第一步是对原有字符串进行预处理。

在任意两个字符之间及字符串的首尾位置插入一个在所求字符串中不会出现的字符,一般使用“#”,在整个字符串的首位置插入一个其他字符,这里使用“@”(因为在后面的while循环中在匹配字符的时候可能会越界)(预处理之后的字符串命名为str)
eg:原串:abba
操作之后:@a#b#b#a#

2. 第二部是对新得到的字符串进行操作。

用id记录上一次操作的中间字符的位置,mx标记上一次的最长子串的最右端,用一个Len[i]数组去存第i个位置到mx位置的长度,然后依次去递推。
首先把mx初始化0,然后开始从i=1开始遍历str。
(1) 当i>=mx的时候,也就是i在mx前面的时候,就让Len[i] = 1,表示在i之前没有回文串出现,所以让Len[i]从1开始,然后进入while循环,查找以当前点为中心的最长回文串的长度,依次比较id左边1,2,3…和右边1,2,3…是否相等,相等的话就让Len[i]++,否则就跳出。然后比较i+len[i]与mx的位置,超过的话就更新一下mx的值,因为mx是此次操作的最长回文子串的最右端。
(2) j是以id为中点的i的对称点,因为j在之前操作中就求出来了,所以Len[j]是已知的又因为回文串的性质可以知道Len[j]<=Len[i]的,my是mx的对称点。
(i) 对于i<mx并且Len[j]<mx-i的,表示j位置的最长回文子串的长度是在id位置的最长回文子串的范围内的,就i的Len[2*id - i]也就是Len[mx-i]的地方开始匹配
在这里插入图片描述
(ii) 当i<mx且Len[j] > mx - i的时候,说明i位置的子串长度超过了mx,但mx以外的地方还没有遍历到,所以我们就从mx-i也就是mx的位置开始对i匹配。
在这里插入图片描述
模板

#include <iostream>
#include <cstdio>
#include <cstring>
#define Min(a,b) a>b?b:a
#define Max(a,b) a>b?a:b
using namespace std;
int Len[3000005];
char str[3000005],s[3000005];
int n,mx,id,len;

void init()
{
    memset(str,0,sizeof(str));
    int k=0;
    str[k++] = '$';
    for(int i=0; i<len; i++)
    {
        str[k++]='#';
        str[k++]=s[i];
    }
    str[k++]='#';
    len=k;
}///预处理操作,得到str数组

int manacher()
{
    Len[0] = 0;
    int sum = 0;
    mx = 0;
    for(int i=1; i<len; i++)
    {
        if(i < mx)
            Len[i] = min(mx - i, Len[2 * id - i]);///当i<mx的时候,直接取mx-i和Len[2 * id - i]的最小值
        else
            Len[i] = 1;///i>=mx的时候代表当前位置已经不被上一状态的最长回文串涵盖在内,需要更新
        while(str[i - Len[i]]== str[i + Len[i]])
            Len[i]++;///寻找最长的Len[i]
        if(Len[i] + i > mx)
        {
            mx = Len[i] + i;///更新mx
            id = i;///更新id
            sum = max(sum, Len[i]);
        }
    }
    return (sum - 1);///sum是处理后的数组的最长回文串的一半,(sum-1)为原字符串的最长回文串长度
}

int main()
{
    scanf("%d",&n);
    while(n--)
    {
        scanf("%s",s);
        len = strlen(s);
        init();
        int ans = manacher();
        printf("%d\n",ans);
    }
    return 0;
}

资源链接:
https://blog.csdn.net/Charles_Zaqdt/article/details/79747073

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值