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算法

 

代码：

 1 #include<stdio.h>

 2 #include<string.h>

 3 #include<stdlib.h>

 4 #define maxn 110020

 5 char str[maxn*2];

 6 int rad[maxn*2];

 7 int Min(int a ,int b){

 8   return a<b?a:b;

 9 }

10 void init(char s[],int len){

11     int i,j=1;

12   for(i=(len<<1)+1,s[i+1]='\0';i>0;i--){

13         if(i&1)  s[i]='#';

14         else {

15        s[i]=s[len-j];

16        j++;

17         }

18   }

19   s[i]='$';

20 }

21 int manacher(int len){

22   memset(rad,0,sizeof(int)*(len+1));

23   int id,mx=0,i,ans=0;

24   for(i=1;i<len;i++){

25       if(mx>i) rad[i]=Min(rad[id*2-i],mx-i);

26     else rad[i]=1;

27     while(str[i+rad[i]]==str[i-rad[i]])rad[i]++;

28     if(i+rad[i]>mx){

29          mx=i+rad[i];

30          id=i;

31     }

32     if(ans<rad[i])ans=rad[i];

33   }

34   return ans-1;

35 }

36 int main(){

37    int len;

38   while(scanf("%s",str)!=EOF){

39    len=strlen(str);

40    init(str,len);

41    printf("%d\n",manacher(len*2+1));

42  }

43  return 0;

44 }

View Code