Hduoj3068【manacher算法】

/*最长回文
Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 8840    Accepted Submission(s): 3042

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 
*/ 
#include<stdio.h>
#include<string.h> 
#define min(A,B)   (A)<(B)?(A):(B) 
char s[110010], st[220020];
int p[220020]; 
int main()
{
	while(scanf("%s", s) != EOF)
	{
		int i, j, l;
		st[0] = '@';
		for(i = 0; s[i]; i++)//求出新的数组
		{
			st[i*2+1] = '#';
			st[i*2+2] = s[i];
		} 
		st[i*2+1] = '#';
		st[i*2+2] = '\0';
		int mx = 0, id, ans = 0;
		for(j = 1; st[j]; j++)//求每个字符为中心的最长回文串
		{
			if(mx > j)   //若回文串最右端大于当前位置
			p[j] = min(p[2*id - j] , mx-j);//则当前位置有最小半径
			else
			p[j] = 1;                   
			while(st[p[j]+j] == st[j-p[j]])//扩展最小半径两边的字符，若相同则半径长度++
			p[j]++;
			if(p[j] + j > mx)//若新的最右端大于原来的更新中心及最右端
			{
				id = j;
				mx = p[j] + j;
			}
			if(p[j] > ans)//保存最大值
			ans = p[j];
		}
		printf("%d\n", ans-1);
	}
}

题意：题意很清楚，就是求一个字符串的最大长度的回文子串。

思路：这道题用普通的方法都会超时，所以得用manacher算法--附链接http://blog.csdn.net/ggggiqnypgjg/article/details/6645824

体会：注意str【0】一定要设置成不同的字符，不能和分隔符一样，否则有可能数组下标跑到负数出现位置错误。