最长回文子串（O(n) )Manacher算法

http://www.spoj.com/problems/MSUBSTR/
As we all know Utkarsh is very good at solving number based problems, this time Arpit thinks smartly and gives Utkarsh to solve a problem on Strings. Arpit gives Utkarsh a string and challenges him to find the length of largest substring that have its mirror string same as its original one and number of such substrings. Now Utkarsh is buzy at preparing Avishkar papers so he asks you to help him in doing this task.

Eg. for mirror string : Consider string “lalit” then its mirror string will be “tilal”.

Input :

There are t numbers of test cases (t<=200) followed t lines where each line contains a character string of lower case characters (a-z) of length l (1>=l<=3000).

Output :

There will be two integers per line separated by space indicating the length of largest substring which have its mirror string same and number of such substrings.

输出最长回文子串的长度和个数

Sample test cases :

Input :

3

lalit

abedcdetr

abcde

Output :

3 1

5 1

1 5
数据范围明显要求我们用O(n)来解决，就是Mancher算法。
这个算法有一个很巧妙的地方，它把奇数的回文串和偶数的回文串统一起来考虑了。这一点一直是在做回文串问题中时比较烦的地方。这个算法还有一个很好的地方就是充分利用了字符匹配的特殊性，避免了大量不必要的重复匹配。
算法大致过程是这样。先在每两个相邻字符中间插入一个分隔符，当然这个分隔符要在原串中没有出现过。一般可以用‘#’分隔。这样就非常巧妙的将奇数长度回文串与偶数长度回文串统一起来考虑了（见下面的一个例子，回文串长度全为奇数了），然后用一个辅助数组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值。
if( mx > i)
p[i]=MIN( p[2*id-i], mx-i);

就是当前面比较的最远长度mx>i的时候，P［i］有一个最小值。这个算法的核心思想就在这里.

#include <cstdio>
#include <iostream>
#include <cstring>
const int  maxn = 1e6+10;
using namespace std;
char s[maxn],t[maxn*2];
int p[maxn*2];
int main(){
    int T;
    scanf("%d",&T);
    while(T--){
        scanf("%s",s);
        memset(p,0,sizeof(p));
        memset(t,'@',sizeof(t));
        int n = strlen(s),res = 1;
       // printf("n = %d\n",n);
        t[0] = t[1] = '#';
        for(int i = 0 ; i < n ;i++){
            t[2*(i+1)] = s[i];
            t[2*(i+1)+1] = '#';
        }
        int ans=0,maxx = 0;
        for(int i = 2 ; i <= 2*n+1 ; i++){
            if(maxx+p[maxx] > i)    p[i] = min(p[2*maxx-i],p[maxx]+maxx-i);
            else    p[i] = 1;
            while(t[i+p[i]] == t[i-p[i]])   p[i]++;
            if(p[i]+i > maxx+p[maxx]){
                maxx = i;
            }
            if(ans == p[i] )    res++;
            if(ans < p[i]){
                ans = p[i];
                res = 1;
            }
        }
        printf("%d %d\n",ans-1,res);
    }
}