Manacher算法 学习笔记

Manacher算法

Manacher算法主要是用来解决字符串回文串的问题的。
Manacher算法把长度为奇数和偶数的字符串一起考虑

算法流程

假设我们现在要求出字符串$s$的最长回文子串.
我们先对原串进行处理，把长度为奇数的和偶数的放在一起考虑。
其实就是在每两个字符串中间加入一个特殊的字符（只要是没有在原串中出现过的就可以，一般是#）

int Init()
{
    int len=strlen(s);
    int j=2;
    T[0]='$';
    T[1]='#';
    for(int i=0;i<len;i++)
    {
        T[j++]=s[i];
        T[j++]='#';
    }
    T[j]='\0';
    return j;
}

其中$j$是经过处理后的字符串的长度.
然后我们就可以开始对处理得到的字符串$T$进行处理。
Manacher算法还要求一个辅助数$Len[]$,$Len[i]$表示以字符$T[i]$为中心的最长回文字串的最右字符到$T[i]$的长度，比如以$T[i]$为中心的最长回文字串是$T[l,r]$,那么$Len[i]=r-i+1$。
显然对于$Len[i]$，我们有一个非常好的性质：
那就是$Len[i]-1$就是以$i$为中心的原串的最长回文子串的长度，（我们不是在原串中每隔两个字符又加了一个字符么）
既然$Len[]$有这么好的性质，那我们计算出$Len[]$数组就可以求出最长回文子串的长度了。
怎样求$Len[]$呢?
我们从左向右开始计算。
我们设$id$表示当前得到的最长回文子串的中心，$mx$为当前得到的最长回文子串的右端点。
对于$Len$数组的计算分为两种情况：
一.当前要计算的$i$点小于$mx$.
那么对于$i$的最长回文子串只有两种情况：一种是$i$关于id的对称点$j$的最长回文子串的长度，另一种是$mx-i$的长度。
对于第一种情况：

显然因为$mx`-mx$是以$id$为中心的最长回文子串，那么以$i$为中心的最长回文子串的长度一定等于，$i$关于$id$的对称点$j$的最长回文子串的长度。
对于第二种情况：
为什么还要$min(p[j],mx-i)$呢？
这是因为$mx`$左边和$mx$右边的字符是不一样的，有可能存在这样一种情况：

二.当前要计算的$i$点大于$mx$.
那么只能一位一位的计算了：

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
const int MAXN=11000300;
using namespace std;
int p[MAXN<<1];
char T[MAXN<<1],s[MAXN<<1];
int Init()
{
    int len=strlen(s);
    int j=2;
    T[0]='$';
    T[1]='#';
    for(int i=0;i<len;i++)
    {
        T[j++]=s[i];
        T[j++]='#';
    }
    T[j]='\0';
    return j;
}
int Manacher()
{
    int len=Init();
    int maxx=-1;
    int mx=0,id;
    for(int i=1;i<len;i++)
    {
        if(i<mx)
            p[i]=min(mx-i,p[2*id-i]);
        else
            p[i]=1;
        while(T[i+p[i]]==T[i-p[i]])
            p[i]++;
        if(i+p[i]>mx)
        {
            id=i;
            mx=i+p[i];
        }
        maxx=max(maxx,p[i]-1);
    }
    return maxx; 
}
int main()
{
    scanf("%s",s);
    printf("%d",Manacher());
    return 0;
}