hdu3068 最长回文--Manacher算法

原题链接： http://acm.hdu.edu.cn/showproblem.php?pid=3068

一：原题内容

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
  
  

二：分析理解

参见： http://blog.csdn.net/laojiu_/article/details/50931672

三：AC代码

#include<iostream>
#include<string.h>
#include<algorithm>
#include<utility>
using namespace std;

char s[110010];
char s_new[220020];
int p[220020];

int Init()
{
	int len = strlen(s);
	s_new[0] = '$';
	s_new[1] = '#';
	int j = 2;

	for (int i = 0; i < len; i++)
	{
		s_new[j++] = s[i];
		s_new[j++] = '#';
	}

	s_new[j] = '\0';//别忘了哦

	return j;//返回s_new的长度
}

int Manacher()
{
	int len = Init();//取得新字符串长度并完成向s_new的转换
	int maxLen = -1;//最长回文长度
	
	int id;
	int mx = 0;

	for (int i = 1; i < len; i++)
	{
		if (i < mx)
			p[i] = min(p[2 * id - i], mx - i);//需搞清楚上面那张图含义,mx,2 * id - i的含义
		else
			p[i] = 1;

		while (s_new[i - p[i]] == s_new[i + p[i]])//不需边界判断，因为左有'$',右有'\0'
			p[i]++;


		if (mx < i + p[i])//我们每走一步i，都要和mx比较，我们希望mx尽可能的远，这样才能更有机会执行if (i < mx) 这句代码，从而提高效率
		{
			id = i;
			mx = i + p[i];
		}

		maxLen = max(maxLen, p[i] - 1);
	}

	return maxLen;
}

int main()
{
	while (scanf("%s", s)==1)
		printf("%d\n", Manacher());
	
	return 0;
}