Manacher算法学习笔记

本文深入讲解了Manacher算法,一种高效解决回文问题的算法。通过预处理和核心思想的阐述,介绍了如何在O(n)的时间复杂度内找出字符串中的最长回文子串。文章详细解释了算法的实现步骤,包括预处理过程、核心思想和代码实现。

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

前言

ManacherManacherManacher算法,又称马拉车算法,是解决回文问题的利器,在字符串题中也非常常用。

最重要的是,它是个简单易懂的算法。


让我们从暴力开始

如何求出一个字符串中最长回文串的长度

废话,当然是暴力。

我们可以枚举字符串中的每一个或两个字符作为回文串的中心,然后向外扩展直至两端字符不同,这样就可以在O(n2)O(n^2)O(n2)的时间复杂度内求出答案了。

而且,有些出题人比较懒,造的数据比较水,暴力可能跑得飞快。

但是,如果出题人存心要卡你呢?例如下面这个数据:

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa......aaaaaaaaaaaaaaaaaaaaaaaaaaaaaa(50000个a)

这时候,就要用ManacherManacherManacher算法了。


一个简单的预处理

在暴力的过程中,我们可以发现一个很烦人的问题:回文串长度的奇偶性。这样就需要操作两次,非常麻烦。

怎么办呢?

我们就需要进行一波预处理。

以一个字符串"HLAKNOI"为例,为了解决回文串长度的奇偶性问题,我们可以在该字符串的每两个字符中间插入一个不可能在字符串中出现的字符(随便什么字符都可以,推荐使用’%’,毕竟’%’(膜)可以带来好运),并在首尾各加入一个不同的且同样不可能在字符串中出现的字符(这样一来就不需要判边界了)。

于是,原来的字符串就变成了这样:"!%H%L%A%K%N%O%I%?"。

这样,再去求回文就方便了许多。


一个简单的小性质

我们可以用ppp数组来记录从预处理后每个字符出发的回文半径,则对于一个字符串"abcbad",它的ppp数组如下:

!%a%b%c%b%a%d%?
ppp112121612121211

仔细观察可以发现,每个字母ppp值减111,就是以这个字母为中心所能得到的最长回文串长度

根据这个性质,不难想到,只要求出ppp数组,我们就能快速地求出整个字符串中最长的回文串长度了。


ManacherManacherManacher算法的核心思想

使用了一个简单的预处理,了解了一个简单的小性质,我们就可以开始理解 ManacherManacherManacher算法的核心思想 了。

首先,我们要考虑一个问题:暴力为什么跑得慢?

答案似乎很简单,因为暴力在求解的过程中,对很多位置会重复访问多次。

那么,能不能只访问一次呢?

ManacherManacherManacher算法正是建立于这个想法之上的。


ManacherManacherManacher算法的大致思路

第一步,我们用一个变量MaxMaxMax来记录已经求解过的位置所产生的回文串能到达的最大的右边界,并记录下该回文串的中心位置ididid

现在,假设我们要对一个新位置xxx求解,由于我们是从左到右操作的,因此可以保证x>idx>idx>id,但是,xxxMaxMaxMax的大小关系就难以确定了,因此我们要对其进行分类讨论:

  • x>Maxx>Maxx>Max
    • 对于这种情况,说明以xxx为中心的回文串还没有被求解过,因此我们就对xxx暴力求解。
    • 虽然这种情况的操作时间复杂度较高,但是要注意,这种情况出现次数越多,平均操作时间就越少,大致是成反比的,因此时间复杂度还是可以保证的。
  • x&lt;Maxx&lt;Maxx<Max
    • 对于这种情况,我们可以找到xxx关于ididid的对称点2∗id−x2*id-x2idx,然后对2∗id−x2*id-x2idx为中心的最长回文串的长度再次进行分类讨论:
    • 2∗id−x2*id-x2idx为中心的最长回文串的长度较小,这个回文串在ididid为中心的最长回文串的左半部分范围之内
      • 由于回文串的性质,我们可以确定px≥p2∗id−xp_x≥p_{2*id-x}pxp2idx
      • 因此,我们可以将pxp_xpx赋值为p2∗id−xp_{2*id-x}p2idx,然后继续向两边扩展
    • 2∗id−x2*id-x2idx为中心的最长回文串的长度较小,这个回文串在ididid为中心的最长回文串的左半部分范围之外
      • 这时,我们只能保证在以ididid为中心的最长回文串的右半部分范围之内的以xxx为中心的串是回文串
      • 因此,我们可以将pxp_xpx赋值为pid+id−xp_{id}+id-xpid+idx
    • 其实上面两种情况我们可以直接概括为将pxp_xpx赋值为min(p2∗id−x,pid+id−x)min(p_{2*id-x},p_{id}+id-x)min(p2idx,pid+idx)

ManacherManacherManacher算法的大致思路就是这样,如果还是有地方不懂,可以自己再多想想理解一下。


代码实现
class Class_Manacher//Manacher算法
{
	private:
		int len,p[(N<<1)+5];
		char ns[(N<<1)+5];
		inline void Init(string st)//一个简单的初始化 
		{
			register int i,l=st.length();
			for(ns[i=0]='!',ns[len=1]='%';i<l;++i) ns[++len]=st[i],ns[++len]='%';
			ns[++len]='?';
		}
	public:
		inline int GetAns(string st)//求解
		{
			register int i,ans=0,id,Max=0;
			for(Init(st),i=1;i<len;++i)
			{
				p[i]=i<=Max?min(p[(id<<1)-i],p[id]+id-i):1;//分情况讨论,求出p[i]的初值
				while(!(ns[i-p[i]]^ns[i+p[i]])) ++p[i];//向外扩展
				if(i+p[i]>Max) Max=i+p[id=i];//更新所能到达的最大右边界
			}
			for(i=1;i<len;++i) ans=max(ans,p[i]);//求出最大的p[i]
			return ans-1;//最大的p[i]减1才是答案
		}
}Manacher;

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值