hdu----（3068)最长回文（manacher）

最长回文

Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 8253    Accepted Submission(s): 2825

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

 

 

Source

2009 Multi-University Training Contest 16 - Host by NIT

 求最长回文子串.....manacher算法

 

代码：

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#define maxn 110020
char str[maxn*2];
int rad[maxn*2];
int Min(int a ,int b){
  return a<b?a:b;
}
void init(char s[],int len){
    int i,j=1;
  for(i=(len<<1)+1,s[i+1]='\0';i>0;i--){
        if(i&1)  s[i]='#';
        else {
       s[i]=s[len-j];
       j++;
        }
  }
  s[i]='$';
}
int manacher(int len){
  memset(rad,0,sizeof(int)*(len+1));
  int id,mx=0,i,ans=0;
  for(i=1;i<len;i++){
      if(mx>i) rad[i]=Min(rad[id*2-i],mx-i);
    else rad[i]=1;
    while(str[i+rad[i]]==str[i-rad[i]])rad[i]++;
    if(i+rad[i]>mx){
         mx=i+rad[i];
         id=i;
    }
    if(ans<rad[i])ans=rad[i];
  }
  return ans-1;
}
int main(){
   int len;
  while(scanf("%s",str)!=EOF){
   len=strlen(str);
   init(str,len);
   printf("%d\n",manacher(len*2+1));
 }
 return 0;
}