首页 > 代码库 > 最长回文
最长回文
回文串包括奇数长的和偶数长的,一般求的时候都要分情况讨论,这个算法做了个简单的处理把奇偶情况统一了。算法的基本思路是这样的,把原串每个字符中间用一个串中没出现过的字符分隔开来(统一奇偶),用一个数组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; }
最长回文