【代码超详解】HDU 4513 吉哥系列故事——完美队形II（Manacher 算法）

一、题目描述

二、算法分析说明与代码编写指导

把 Manacher 算法 改一下，向两端查找更长的回文子串时，必须满足两头的不大于中间的才可以被选入更长的回文子串中。

三、AC 代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#pragma warning(disable:4996)
using namespace std;
const unsigned nmax = 100001;
unsigned t, n, h[nmax], d1[nmax], d2[nmax];
template<class _Ty> inline void manacher(const _Ty* const s, const _Ty& L, _Ty* const d1, _Ty* const d2) {
	_Ty i, k, l, r;
	for (i = 1, l = 1, r = 0; i <= L; ++i) {
		k = (i > r) ? 1 : min(d1[l + r - i], r - i);
		while (i - k > 0 && i + k <= L && s[i - k] == s[i + k] && s[i - k] <= s[i - k + 1])++k;
		d1[i] = k; --k;
		if (i + k > r) { l = i - k; r = i + k; }
	}
	for (i = 1, l = 1, r = 0; i <= L; ++i) {
		k = (i > r) ? 0 : min(d2[l + r - i + 1], r - i + 1);
		while (i - k - 1 > 0 && i + k <= L && s[i - k - 1] == s[i + k] && s[i - k - 1] <= s[i - k])++k;
		d2[i] = k; --k;
		if (i + k > r) { l = i - k - 1; r = i + k; }
	}
}
int main() {
	scanf("%u", &t); ++t;
	while (--t) {
		scanf("%u", &n); for (unsigned i = 1; i <= n; ++i)scanf("%u", h + i);
		manacher(h, n, d1, d2);
		printf("%u\n", max(*max_element(d1 + 1, d1 + 1 + n) * 2 - 1, *max_element(d2 + 1, d2 + 1 + n) * 2));
	}
	return 0;
}