POJ-3974-Palindrome- Manacher 马拉车算法（On寻找最长回文串）

http://poj.org/problem?id=3974

马拉车算法的讲解，这个地方讲得很详细：http://www.cnblogs.com/grandyang/p/4475985.html

很简单的又很有威力的一个算法。

实现非常简单

该算法主要步骤就两步：

1、通过加不影响结果的分隔符#,使得不管奇偶串通通变成 奇数串

2、通过利用对称的关系  线性复杂度求得P数组

则最长回文串便是  max_p[i]-1了 （p[i]记录的是以i为中心的最长回文半径，但是由于处理后的数组中存在的#可以证得p[i]-1刚好是 去掉#后的回文串长度）

如果最长回文串的中心为#，表示这是一个偶数串，否则为奇数串

#include <cstdio>
#include <cmath>
#include <cstring>
#include <string>
#include <algorithm>
#include <queue>
#include <map>
#include <stack>
#include <vector>
#include <iostream>
using namespace std;
int min(int a,int b){return a<b?a:b;}

char tm[1000000+5];
char t[2000000+5];
int p[2000000+5];
int ok=0;

int Manacher(char *t)
{
int ans_len=1,ans_center=0;
int mx=0,id,i;
for (i=1;i<=ok;i++)
{
if (mx>i)
p[i]=min(p[2*id-i],mx-i);
else
p[i]=1;
while(t[i-p[i]]==t[i+p[i]])
p[i]++;
if (i+p[i]>mx)
{
mx=i+p[i];id=i;
if (p[i]>ans_len)
ans_len=p[i],ans_center=i;
}
}
return ans_len-1;
}

int main()
{
int i,j,k;
int cnt=1;
while(scanf("%s",tm+1)!=EOF)
{
if (strcmp(tm+1,"END")==0) break;
int len=strlen(tm+1);
ok=0;
t[0]='$';
for (i=1;i<=len;i++)
{
t[++ok]='#';
t[++ok]=tm[i];
}
t[++ok]='#';
t[ok+1]=0;
int ans=Manacher(t);
printf("Case %d: %d\n",cnt++,ans);

}

return 0;

}