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


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

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

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

[align=left]Sample Input[/align]

aaaa

abab

[align=left]Sample Output[/align]

4
3

[align=left]Source[/align]
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;
}

View Code