HDU 3068 manacher算法(O(n)的最长回文串算法)

宿舍的交换机坏了，因为之前还没开学，所以就没人修,然后就拖到现在寝室才有网。然后这几天因为省赛选拔赛所以心好累
/**************************************************************

manacher是很久之前就听说的想学的，今天晚上因为选拔完了就闲了看的

数组:
s 字符原串
str 通过在每两个数组元素之间添加’#’形成的新串
p[i] 代表以str[i]为中心的回文串最大长度
元素：
mx 当前位置之前的回文串能达到最右边的位置
id 当前形成mx的回文串的中心点

/此处将末尾改为0纯粹是因为想跑FOR方便而已/
众所周知，回文串长度可以为奇，可以为偶，所以进行回文串判断过程中就会重复一些不必要的过程，manacher算法通过每两个元素之间添加新无关元素来实现长度为 奇+偶=奇来统一长度减少时间复杂度

一般方法，跑到当前i节点然后向两边扩展，最坏时间复杂度为o(n^2)
manacher算法通过运用回文串的对称性来达到跑到i节点，可以从k长度开始扩展,k=与i点关于id点对称的j点的p[j],但是考虑到一种情况

所以,k=min(j-x,p[j])

#include<stdio.h>#include<string.h>#include<iostream>#include<algorithm>#include<math.h>#include<queue>#include<stack>#include<string>#include<vector>#include<map>#include<set>using namespace std;#define rfor(i,a,b) for(i=a;i<=b;++i)#define lfor(i,a,b) for(i=a;i>=b;--i)#define sfor(i,a,h) for(i=h[a];i!=-1;i=e[i].next)#define mem(a,b) memset(a,b,sizeof(a))#define mec(a,b) memcpy(a,b,sizeof(b))#define cheak(i) printf("%d ",i)#define min(a,b) (a>b?b:a)#define max(a,b) (a>b?a:b)#define inf 0x7ffffffftypedef long long LL;#define maxn 500005char str[maxn],s[maxn];int p[maxn];int manacher(char* s,int len){    int i,mx=0,l,id,ans=0;    str[0]='@';str[1]='#';    for(i=1;i<=len;++i)    str[i<<1]=s[i],str[(i<<1)+1]='#';    str[len*2+1]='#';str[len*2+2]=0;    for(i=1;str[i];++i)    {        if(mx>i) p[i]=min(p[2*id-i],mx-i);        else p[i]=1;        while(str[i+p[i]]==str[i-p[i]])        p[i]++;        if(i+p[i]>mx)        {            mx=i+p[i];            id=i;        }        ans=max(p[i],ans);    }    return ans-1;}int main(){    int len,i;    while(~scanf("%s",s+1))    {        len=strlen(s+1);        int ans=manacher(s,len);        printf("%d\n",ans);    }    return 0;}