【代码超详解】POJ 3974 Palindrome（Manacher 算法 · 模板）

最新推荐文章于 2022-12-24 18:58:07 发布

原创最新推荐文章于 2022-12-24 18:58:07 发布 · 254 阅读

1 ·

CC 4.0 BY-SA版权

文章标签：

#详解

ACM-ICPC 专栏收录该内容

238 篇文章

订阅专栏

博客详细解析了POJ 3974 Palindrome问题，重点介绍了Manacher算法的原理和应用，强调fgets在读取字符串时的注意事项，以确保正确读取包含换行符的输入。

一、题目描述

在这里插入图片描述

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

Manacher 算法
注意：fgets 的第二个参数给出的长度需要至少比字符串的最长长度大 2，因为 fgets 会保存换行符 ‘\n’，并在读取到 ‘\n’ 后写 0 结束字符串。如果只比字符串的最长长度大 1，那么当字符串的长度达到最大时，紧跟其后的 ‘\n’ 将不会被读取。下一次 fgets 时会因为读取到了 ‘\n’ 而直接结束读取，导致 WA。

三、AC 代码

#include<cstdio>
#include<algorithm>
#include<cstring>
#pragma warning(disable:4996)
using namespace std;
template<class _Ty> inline void manacher(const char* 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])++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])++k;
		d2[i] = k; --k;
		if (i + k > r) { l = i - k - 1; r = i + k; }
	}
}
const unsigned lmax = 1000003; char s[lmax]; unsigned d1[lmax], d2[lmax], c, l;
int main() {
	for (;;) {
		fgets(s + 1, lmax - 1, stdin); l = strlen(s + 1);
		if (s[1] == 'E')return 0;
		manacher(s, l, d1, d2);
		printf("Case %u: %u\n", ++c, max(*max_element(d1 + 1, d1 + l + 1) * 2 - 1, *max_element(d2 + 1, d2 + l + 1) * 2));
	}
}