hihocoder1032 最长回文子串

传说中的Manacher算法，算法的核心就在这一句上了：p[i] = min(p[2*id-i], p[id] + id - i);

这个算法要解决的就是一个字符串中最长的回文子串有多长。这个算法可以在O（n）的时间复杂度内既线性时间复杂度的情况下，求出以每个字符为中心的最长回文有多长，这个算法有一个很巧妙的地方，它把奇数的回文串和偶数的回文串统一起来考虑了。这一点一直是在做回文串问题中时比较烦的地方。这个算法还有一个很好的地方就是充分利用了字符匹配的特殊性，避免了大量不必要的重复匹配。

算法大致过程是这样。先在每两个相邻字符中间插入一个分隔符，当然这个分隔符要在原串中没有出现过。一般可以用‘#’分隔。这样就非常巧妙的将奇数长度回文串与偶数长度回文串统一起来考虑了（见下面的一个例子，回文串长度全为奇数了），然后用一个辅助数组P记录以每个字符为中心的最长回文串的信息。P［id］记录的是以字符str［id］为中心的最长回文串，当以str［id］为第一个字符，这个最长回文串向右延伸了P［id］个字符。

原串： w aa bwsw f d

新串： # w# a # a # b# w # s # w # f # d #

辅助数组P： 1 2 1 2 3 2 1 2 1 2 1 4 1 2 1 2 1 2 1

这里有一个很好的性质，P［id］-1就是该回文子串在原串中的长度（包括‘#’）。如果这里不是特别清楚，可以自己拿出纸来画一画，自己体会体会。当然这里可能每个人写法不尽相同，不过我想大致思路应该是一样的吧。

好，我们继续。现在的关键问题就在于怎么在O（n）时间复杂度内求出P数组了。只要把这个P数组求出来，最长回文子串就可以直接扫一遍得出来了。

由于这个算法是线性从前往后扫的。那么当我们准备求P［i］的时候，i以前的P［j］我们是已经得到了的。我们用mx记在i之前的回文串中，延伸至最右端的位置。同时用id这个变量记下取得这个最优mx时的id值。（注：为了防止字符比较的时候越界，我在这个加了‘#’的字符串之前还加了另一个特殊字符‘$’，故我的新串下标是从1开始的）

if( mx > i)

p[i]=MIN( p[2*id-i], mx-i);

就是当前面比较的最远长度mx>i的时候，P［i］有一个最小值。这个算法的核心思想就在这里，为什么P数组满足这样一个性质呢?

（下面的部分为图片形式）

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define M 1000005
string s;
int f[M*2+4];
int solve(string s)
{
string s1;
s1.resize(s.length()*2+2);
s1[0]='$';
s1[1]='#';
for(int i=0;i<s.length();i++)
{
s1[(i+1)<<1]=s[i];              //从1开始的
s1[(i+1)*2+1]='#';
}
memset(f,0,sizeof(f));
int ans=0,mx=0;                           //mx即为上面说的最右端的位置
for(int id=0,i=1;i<s1.length();i++)
{
if(mx>i)
f[i]=min(f[id*2-i],f[id]-i+id);
else
f[i]=1;
while(s1[i+f[i]]==s1[i-f[i]])           //更新i位置的f[i]
f[i]++;
if(f[i]+i>mx)
{
mx=f[i]+i;
id=i;
}
ans=max(f[i],ans);
}
return ans-1;
}
int main()
{
int t;
scanf("%d",&t);
while(t--)
{
cin >> s;
printf("%d\n",solve(s));
}
return 0;
}