POJ3974-Manacher算法--O（n）回文子串算法

题目：题目链接

题目意思：题意应该很好理解，就是求一个字符串中，最长的回文字串的长度是多少？字符串长度 1000000 。暴力T。看了神牛的代码，知道了这个Manacher算法--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开始的）

对于字符串长度大的，我们就尽量要使用O（n）的算法：

#include <iostream>
#include <cstdio>
#include <string>
#include <string.h>
#include <map>
#include <vector>
#include <cstdlib>
#include <cmath>
#include <algorithm>
#include <queue>
#include <set>
#include <stack>
#define M 10000005
using namespace std;

string s,tmp;
int casen,ans,rad[M];
int main()
{
    while(cin>>tmp)
    {
        if("END"==tmp)
            break;
        cout<<"Case "<<++casen<<": ";
        int n = tmp.size();
        s="";
        memset(rad, 0, sizeof(rad));
        for(int i=0; i<n; i++)//填充特殊字符
        {
            s += "#";
            s += tmp[i];
        }
        s += '#';
        n = n*2+1;
        int i=0, j=1, k;
        while(i < n)
        {
            while(i - j >= 0 && i+j < n && s[i-j]==s[i+j])
                j++;
            rad[i] = j-1;//rad[i]记录的是以i为中心左右做多可以扩展rad[i]个字符。
            k=1;
            while(k<=rad[i] && rad[i]-k!=rad[i-k])
            {
                rad[i+k] = min(rad[i-k],rad[i]-k);//DP
                k++;
            }
            i+=k;
            j=max(j-k,0);
        }
        for(i = 0; i < n; i++)
            ans = max(ans, rad[i]);
        cout<< ans <<endl;
        ans = 0;
    }
    return 0;
}