Manacher回文串算法 hdu3068 hdu4513

学习资料：http://blog.csdn.net/ggggiqnypgjg/article/details/6645824

时间复杂度O(N)

模板：

void Manacher(int n)
{
    int i;
    int mx=0;//记录前面回文串最长影响到的地方。不一定是最长回文串造成的。
    int id;//最长影响串的ID;
    p[0]=0;
    for(i=1;i<n;i++)
    {
        p[i]=1;//至少是1
        if(mx>i)//i受到影响即，id+p[id]-1>=i;
        {
            p[i]=p[2*id-i];//2*id-i是i关于id的对称点相当于是id-(i-id);
            if(mx-i<p[i])p[i]=mx-i;
            //由于对称点的回文半径可能超过mx-i,因为mx后面的还没有配过
            //所以要取小的。等等继续配
        }
  
        //向两端配匹
        while(str[i-p[i]]==str[i+p[i]])p[i]++;
  
        if(i+p[i]>mx)
        {
            mx=i+p[i];
            id=i;
        }
    }
}

hdu3068：求一个字符串的最长回文子串,就是一个满足本身是回文的最长的子串

#include<cstdio>
#include<cstring>
const int M = 110010;
char s[M],str[M*2];
int p[M*2];
int n;
int max(int a,int b){
    return a>b?a:b;
}
int min(int a,int b){
    return a<b?a:b;
}
void init(){
    int i,j=0;
    n=strlen(s);
    str[j++]='$';str[j++]='#';
    for(i=0;i<n;i++){
        str[j++]=s[i];
        str[j++]='#';
    }
    str[j]='\0';
    n=j;
}
int solve(){
    int i,id,mx=0;
    for(i=1;i<n;i++){
        if(mx>i) p[i]=min(mx-i,p[2*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;
        }
    }
    int res=0;
    for(i=1;i<n;i++)
        res=max(res,p[i]);
    return res;
}
int main(){
    while(scanf("%s",&s)!=EOF){
        init();
        printf("%d\n",solve()-1);
    }
    return 0;
}

hdu4513：

#include<cstdio>
#include<cstring>
const int M = 100010;
int s[M],str[M*2];
int p[M*2];
int T,n;
int max(int a,int b){
    return a>b?a:b;
}
int min(int a,int b){
    return a<b?a:b;
}
void init(){
    int i,j=0;
    str[j++]=-1;str[j++]=0;
    for(i=0;i<n;i++){
        str[j++]=s[i];
        str[j++]=0;
    }
    n=j;
}
int solve(){
    int i,id,mx=0;
    for(i=1;i<n;i++){
        if(mx>i) p[i]=min(mx-i,p[2*id-i]);
        else p[i]=1;
    while(str[i+p[i]]==str[i-p[i]]){
        if(str[i+p[i]]<10) p[i]++;
        else {
            int a=i-p[i];
                    int b=i+p[i];
  
                    if(a+2<=b-2&&str[a]<=str[a+2]&&str[b-2]>=str[b]){
                            p[i]++;
                            continue;
                    }
                    else if(a==b){
                            p[i]++;
                            continue;
                       }
            else if(a+2>b-2){
                            p[i]++;
                            continue;
                    }
                    break;
                }
    }        
        if(p[i]+i>mx){
            mx=p[i]+i;
            id=i;
        }
    }
    int res=0;
    for(i=2;i<n;i++)
        res=max(res,p[i]);
    return res;
}
int main(){
    int i;
    scanf("%d",&T);
    while(T--){
        scanf("%d",&n);
        for(i=0;i<n;i++){
            scanf("%d",&s[i]);
        }
        init();
        printf("%d\n",solve()-1);
    }
    return 0;
}