首页 > 代码库 > 最长回文

最长回文

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


最长回文