回文串包括奇数长的和偶数长的,一般求的时候都要分情况讨论,这个算法做了个简单的处理把奇偶情况统一了。算法的基本思路是这样的,把原串每个字符中间用一个串中没出现过的字符分隔开来(统一奇偶),用一个数组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数组的性质,新串中以str[i]为中间字符的回文串的长度为p[i]-1,以#为中间字符的就是长度为偶数的,以非#号为中间字符的就是长度为奇数的,那么怎么求p[ ]数组呢。
从左到右计算,也就是计算p[i]时 p[0.....i-1] 都以计算出,并且用一个变量mx记录 max{ k+p[ k ] } (k=0.....i-1),用id记录取最大值时的k, 则 p[ i ]= min( p[2*id - i ], mx - i )
对于第一幅图以i为中间字符的回文串被以id为中间字符的回文串所覆盖,由对称性,p[ i ] = p[ 2*id - i ] 。对于第二幅图没有完全被覆盖,所以对于k>mx的字符,要一个一个匹配,才能确定p [ i ]。
以下是Hiho上通过了的代码:
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#define MAXSIZE 1000002+1002
#define MIN(a,b) ((a)>(b)?(b):(a))
char str[MAXSIZE];
char newstr[2*MAXSIZE];
int p[2*MAXSIZE];
int count_length();
int MIN(int a,int b);
int main()
{
int num,i=1;
freopen("a.txt","r",stdin);
scanf("%d",&num);
p[0]=-1;
while(i<=num)
{
scanf("%s",str);
if(i==num)
printf("%d",count_length());
else
printf("%d\n",count_length());
i++;
}
}
int count_length()
{
int num=0,mx=0,id,i,j;
newstr[0]='$';
for(i=0,j=1;str[i]!='\0';i++)
{
newstr[j++]='#';
newstr[j++]=str[i];
}
newstr[j++]='#';
newstr[j]='\0'; //这句话不能掉,或者每次调用之前用menset将内存全部赋值为0
for(i=1;newstr[i]!='\0';i++)
{
if(mx>i)
p[i]=MIN(p[2*id-i],mx-i);
else
p[i]=1;
for(;newstr[i+p[i]]==newstr[i-p[i]];p[i]++)
;
if(i+p[i]>mx)
{
mx=i+p[i];
id=i;
}
if(p[i]>num)
num=p[i];
}
return num-1;
}