HDU 3068 最长回文 Manacher算法O(n)回文子串算法

http://acm.hdu.edu.cn/showproblem.php?pid=3068

最长回文

Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)

Total Submission(s): 21482    Accepted Submission(s): 7772

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
  
  

题意：找出字符串中最长的回文串，输出其长度。

思路：Manacher算法O(n)回文子串算法，

http://blog.csdn.net/xingyeyongheng/article/details/9310555这篇博客讲的很清楚，虽然我刚开始也是一脸迷茫，不知道他写的是啥，要想学会这个算法，你就最好自己模拟代码，博客里有栗子，自己手推吧。

http://blog.csdn.net/ggggiqnypgjg/article/details/6645824这篇博客也挺好的

有一点开始的时候一直不明白，这里我要点一下。

例如：串aa（偶串），他的对称中心没有字母，串aaa（奇串），他的对称中心有字母。为了统一处理，所以加上了符号#，这样奇串偶串都有了对称中心的字母。

  原串：    w a a b w s w f d

   i      1 2 3 4 5 6 7 8 9 10111213141516171819
  新串：   # w # a # a # b # w # s # w # f # d #
辅助数组P: 1 2 1 23 2 1 2 1 2 1 4 1 2 1 2 1 2 1

p[id]代表以id为中心，p[id]为半径的回文串（这个回文串包括#）

例如：p[5]=3表示以i=5为中心，有长度为2*p[5]-1回文串。

p数组有一个性质，就是p[i]-1就表示原串中（以i为中心）回文串的长度

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<algorithm>
using namespace std;

const int MAX=110000+10;
char s[MAX*2];
int p[MAX*2];

int main()
{
    while(~scanf("%s",s)){
        int len=strlen(s),id=0,maxlen=0;
        for(int i=len;i>=0;--i){//插入'#'
            s[i+i+2]=s[i];
            s[i+i+1]='#';
        }//插入了len+1个'#',最终的s长度是1~len+len+1即2*len+1,首尾s[0]和s[2*len+2]要插入不同的字符
        s[0]='*';//s[0]='*',s[len+len+2]='\0',防止在while时p[i]越界
        for(int i=2;i<=2*len;++i)
        {
            if(p[id]+id>i) p[i]=min(p[2*id-i],p[id]+id-i);//这是算法避免重复计算的核心（重点模拟找规律）
            else p[i]=1;
            while(s[i-p[i]] == s[i+p[i]])++p[i];
            if(id+p[id]<i+p[i])id=i;
            if(maxlen<p[i])maxlen=p[i];//统计最大值
        }
        cout<<maxlen-1<<endl;
    }
    return 0;
}