【最长回文子串】HDU3068最长回文【Manacher算法】

一张图领悟Manacher算法,计算字符串最长回文子串

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3068

Problem Description

给出一个只由小写英文字符a,b,c...y,z组成的字符串S,求S中最长回文串的长度.
回文就是正反读都是一样的字符串,如aba, abba等

 

Input

输入有多组case,不超过120组,每组输入为一行小写英文字符a,b,c...y,z组成的字符串S
两组case之间由空行隔开(该空行不用处理)
字符串长度len <= 110000

 

Output

每一行一个整数x,对应一组case,表示该组case的字符串中所包含的最长回文长度.

 

Sample Input

   
   
   
   

    
    
    
    aaaa

abab
   
   
   
   

 

Sample Output

   
   
   
   

    
    
    
    4
3
   
   
   
   
   
   
   
   


   
   
   
   

代码：

#include<iostream>
#include<cstring>
using namespace std;
string s;
int p[110050*2];    //  这里需要注意,数组的长度大于字符串的长度的两倍;
//  将字符串转换成奇字符串
void init()
{
    memset(p,0,sizeof(p));
    string s1="$";
    for(int i=0;i<s.size();i++){
        s1+='#';
        s1+=s[i];
    }
    s1+='#';
    s=s1;
}
//  计算p[id]表示id位置的回文半径
void Manacher()
{
    int id=0;       //  id位置的右边界最靠右(在已知范围内);
    int mx=0;       //  id的右边界位置;
    for(int i=1;i<s.size();i++){
        if(mx>i)
            p[i]=min(p[2*id-i],mx-i);
        else
            p[i]=1;
        for(;s[i-p[i]]==s[i+p[i]];p[i]++);
        if(i+p[i]>mx){
            mx=i+p[i];
            id=i;
        }
    }
}
int main()
{
    cin.sync_with_stdio(false);
    while(cin>>s){
        init();
        Manacher();
        int ans=0;
        for(int i=1;i<s.size();i++){
            if(ans<p[i])
                ans=p[i];
        }
        cout<<ans-1<<endl;
    }
    return 0;
}