求解最长回文子串 之Manacher算法

回文串包括奇数长的和偶数长的，一般求的时候都要分情况讨论，这个算法做了个简单的处理把奇偶情况统一了。算法的基本思路是这样的，把原串每个字符中间用一个串中没出现过的字符分隔开来(统一奇偶)，用一个数组p[ i ]记录以 str[ i ] 为中间字符的回文串向右能匹配的长度。

原串：w  a   a   b   w   s   w   f   d

新串：   #   w   #   a   #   a   #   b   #   w   #   s    #   w    #     f    #    d     #

             0   1   2   3   4    5   6   7   8   9  10  11 12  13  14   15  16   17   18

p数组：1    2   1   2   3    2   1   2   1   2   1    4   1    2   1     2    1    2    1

可以证明P[i]-1即为原串中str中 最长回文子串

因此问题归结到怎么求出P数组

由p数组的性质，新串中以str[i]为中间字符的回文串的长度为p[i]-1，以#为中间字符的就是长度为偶数的，以非#号为中间字符的就是长度为奇数的，那么怎么求p[ ]数组呢。

从左到右计算，也就是计算p[i]时 p[0.....i-1] 都以计算出，并且用一个变量mx记录 max{ k+p[ k ] } (k=0.....i-1),用id记录取最大值时对于第一幅图以i为中间字符的回文串被以id为中间字符的回文串所覆盖，由对称性，p[ i ] = p[ 2*id - i ] 。对于第二幅图没有完全被覆盖，所以对于k>mx的字符，要一个一个匹配，才能确定p [ i ]。的k， 则 p[ i ]= min( p[2*id - i ], mx - i )

p[ i ]= min( p[2*id - i ], mx - i )这处不是很容易理解：

//记j = 2 * id - i，也就是说 j 是 i 关于 id 的对称点。
if (mx - i > P[j]) 
    P[i] = P[j];
else /* P[j] >= mx - i */
    P[i] = mx - i; // P[i] >= mx - i，取最小值，之后再匹配更新

对于第一幅图以i为中间字符的回文串被以id为中间字符的回文串所覆盖，由对称性，p[ i ] = p[ 2*id - i ] 。

对于第二幅图没有完全被覆盖，所以对于k>mx的字符，要一个一个匹配，才能确定p [ i ]。

所有会有p[ i ]= min( p[2*id - i ], mx - i )

//str 为插入字符后的串，在串首加入个没有出现的字符如‘￥’，防止比较的时候指针小于0
//所以实际str是从下标1开始的
void calp(int n)
{
	int i,mx=0,id;
	for(i=1;i<=n;i++)
	{
		if(mx>i)
			p[i]=min(p[2*di-i],mx-i);
		else
			p[i]=1;
		for(;str[i+p[i]]==str[i-p[i]];p[i]++)  //没完全覆盖
			;
		if(p[i]+i>mx)
		{
			mx=p[i]+i;
			id=i;
		}
	}
}

<span style="color: rgb(51, 51, 51); font-family: Arial; font-size: 18px; line-height: 26px;">因为最多比较字符n次，所以复杂度为O（n）</span>

<pre name="code" class="cpp">#include<vector>
#include<iostream>
using namespace std;

const int N=300010;
int n, p[N];
char s[N], str[N];

#define _min(x, y) ((x)<(y)?(x):(y))

void kp()
{
    int i;
    int mx = 0;
    int id;
    for(i=n; str[i]!=0; i++)
        str[i] = 0; //没有这一句有问题。。就过不了ural1297，比如数据：ababa aba
    for(i=1; i<n; i++)
    {
        if( mx > i )
            p[i] = _min( p[2*id-i], p[id]+id-i );
        else
            p[i] = 1;
        for(; str[i+p[i]] == str[i-p[i]]; p[i]++)
            ;
        if( p[i] + i > mx )
        {
            mx = p[i] + i;
            id = i;
        }
    }
}

void init()
{
 int i, j, k;
 str[0] = '$';
 str[1] = '#';
 for(i=0; i<n; i++)
 {
  str[i*2+2] = s[i];
  str[i*2+3] = '#';
 }
 n = n*2+2;
 s[n] = 0;
}

int main()
{
 int i, ans;
 while(scanf("%s", s)!=EOF)
 {
  n = strlen(s);
  init();
  kp();
  ans = 0;
  for(i=0; i<n; i++)
   if(p[i]>ans)
    ans = p[i];
  printf("%d\n", ans-1);
 }
 return 0;
}

原文 http://blog.163.com/kevinlee_2010/blog/static/169820820201111210740329/