最长回文串O（n）算法 hdu3068

//http://bbs.dlut.edu.cn/bbstcon.php?board=Competition&gid=23474 一个奇妙的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就是该回文子串在原串中的长度（包括‘#’）。

 另外法是线性从前往后扫的。那么当我们准备求P［i］的时候，i以前的P［j］我们是已经得到了的。我们用mx记在i之前的回文串中，延伸至最右端的位置。同时用id这个变量记下取得这个最优mx时的id值。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;
#define N 320500

char str[N],a[N];
int p[N];
int main ()
{
    while(scanf("%s",str)!=EOF)
    {
        a[0]='$';a[1]='#';
        int n=strlen(str);
        for(int i=0;i<n;++i)
        {
            a[2*i+2]=str[i];
            a[2*i+3]='#';
        }
        //memset(p,0,sizeof(p));
        n=2*n+2;
        int mx=1,id=0,ans=0;
        for(int i=0;i<n;++i)
        {
            if(mx>i)
                p[i]=min(p[2*id-i],mx-i);
            else p[i]=1;
            while(a[i-p[i]]==a[i+p[i]]) p[i]++;
            if(i+p[i]>mx)
            {
                id=i;
                mx=i+p[i];
            }
            ans=max(ans,p[i]);
        }
        printf("%d\n",ans-1);
    }
    return 0;
}