计算最长回文子串

以下理论性知识转自http://www.cnblogs.com/wuyiqi/archive/2012/06/25/2561063.html

友情链接（总结的很全很好）：https://segmentfault.com/a/1190000003914228

    其实原文说得是比较清楚的，只是英文的，我这里写一份中文的吧。
    首先：大家都知道什么叫回文串吧，这个算法要解决的就是一个字符串中最长的回文子串有多长。这个算法可以在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开始的）

下面是自己弄的模板：

#include<iostream>
#include<cstdio>
#define Min(x,y) x>y?y:x
using namespace std;
int p[150];
string str1,str2;
int len;
void pre()     //预处理
{
    str2[0]='$';
    str2[1]='#';
    len=str1.size();
    for(int i=0;i<len;i++)
    {
        str2[i*2+2]=str1[i];
        str2[i*2+3]='#';
    }
    len=len*2+2;
    str2[len]=0;
}

void pk()      //计算回文子串
{
    int mx=0,id;
    for(int i=1;i<len;i++)
    {
        if(mx>i)
        {
            p[i]=Min(p[id*2-i],mx-i);
        }
        else p[i]=1;
        for(;str2[i-p[i]]==str2[i+p[i]];p[i]++) {}
        if(mx<p[i]+i)
        {
            mx=p[i]+i;
            id=i;
        }
    }
}
int checkmax()    //找最长回文子串
{
    int ans=0,flag=0;
    for(int i=1;i<len;i++)
    {
        if(ans<p[i])
        {
            ans=p[i]; flag=i;
        }
    }
    return flag;
}
int main()
{
    int n;
    cin>>n;
    while(n--)
    {
        cin>>str1;
        pre();
        pk();
        int flag=checkmax();
        int l=p[flag]-1,count=0;
        for(int i=flag-l;count<l*2;i++)
        {
            count++;
            if(str2[i]=='#') continue;
            printf("%c",str2[i]);
        }
        puts(" ");
    }
    return 0;
}