回文子串——Manacher马拉车算法

2024-02-04 12:08

本文主要是介绍回文子串——Manacher马拉车算法,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

我们有些时候会遇到求最大回文字符串的问题。对于回文字符串,我们首先想到的就是N^2的暴力,对于每一个点都向左或者向右延伸。但这样做显然不是特别优秀。而今天介绍的马拉车算法的最大贡献就是将时间复杂度提升到了线性,这是非常了不起的。
回文字符串有两种,一种长度为奇数,如aba。另一种长度为偶数,如abba。对于奇数,只要询问中间每一个点即可,而偶数不存在中间点,所以我们得创造一个中间点。我们将每两个字符间插入一个字符#。即可完成。例如aba变为#a#b#a#,而abba变为#a#b#b#a#。
可见奇数长度回文子串的中间字符不为’#’,而偶数长度子串的字符为#。
我们定义p[i],表示从i开始向右走和向左走字符一样的字符长度。这么讲大家可能有点懵,所以我们演示一下。
对于字符串abbababaaab,加完#后就是#a#b#b#a#b#a#b#a#a#a#b#, 对于p[11](也就是第3个a),其值为6,因为我们扩展出的子串为a#b#a#(包括它自己本身)。
那么问题来了,怎么求数组p呢?在求解之前,我们先规定id为最大回文子串中心的位置,mx是回文串能延伸到的最右端的位置。我先把公式搁这儿:
p[i]=mx>i?min(p[id*2-i],mx-i):1;
为什么是这样的呢?我们可以先进行讨论,如果i>mx,也就是i在我们目前已知的回文串之外,那我们就管不了它了,先让它维护它自己(也就是1),到之后再去暴力枚举,更新mx和id。
而如果i< mx,我们又要进行讨论,若i关于id的对应点j的值p[j]< mx-i,说明i的p值也一定小于mx-i,应为对于子串jid-1和id+1i,它们两两相同。

这里写图片描述

而如果p[j]>mx-i,说名p[i]也可以大于mx-i,所以我们先将p[i]赋值为mx-i,然后再去暴力枚举更新p[i],并更新mx,id。
这样一直更新就是我们的马拉车算法。那么这个算式也就十分容易理解了。通过平均数定力x=(x1+x2)/2得,2*id-i也就是对应点j,所以其实这个算式的意思是
if(mx>i){if(mx-i>p[j]) p[i]=p[j];else p[i]=mx-i;
}
else p[i]=1;

传送门:马拉车算法的练习题

这道题是一道简单的马拉车算法题,因为只需要求技术回文子串,所以不需要添加#(我的程序还是加了),我们只需要用前缀和来维护序列的长度。由于若一个子串长度为l,则这个串里一定至少有一个长为l-2,l-4,l-6……的串(只需要去掉头尾,所以一直-2),所以对于长度为l的串的总数tot[l],只要加上tot[l+2],即可维护总共的长度为l的串的值。
这题还得写快速幂,因为数据较大,直接for来乘可能会超时
#include<bits/stdc++.h>
#define ll long long
#define md 19930726
using namespace std;
ll n,m,k,p[5000005],mx,id,tot[1000005],ans=1;
string s,t;
ll read(){char c;ll x;while(c=getchar(),c<'0'||c>'9');x=c-'0';while(c=getchar(),c>='0'&&c<='9') x=x*10+c-'0';return x;
}
void manacher(){t="@#";   //在串头随意添加一个字符,如@,防止其一直匹配空串for(ll i=0;i<n;i++){t.push_back(s[i]);t.push_back('#');}for(ll i=1;i<t.size();i++){p[i]=mx>i?min(p[id*2-i],mx-i):1;while(t[i+p[i]]==t[i-p[i]]) p[i]++;if(mx<i+p[i]){mx=i+p[i];id=i;}}
}
ll pows(ll a,ll b){ll base=1;while(b){if(b&1) base=base*a%md;a=a*a%md;b/=2;}return base;
}
int main()
{n=read();k=read();cin>>s;manacher();for(ll i=1;i<t.size();i++){if(t[i]=='#') continue;tot[p[i]-1]++;}for(int i=n;i>=1;i--)if(i%2==1) tot[i]+=tot[i+2];for(ll i=n;i>=1;i--){if(!tot[i]) continue;if(k-tot[i]>=0){k-=tot[i];ans=(ans*pows(i,tot[i]))%md;}else{ans=(ans*pows(i,k))%md;break;}}printf("%lld",ans%md);return 0;
}

这篇关于回文子串——Manacher马拉车算法的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/677438

相关文章

SpringBoot实现MD5加盐算法的示例代码

《SpringBoot实现MD5加盐算法的示例代码》加盐算法是一种用于增强密码安全性的技术,本文主要介绍了SpringBoot实现MD5加盐算法的示例代码,文中通过示例代码介绍的非常详细,对大家的学习... 目录一、什么是加盐算法二、如何实现加盐算法2.1 加盐算法代码实现2.2 注册页面中进行密码加盐2.

Java时间轮调度算法的代码实现

《Java时间轮调度算法的代码实现》时间轮是一种高效的定时调度算法,主要用于管理延时任务或周期性任务,它通过一个环形数组(时间轮)和指针来实现,将大量定时任务分摊到固定的时间槽中,极大地降低了时间复杂... 目录1、简述2、时间轮的原理3. 时间轮的实现步骤3.1 定义时间槽3.2 定义时间轮3.3 使用时

如何通过Golang的container/list实现LRU缓存算法

《如何通过Golang的container/list实现LRU缓存算法》文章介绍了Go语言中container/list包实现的双向链表,并探讨了如何使用链表实现LRU缓存,LRU缓存通过维护一个双向... 目录力扣:146. LRU 缓存主要结构 List 和 Element常用方法1. 初始化链表2.

C++实现回文串判断的两种高效方法

《C++实现回文串判断的两种高效方法》文章介绍了两种判断回文串的方法:解法一通过创建新字符串来处理,解法二在原字符串上直接筛选判断,两种方法都使用了双指针法,文中通过代码示例讲解的非常详细,需要的朋友... 目录一、问题描述示例二、解法一:将字母数字连接到新的 string思路代码实现代码解释复杂度分析三、

golang字符串匹配算法解读

《golang字符串匹配算法解读》文章介绍了字符串匹配算法的原理,特别是Knuth-Morris-Pratt(KMP)算法,该算法通过构建模式串的前缀表来减少匹配时的不必要的字符比较,从而提高效率,在... 目录简介KMP实现代码总结简介字符串匹配算法主要用于在一个较长的文本串中查找一个较短的字符串(称为

通俗易懂的Java常见限流算法具体实现

《通俗易懂的Java常见限流算法具体实现》:本文主要介绍Java常见限流算法具体实现的相关资料,包括漏桶算法、令牌桶算法、Nginx限流和Redis+Lua限流的实现原理和具体步骤,并比较了它们的... 目录一、漏桶算法1.漏桶算法的思想和原理2.具体实现二、令牌桶算法1.令牌桶算法流程:2.具体实现2.1

Python中的随机森林算法与实战

《Python中的随机森林算法与实战》本文详细介绍了随机森林算法,包括其原理、实现步骤、分类和回归案例,并讨论了其优点和缺点,通过面向对象编程实现了一个简单的随机森林模型,并应用于鸢尾花分类和波士顿房... 目录1、随机森林算法概述2、随机森林的原理3、实现步骤4、分类案例:使用随机森林预测鸢尾花品种4.1

不懂推荐算法也能设计推荐系统

本文以商业化应用推荐为例,告诉我们不懂推荐算法的产品,也能从产品侧出发, 设计出一款不错的推荐系统。 相信很多新手产品,看到算法二字,多是懵圈的。 什么排序算法、最短路径等都是相对传统的算法(注:传统是指科班出身的产品都会接触过)。但对于推荐算法,多数产品对着网上搜到的资源,都会无从下手。特别当某些推荐算法 和 “AI”扯上关系后,更是加大了理解的难度。 但,不了解推荐算法,就无法做推荐系

康拓展开(hash算法中会用到)

康拓展开是一个全排列到一个自然数的双射(也就是某个全排列与某个自然数一一对应) 公式: X=a[n]*(n-1)!+a[n-1]*(n-2)!+...+a[i]*(i-1)!+...+a[1]*0! 其中,a[i]为整数,并且0<=a[i]<i,1<=i<=n。(a[i]在不同应用中的含义不同); 典型应用: 计算当前排列在所有由小到大全排列中的顺序,也就是说求当前排列是第

csu 1446 Problem J Modified LCS (扩展欧几里得算法的简单应用)

这是一道扩展欧几里得算法的简单应用题,这题是在湖南多校训练赛中队友ac的一道题,在比赛之后请教了队友,然后自己把它a掉 这也是自己独自做扩展欧几里得算法的题目 题意:把题意转变下就变成了:求d1*x - d2*y = f2 - f1的解,很明显用exgcd来解 下面介绍一下exgcd的一些知识点:求ax + by = c的解 一、首先求ax + by = gcd(a,b)的解 这个