题解/算法-manacher算法

题目链接:

luogu

题目描述:

给出一个只由小写英文字符a,b,c…y,z组成的字符串S,求S中最长回文串的长度.字符串长度为n

字符串长度len <= 11000000

题解:

用处:O(n)时间复杂度查最长回文子串

建议先看一遍代码

#include <iostream>  
#include <cstring>
#include <algorithm>  
using namespace std;

char s[1000]; //原数组
char s_new[2000]; //插入分隔字符后的新数组
int p[2000]; //Palindrome,记录每个位置的最长回文

int Init(){ //插入分隔字符#
    int len = strlen(s);
    s_new[0] = '$';
    s_new[1] = '#';
    int j = 2;

    for (int i = 0; i < len; i++){
        s_new[j++] = s[i];
        s_new[j++] = '#';
    }

    s_new[j] = '\0';  //字符串结束
    return j;  // 返回 s_new 的长度
}

int Manacher(){
    int len = Init();  //取得新字符串长度并完成向 s_new 的转换
    int max_len = -1;  //最长回文长度

    int mid;
    int mx = 0; //mx->range

    for (int i = 1; i < len; i++){
        if (i < mx)
            p[i] = min(p[2 * mid - i], mx - i);
        else
            p[i] = 1;

        //往两端扩展,判断是否为回文
        while (s_new[i - p[i]] == s_new[i + p[i]])  // 不需边界判断,因为左有'$',右有'\0'
            p[i]++;

        //每走一步i,都要和mx比较,希望mx尽可能的远,这样才能更有机会执行 if(i < mx)这句代码,从而提高效率
        if (mx < i + p[i]){
            mid = i;
            mx = i + p[i];
        }

        max_len = max(max_len, p[i] - 1);
    }

    return max_len;
}

int main(){
        scanf("%s", s);
        printf("%d\n", Manacher());
    }
    return 0;
}

由于回文串的中点可能在字符上也可能在字符中间,所以我们用'#’占满字符中间的位置

int Init(){ //插入分隔字符#
    int len = strlen(s);
    s_new[0] = '$';
    s_new[1] = '#';
    int j = 2;

    for (int i = 0; i < len; i++){
        s_new[j++] = s[i];
        s_new[j++] = '#';
    }

    s_new[j] = '\0';  //字符串结束
    return j;  // 返回 s_new 的长度
}

然后就很难理解了

我们从s[1]遍历到s[s_new.size],循环变量为i;
设置一个辅助变量mx,表示已经触及到的最右边的字符,还要设置一个辅助变量mid,表示包含maxright的回文串的中点所在的位置。
我们可以计算出在这个回文串中i的对称点是2 * mid - i

如果i<mx,可以先把p[i]更新为

p[i] = min(p[2 * mid - i], mx - i);

mx-i就是i到mx的距离
由于i的回文串可能会扩展到mx的右边,我们还要再加上这部分
而当i>mx时,我们无法得知关于p[i]的信息,只好从1开始遍历

 //往两端扩展,判断是否为回文
        while (s_new[i - p[i]] == s_new[i + p[i]])  // 不需边界判断,因为左有'$',右有'\0'
            p[i]++;

最后再更新mx和 max_len就可以了

 //每走一步i,都要和mx比较,希望mx尽可能的远,这样才能更有机会执行 if(i < mx)这句代码,从而提高效率
        if (mx < i + p[i]){
            mid = i;
            mx = i + p[i];
        }
        
 max_len = max(max_len, p[i] - 1);

如果没看懂,推荐一篇不错的题解传送门

写作时间:

2019-7-25

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值