manacher(马拉车)

作用

可以在O(n)的时间计算出以每个字符或者两个字符为中心的最长回文串的长度,也可以凭借此求出最长回文串等.
相比于回文自动机,manacher侧重于求以每个或两个字符为中心的最长回文串,而回文自动机则侧重于求以每个字符为左右端点的最长回文串,而且manacher更加好写.

实现方法

首先我们在字符串两两字符之间加入永远不可能出现的字符,比如:!@#$%^&*,这样原串中的所有回文串都有一个中心.
为了后面操作方便,再在头尾加上另外两个永远不会出现的不同字符.
利用回文串的性质我们可以发现,若存在回文串l~r,中心为mid,则对于i(mid < i < r)位置的字符,以其为中心的回文串长度至少为min(ans[mid*2-i],(r-i) *2+1),这样就可以快速求解了.

代码(以POJ - 3974 Palindrome为例)

#include<iostream>
#include<cstdio>
#include<cstring>
#define N 1001000
using namespace std;

int n,ans,r,mid,tmp[N<<1],T;
char str[N],a[N<<1];

int main()
{
    int i,j,t;
    for(;;)
    {
        scanf("%s",str+1);
        if(str[1]=='E') return 0;
        for(i=1,t=strlen(str+1);i<=t;i++)
        {
            a[i<<1]=str[i];
            a[(i<<1)+1]='*';
        }
        n=2*t+1;
        a[0]='&',a[1]='*',a[n+1]='#';
        ans=r=0;
        for(i=1;i<=n;i++)
        {
            t= i<r?min(r-i+1,tmp[mid*2-i]):1;
            for(;a[i+t]==a[i-t];t++);
            tmp[i]=t;
            ans=max(ans,t-1);
            if(i+t-1>r)
            {
                r=i+t-1;
                mid=i;
            }
        }
        printf("Case %d: %d\n",++T,ans);
    }
}