求最长回文字串

1.暴力求解

string findLongestPalindrome(string &s)  
{  
    int length=s.size();//字符串长度  
    int maxlength=0;//最长回文字符串长度  
    int start;//最长回文字符串起始地址  
    for(int i=0;i<length;i++)//起始地址  
        for(int j=i+1;j<length;j++)//结束地址  
        {  
            int tmp1,tmp2;  
            for(tmp1=i,tmp2=j;tmp1<tmp2;tmp1++,tmp2--)//判断是不是回文  
            {  
                if(s.at(tmp1)!=s.at(tmp2))  
                    break;  
            }  
            if(tmp1>=tmp2&&j-i>maxlength)  
            {  
                maxlength=j-i+1;  
                start=i;  
            }  
        }  
        if(maxlength>0)  
            return s.substr(start,maxlength);//求子串  
        return NULL;  
}  

2.动态规划

回文字符串的子串也是回文，比如P[i,j]（表示以i开始以j结束的子串）是回文字符串，那么P[i+1,j-1]也是回文字符串。这样最长回文子串就能分解成一系列子问题了

P[i,j]=0表示子串[i,j]不是回文串。P[i,j]=1表示子串[i,j]是回文串。

string findLongestPalindrome(string &s)  
{  
    const int length=s.size();  
    int maxlength=0;  
    int start;  
    bool P[50][50]={false};  
    for(int i=0;i<length;i++)//初始化准备  
    {  
        P[i][i]=true;  
        if(i<length-1&&s.at(i)==s.at(i+1))  
        {  
            P[i][i+1]=true;  
            start=i;  
            maxlength=2;  
        }  
    }  
    for(int len=3;len<length;len++)//子串长度  
        for(int i=0;i<=length-len;i++)//子串起始地址  
        {  
            int j=i+len-1;//子串结束地址  
            if(P[i+1][j-1]&&s.at(i)==s.at(j))  
            {  
                P[i][j]=true;  
                maxlength=len;  
                start=i;  
            }  
        }  
    if(maxlength>=2)  
        return s.substr(start,maxlength);  
    return NULL;  
}  

3.中心扩展

中心扩展就是把给定的字符串的每一个字母当做中心，向两边扩展，这样来找最长的子回文串。算法复杂度为O(N^2)。

但是要考虑两种情况：

1、像aba，这样长度为奇数。

2、想abba，这样长度为偶数。

string findLongestPalindrome(string &s)  
{  
    const int length=s.size();  
    int maxlength=0;  
    int start;  
  
    for(int i=0;i<length;i++)//长度为奇数  
    {  
        int j=i-1,k=i+1;  
        while(j>=0&&k<length&&s.at(j)==s.at(k))  
        {  
            if(k-j+1>maxlength)  
            {  
                maxlength=k-j+1;  
                start=j;  
            }  
            j--;  
            k++;  
        }  
    }  
  
    for(int i=0;i<length;i++)//长度为偶数  
    {  
        int j=i,k=i+1;  
        while(j>=0&&k<length&&s.at(j)==s.at(k))  
        {  
            if(k-j+1>maxlength)  
            {  
                maxlength=k-j+1;  
                start=j;  
            }  
            j--;  
            k++;  
        }  
    }  
    if(maxlength>0)  
        return s.substr(start,maxlength);  
    return NULL;  
}  

4.Manachar

#include<stdio.h>
#include<string.h>
#include<iostream>
using namespace std;
#define N 3000000
char str[N],ch[N*2];
int p[N*2];
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%s",str);
        int len=strlen(str);
        len=len*2+2;
        ch[0]='#';
        ch[1]='#';
        for(int i=0;i<len;i++)
        {
            ch[i*2+2]=str[i];
            ch[i*2+3]='#';
        }
        ch[len+1]='#';
        int ans=0,mx=0,id=0;
        for(int i=1;i<len;i++)
        {
            if(mx>i)
            {
                p[i]=min(p[2*id-i],mx-i);
            }else {
            p[i]=1;
            }
            for( ;ch[i+p[i]]==ch[i-p[i]];p[i]++)
            {
                if(mx<i+p[i])
                {
                    mx=i+p[i];
                    id=i;
                }
            }
        }
        for(int i=1;i<len;i++)
            ans=max(ans,p[i]-1);
        printf("%d\n",ans);
    }
}