[hihoCoder] #1032 : 最长回文子串

这是我看到的做回文子串最神奇的一个方法（其实之前只看过暴力）

详解：http://www.cnblogs.com/easonliu/p/4454213.html

用O(n)的时间复杂度解决这个问题，突然有种做kmp的感觉。

算法精髓是p[i] = min（p[2*id-i],p[id]+id-i）

  算法大致过程是这样。先在每两个相邻字符中间插入一个分隔符，当然这个分隔符要在原串中没有出现过。一般可以用‘#’分隔。这样就非常巧妙的将奇数长度回文串与偶数长度回文串统一起来考虑了（见下面的一个例子，回文串长度全为奇数了），然后用一个辅助数组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 <bits/stdc++.h>using namespace std;int N;string s;void solve(){string s1;s1.resize(2*s.size()+2);s1[0]='$';s1[1]='#';int len=s.size();for(int i=0;i<len;i++){    s1[(i+1)<<1]=s[i];    s1[((i+1)<<1)+1]='#';} vector<int> p(s1.size(), 0);int res=0;len=s1.size();for(int id=0,i=1;i<len;i++){    if(p[id]+id>i)p[i]=min(p[2*id-i],p[id]+id-i);    else p[i]=1;    while(s1[i+p[i]]==s1[i-p[i]])p[i]++;    if(i+p[i]>id+p[id])id=i;    res=max(res,p[i]);}printf("%d\n",res-1);}int main(){cin>>N;getchar();while(N--){    cin>>s;    solve();}    return 0;}