hdu 3068 最长回文（Manacher算法入门模板题）+解题套路

题目链接：【hdu 3068】

输入有a~z组成的字符串，输出最长回文子串的长度

*Manacher算法模板  ==> 一个专用于求最长回文串的算法

1、转化字符串：abcba  ==>  $#a#b#c#b#a#（前面加$防止数组溢出）

2、s:  $  #  a  #  b  #  c  #  b  #  a  #

     p:     1   2  1  2  1  6  1  2  1  2  1

     p[i]的含义：以s[i]为中心的回文串右扩张的最长长度，p[i] - 1就是转化前的字符串以i为中心的回文子串的长度

3、求p[]，设id为回文子串的中心，mx为边界，则mx = id + p[id]，i与j关于id对称

     (1)id<i<mx

          ==>p[j] < mx-i 

     

         整张图是关于id对称的，所以p[i] = p[i] 

        ==>p[j]>=mx-i

      

          此时p[i] = mx-i

     (2)i>=mx

         这里的p[i]初始值只能是1，通过后续计算慢慢累加

void getp()
{
	memset(p, 0, sizeof(p));
	int id = 0, mx = 0;
	for(int i=1; i<len; i++)
	{
		int j = 2*id-i;
		p[i] = mx>i ? min(p[j], mx-i):1;
		while(s[i+p[i]] == s[i-p[i]]) p[i]++;
		if(i+p[i] > mx)
		{
			mx = i+p[i];
			id = i;
		}
	}
}

所以本题的完整代码就是：

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <string>
using namespace std;
#define inf 310010
char s[inf], str[inf];
int p[inf], len;
void gets()
{
	int li = strlen(str);
	len = 0;
	s[len++] = '&';
	s[len++] = '#';
	for(int i=0; i<li; i++)
	{
		s[len++] = str[i];
		s[len++] = '#';
	}
}
void getp()
{
	gets();
	memset(p, 0, sizeof(p));
	int id = 0, mx = 0;
	for(int i=1; i<len; i++)
	{
		int j = 2*id-i;
		p[i] = mx>i ? min(p[j], mx-i):1;
		while(s[i+p[i]] == s[i-p[i]]) p[i]++;
		if(i+p[i] > mx)
		{
			mx = i+p[i];
			id = i;
		}
	}
}
int main()
{
	while(~scanf("%s", str))
	{
		getp();
		int ans=1;
		for(int i=0; i<len; i++) ans = max(ans, p[i]-1);
		printf("%d\n", ans);
	}
	return 0;
}