回文串的Manacher算法

Manacher算法较传统算法的优化之处在于它对每个回文中心寻找回文半径的时候并不是都从半径为1开始找的，而是利用前面已经完成的任务，寻找一个初始的开始搜索的半径大小，复杂度是线性的。

参考博客：https://www.cnblogs.com/z360/p/6375514.html

下面附上hdoj的3068模板：

#include<bits/stdc++.h>
using namespace std;
typedef unsigned int ui;
typedef long long ll;
typedef unsigned long long ull;
#define pf printf
#define mem(a,b) memset(a,b,sizeof(a))
#define prime1 1e9+7
#define prime2 1e9+9
#define pi 3.14159265
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define scand(x) scanf("%llf",&x) 
#define f(i,a,b) for(int i=a;i<=b;i++)
#define scan(a) scanf("%d",&a)
#define dbg(args) cout<<#args<<":"<<args<<endl;
#define inf 0x3f3f3f3f
#define maxn 110005
int n,m,t;
char s[maxn];
char tmp[maxn<<1];
int l[maxn<<1];//转换后的字符串以及最大回文半径 
int manacher(char *st,int len) 
{
    int t=0;
    tmp[t++]='$';
    tmp[t++]='#';
    f(i,0,len-1)
    {
        tmp[t++]=s[i];
        tmp[t++]='#';
    }
    tmp[t]=0;
    int mr=0,ans=0,pos=0;//最长回文右端、最长回文长度以及回文中心 
    f(i,0,t-1)
    {
        l[i]=mr>i?min(mr-i,l[2*pos-i]):1;
        //i>=mr时需要以i为新的回文中心开始搜索回文长度  
        while(tmp[i-l[i]]==tmp[i+l[i]])l[i]++;
        if(i+l[i]>mr)//更新最大回文右端点 
        {
            mr=i+l[i];
            pos=i;
        }
        ans=max(ans,l[i]);
    }
    return ans-1;//每次回文串长度为2*k+1时，多余字符为k+1个，有k个属于原串的字符 
}
int main()
{
    //freopen("input.txt","r",stdin);
    //freopen("output.txt","w",stdout);
    std::ios::sync_with_stdio(false);
    while(scanf("%s",s)!=EOF)
    {
        int len=strlen(s); 
        pf("%d
",manacher(tmp,len));
    }
    
 }