码迷,mamicode.com
首页 > 编程语言 > 详细

马拉车算法

时间:2018-06-13 19:38:24      阅读:274      评论:0      收藏:0      [点我收藏+]

标签:字符   暴力   man   printf   har   int   pac   namespace   ace   

马拉车算法

给出一个只由小写英文字符a,b,c...y,z组成的字符串S,求S中最长回文串的长度。字符串长度len<=1e7

马拉车算法是用来求最长回文子串的。如果暴力的话是\(O(n^2)\)。能不能用到什么回文串的性质呢?

我们发现回文串内的任意一个串,都可以找到一个对于回文串对称中心对称的镜像串。用\(p[i]\)表示第i位的最大回文串长度/2。设\(R=max(j+p[j]), j<i\),也就是前面的回文串的最右边的位置,同时找到pos使得\(p[pos]=R\),也就是那个贡献最右回文串的中心位置。

技术分享图片

如果情况是第一种,那么显然\(p[i]=p[j]\)。如果情况是第二种,那么说明,i的回文串长度至少为\(p[j]\),但是会延伸到R以外,因此暴力判断R以外的红色部分即可。

由于R不停往右挪,因此时间复杂度是\(O(n)\)

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

const int maxn=3e7+5;
char s[maxn], s2[maxn];
int n;

int p[maxn], pos, R, ans;
int Manachar(char *s){
    pos=R=ans=0; memset(p, 0, sizeof(p));
    for (int i=1; i<n; ++i){
        if (i<R) p[i]=min(p[pos*2-i], R-i); else p[i]=1;
        for (; i-p[i]>=0&&s[i+p[i]]==s[i-p[i]]; ++p[i]);
        if (p[i]+i>R) R=p[i]+i, pos=i;
    }
    for (int i=0; i<n; ++i) ans=max(ans, p[i]);
    return ans;
}

int main(){
    scanf("%s", s); n=strlen(s);
    for (int i=0; i<n; ++i) s2[i<<1]='O', s2[i<<1|1]=s[i];
    n<<=1; s2[n++]='O';
    printf("%d\n", Manachar(s2)-1);
    return 0;
}

马拉车算法

标签:字符   暴力   man   printf   har   int   pac   namespace   ace   

原文地址:https://www.cnblogs.com/MyNameIsPc/p/9179331.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!