manacher模板

给定一个字符串,要求O(n)时间求出其最长回文串长度。

首先我们是会O(n^2)的暴力的,就是枚举每个字符作为对称中心,再枚举相同的相邻字符作为对称中心,然后求得答案。

能不能优化呢?

我们发现,当一个字符串是回文串时,它满足完全对称,比如:

abadaba

设p[i]为以i为中心的回文半径,则p为:

1 2 1 4 1 2 1

然后发现右半边竟然可以O(1)出解!

然后就可以有第一步优化。

回到上面,按中心分类,回文分为两种。

能不能将他们都归为满足优化1的那一种?

答案是能的。在每个字符前后加没有意义的字符,比如 ' # ';

这就是优化2。

就有了manacher算法:

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
char a[51000050],b[51000050];
int len;
int p[102000050];
void manacher()
{
    int mx = 0,mid =0;
    for(int i=1;i<=((len<<1)|1);i++)
    {
        if(i<mx)p[i] = min(p[2*mid-i],mx-i+1);
        else p[i]=1;
        while(b[i-p[i]]==b[i+p[i]])p[i]++;
        if(i+p[i]-1>mx)
        {
            mx = i+p[i]-1;
            mid = i;
        }
    }
}
int main()
{
    scanf("%s",a+1);
    len = strlen(a+1);
    b[0]='!',b[(len<<1)+2]='^';
    b[1]='#';
    for(int i=1;i<=len;i++)
    {
        b[i<<1] = a[i];
        b[i<<1|1] = '#';
    }
    manacher();
    int ans = 0;
    for(int i=1;i<=((len<<1)|1);i++)
    {
        ans=max(ans,p[i]);
    }
    printf("%d\n",ans-1);
    return 0;
}

 

posted @ 2018-09-18 21:16  LiGuanlin  阅读(124)  评论(0编辑  收藏  举报