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;
}