Manacher算法

此篇文章给出了最长回文串的普通解法：

https://blog.csdn.net/qq_45987579/article/details/121324132?spm=1001.2014.3001.5501
但是这种解法计算每个遍历点时并没有充分利用之前遍历点计算出的结论，时间复杂度为O（n2），本文算法充分利用回文字符串的性质，降低算法时间复杂度，该算法就是马拉车算法。

首先给出一段字符串分析：

“abctftghjhgtftclu”，这个字符串首先可以拆分成3段：
第一段：“ab”
第二段：“ctftghjhgtftc”
第三段：“lu”
第一段和第三段的最长回文字符串长度皆为1，第二段字符串是一串长度为11的回文字符串，首先我们按照通用的方法进行遍历，从’a’到第一个’h’得到的最大回文子串长度为3（abctftghjhgtftclu)，遍历到j时,更新最长回文子串长度为11（abctftghjhgtftclu），继续往后遍历时，会发现"ctftghjhgtftc"的后半部分和前半部部分是对称的，例如"abctftghjhgtftclu"后面的f为中心计算回文子串时，可以参考"abctftghjhgtftclu"前面的f计算出来的结果，同时还需要考虑对称回文串是否在边界的情况。
为了方便算法的实现，如果每个回文子串长度都是奇数例如：“aba”，而不是偶数例如："baab"形式的，将会给我们的计算带来很大的遍历，所以在算法初始时需要做一下特殊操作：在源字符串首尾以及每两个字符之间插入一个特殊标记字符，这样所有的回文子串长度就都是奇数了：#a#b#c#t#f#t#g#h#j#h#g#t#f#t#c#l#u#

马拉车解法实现：

int FindPalindrome(char *s, int left, int right)
{
	while (*(s + left) == *(s + right))
	{
		if ((left == 0) || (right == strlen(s) - 1))
		{
			return right - left + 1;
		}
		--left;
		++right;
	}

	return right - left - 1;
}

char * longestPalindrome(char * s)
{
	int i = 0;
	char *my_s = calloc(1, strlen(s) * 2 + 2);
	int *p = calloc(1, (strlen(s) * 2 + 1) * 4);
	int center = 0, right_max = 0, max_p = 0, max_i = 0;
	
	for (i = 0; i < strlen(s); i++)
	{
		*(my_s + i * 2) = '#';
		*(my_s + i * 2 + 1) = *(s + i);  
	}
	
	*(my_s + i * 2) = '#';

	 for (i = 0; i < strlen(my_s) - 1; i++)
	 {
	 
		if (i >= right_max)
		{
			/* 中心扩散 */
			*(p + i) = (FindPalindrome(my_s, i, i) - 1) / 2;
			right_max = *(p + i) + i;
			center = i;
		}
		else
		{
			if (*(p + (2 * center - i)) < right_max - i)
			{
				/* 使用对称点的结果 */
				*(p + i) = *(p + (2 * center - i));
				continue;
			}
			else if (*(p + (2 * center - i)) > right_max - i)
			{
				*(p + i) = right_max - i;
			}
			else
			{
				/* 中心扩散 */
				*(p + i) = (FindPalindrome(my_s, 2 * i - right_max, right_max) - 1) / 2;
				right_max = *(p + i) + i;
				center = i;
			}
		}

		max_i = *(p + i) > max_p ? i : max_i;
		max_p = *(p + i) > max_p ? *(p + i) : max_p;
	}

	char *ret = calloc(max_p + 1, 1);
	memcpy(ret, s + (max_i - max_p) / 2, max_p);

	free (my_s);
	free (p);
	return ret;
}