最长回文子串

关键词：manacher算法

预处理以每个字符为中心的最长回文半径p[i].mx是前以i个元素为中心回文串的最右边界，id是该回文串的中心下标。

则2*id-i是i关于id的对称点。关键代码如下：

if(mx>i) p[i]=min(mx-i,p[2*id-i]);
else p[i]=1;

但是这种方法只能求长度为奇数的回文串长度，于是将原串扩增，wwaww->@#w#w#a#w#w#\0，用#将原串包围，并添加@和\0到首尾。新串的回文串只能是奇数

结论：新串中的p[i]-1就是原串回文长度

由例子可知，wwaww的回文串最长为5，新串中最长回文串是#w#w#a#w#w# p[6]=6,ans=p[6]-1,符合结论。

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<algorithm>
#include<map>
#include<vector>
#include<queue>
#include<stack>
#define ll long long
#define sf scanf
#define pf printf
#define INF 0x3f3f3f3f
#define maxn 110010
#define mem(a,b) memset(a,b,sizeof(a))
#define lowbit(x) x&(-x)
const ll mol=1000000007;
using namespace std;

char s[maxn<<2],str[maxn<<2];
int p[maxn<<2],ans;

int manacher(){
    mem(p,0),mem(str,0),ans=0;
    int l=strlen(s);
    str[0]='@';
    for(int i=0;i<l;i++){ str[2*i+2]=s[i],str[2*i+1]='#'; }
    str[2*l+1]='#',str[2*l+2]='\0';
    int mx=0,id;
    for(int i=1;str[i];i++){
        if(mx>i) p[i]=min(mx-i,p[2*id-i]);
        else p[i]=1;
        for(;str[i+p[i]]==str[i-p[i]];) p[i]++;
        if(i+p[i]>mx) { mx=i+p[i],id=i; }
        if(p[i]>ans) ans=p[i];
    }
    for(int i=1;str[i];i++) printf("%c %d\n",str[i],p[i]);
    return ans-1;
}

int main(){
    //freopen("a.txt","r",stdin);
    while(scanf("%s",s)!=EOF){
        printf("%d\n",manacher());
    }
}