POJ-3974 Palindrome（manacher算法）

题目链接：点击打开链接

经典的求最长回文子串题目，需要用到manacher算法。

manacher算法的具体介绍：http://blog.csdn.net/ggggiqnypgjg/article/details/6645824

#include<iostream>#include<cstdio>#include<cmath>#include<cstring>#include<algorithm>using namespace std;const int maxn=1000005;char s[maxn];char str[maxn*2];//用来存储s更新后的字符串int p[maxn*2];   //p[i]表示以str[i]为中心的最长回文串的半径int len1,len2;void init()  //初始化操作，更新s到str中，将原字符串两两间隔插入'#'变成奇数串，以便于处理{    len1=strlen(s);    str[0]='$'; //防止数组越界，字符串首端用'$'标记    str[1]='#';    for(int i=0;i<len1;i++)    {        str[i*2+2]=s[i];        str[i*2+3]='#';    }    str[2*len1+2]='\0';//同样是防止越界的操作，我个人习惯用'\0'}void manacher(){    init();    memset(p,0,sizeof p);    int id=0,mx=0;    //id<i，mx=i+p[id]表示在i之前回文串的最大半径距离i的最短距离    len2=strlen(str);    for(int i=1;i<len2;i++)    {        if(mx>i)            p[i]=min(p[2*id-i],mx-i); //2*id-i表示i关于id对称的点，根据对称性可得到此行的关系式        else p[i]=1;        while(str[i+p[i]]==str[i-p[i]]) p[i]++;//从i+p[i]继续向外延伸判断回文串是否能延长        if(i+p[i]>mx)        {            mx=i+p[i];            id=i;        }    }}int main(){    int cnt=1;    while(~scanf("%s",s))    {        if(s[0]=='E') break;        manacher();        int ans=0;        for(int i=0;i<len2;i++)            ans=max(ans,p[i]);        printf("Case %d: %d\n",cnt++,ans-1);    }    return 0;}