Manachar's Algorithm

1、模板

#include<bits/stdc++.h>
using namespace std;
const int MAX=21000020;
char s[MAX],t[MAX<<1];
int p[MAX<<1],cnt=0,mid,mr;
void manachar()
{
    t[cnt++]='$';
    int len=strlen(s);
    for(int i=0;i<len;i++)
    {
        t[cnt++]='#';
        t[cnt++]=s[i];
    }
    t[cnt++]='#';
    t[cnt]='';
    
    for(int i=0;i<cnt;i++)
    {
        p[i]=i<mr?min(p[(mid<<1)-i],mr-i):1;
        for(;t[p[i]+i]==t[i-p[i]] and p[i]+i<cnt;p[i]++)
            if(p[i]+i>mr)
            {
                mid=i;
                mr=p[i]+i;
            }
    }
}
int main()
{
    scanf("%s",s);
    manachar();
    int ans=1;
    for(int i=0;i<cnt;i++)
        ans=max(ans,p[i]);
    printf("%d",ans-1);
    return 0;
}

2、一些结论

　　1、在往后推maxr时，此时的p[i]就是以maxr为结尾的最长回文串（例题：最长双回文串（有坑））

#include<bits/stdc++.h>
using namespace std;
const int MAX=1e5+10;
char pic[MAX],s1[MAX<<1],s2[MAX<<1];
int p1[MAX<<1],p2[MAX<<1],cnt=2;
int ans1[MAX<<1] ,ans2[MAX<<1],len;
void manachar(char * s,int * p,int * ans)
{
    int mid=0,mr=0;
    for(int i=0;i<cnt;i++)
    {
        p[i]=i<mr?min(p[mid*2-i],mr-i):1;
        for(;s[p[i]+i]==s[i-p[i]];p[i]++)
            if(i+p[i]>mr)
            {
                mr=i+p[i];
                ans[mr]=p[i];
                mid=i;
            }
    }
}
int main()
{
    scanf("%s",pic);
    len=strlen(pic);
    s1[0]=s2[0]='$';
    s1[1]=s2[1]='#';
    for(int i=0;i<len;i++)
    {
        s1[cnt]=pic[i];
        s2[cnt]=pic[len-i-1];
        s1[++cnt]='#';
        s2[cnt++]='#';
    }
    cnt++;
    s1[cnt]=s2[cnt]='';
    manachar(s1,p1,ans1);
    manachar(s2,p2,ans2);
    int j=cnt-4,ans=2;
    for(int i=3;i<cnt-2;i++)
    {
        if(s1[i]=='#')
            ans=max(ans,ans1[i]+ans2[j]);
        j--;
    }
    cout<<ans<<endl;
    return 0;
}

　　2、