[Luogu3805] 【模板】manacher算法 [Manacher][PAM]

最新推荐文章于 2024-02-21 15:53:43 发布

_er

最新推荐文章于 2024-02-21 15:53:43 发布

阅读量149

点赞数

本文链接：https://blog.youkuaiyun.com/Estia_/article/details/89683492

版权

本文深入探讨了Manacher算法的原理与应用，通过详细解释并提供了一个完整的C++实现案例，展示了如何利用Manacher算法高效地求解最长回文子串问题。文章特别注意到了边界条件的处理，并附带了一次SIGTRAP异常的记录，为读者提供了宝贵的调试经验。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

Link
Luogu - https://www.luogu.org/problemnew/show/P3805

板子。
第一次踩了 SIGTRAP 纪念一下
但是我并没有写 PAM

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<cmath>
#include<ctime>
#include<cctype>
using namespace std;
#define R register
const int MAXN = 22e6 + 5;
int n, Rad[MAXN];
string S;
void Manacher()
{
	for (R int Right = -1, Mid = -1, Tmp = -1, i = 0; i < n; ++i)
	{
		if (Right <= i)
		{
			while ((i + Rad[i] < n - 1) && (i - Rad[i] > 0) && (S[i - Rad[i] - 1] == S[i + Rad[i] + 1])) ++Rad[i];
			Right = i + Rad[i];
			Mid = i;
		}
		else
		{
			Tmp = i + Rad[(Mid << 1) - i - Rad[i]];
			if (Tmp > Right)
			{
				Rad[i] = Right - i;
			}
			else
			{
				Rad[i] = Tmp - i;
				if (Tmp == Right)
				{
					while ((i + Rad[i] < n - 1) && (i - Rad[i] > 0) && (S[i - Rad[i] - 1] == S[i + Rad[i] + 1])) ++Rad[i];
					Right = i + Rad[i];
					Mid = i;
				}
			}
		}
	}
}
int main()
{
	cin >> S;
	n = S.size();
	S.reserve((n<<1)|1);
	for (R int i = n - 1; i >= 0; --i)
	{
		S[i << 1 | 1] = S[i];
		S[i << 1] = '#';
	}
	n <<= 1;
	++n;
	S[n-1] = '#';
	Manacher();
	int Ans = 0;
	for (R int i = 0; i < n; ++i)
	{
		Rad[i] = ((Rad[i] >> 1) << 1) + (i & 1);
		if (Rad[i] > Ans) Ans = Rad[i];
	}
	cout << Ans;
	return 0;
}