manacher算法 求最长回文串

本题链接:用户登录

题目:

样例:

输入
aa1ABA1b

输出
5

思路:

        根据题目意思,求出最长回文串,我们可以用模板 manacher 算法 求最长回文串。

        manacher算法 求最长回文串 核心有两个步骤。

        一、将字符串转化为  奇数 长度的字符串,方便我们求回文串。

        二、利用回文串对称性的特点,推导最长回文串长度。

        

        根据我们的核心开始一步一步来。

第一步,转化为奇数长度的字符串

        我们把每个字符作为回文串的中心去处理。

        我们只要在每个字符中间加入一个该字符串没有的符号,并且在头尾加入没出现过的不同的符号就行了(如下面这个字符串)。

        根据给出的字符串 aa1ABA1b  ,我们将每个字符分隔起来,假设 用 ‘#’ 作为分隔符,头尾分别用 '$' 和 '^' 作为边界分隔符。

        得到的转化字符串为 :     $a#a#1#A#B#A#1#b^

        这样我们就可以得到长度一定为奇数的字符串了。

        新字符串长度 = 原长度 * 2 + 3;(’#’ 个数始终比原字符个数多 1)。

第二步,利用回文串对称性推导最长回文串

        这里怎么利用回文串对称特性呢?

        我们只要维护一个最长 回文串长度 来一步步推导就可以了。

        我们看一下上面这个字符串   

                新字符串:      $a#a#1#A#B#A#1#b^   

                 原字符串:      aa1ABA1b

         假设我们找到了对称 的    字符串 

         1ABA1   这个长度是   5 

        我们再从新字符串看一下这个长度     #1#A#B#A#1#     这个长度是  11

        这里有一个规律,我们截取一下最右边的一半回文串:    B#A#1#      这个长度是  6

        所以我们可以得到     新回文字符串的一半长度 - 1 == 回文字符串长度

       这就是我们利用回文串对称性推导最长回文串。

        所以我们需要两个变量进行推导   一个是字符串中心点  mid   和一个 最右边的回文串长度 mr

        这两个变量就是我们需要   维护 的一个 最长回文串 的变量

        然后我们再用一个 p  数组进行记录我们遍历整个字符串的过程中 最长回文串长度变化即可。

       

        有人会疑问,如果不是对称的情况呢?就利用不到对称性了。

        这里我们就 暴力向外扩张 求对称性 p 即可。

        这里又有个疑问,如果对称的字符串不在我们的 维护的最长回文串里面呢?

        所以这里有两种情况。

        一种是: p[ i ] 关于 最右回文串的对称字符串在 我们维护的 最长回文串 ml 、 mr 范围内

                         那么 ans[ i ] = ans[ j ] (j 为对称子串的中心)。

                        其中 j  等于 mid * 2 -  i 毕竟 (i + j) / 2 == mid   即  p[ i ] = p[ mid * 2 - i ]

        一种是:不在范围内,因为(ml - mr) 是最右的回文串的边界,所以 p[ i ] 等于其对称子串在 ml-mr 中的最大半径。(因为其对称子串范围超过了 ml-mr,且 ml-mr 无法向右扩张,即 s[sr + 1] != s[sl - 1]。 那么此时 p[ i ] 就等于其对称子串在sl-sr范围内的半径   即    p[ i ] ==  mr - i 

  

代码详解如下:

        

#include <iostream>
#include <cstring>
#include <vector>
#define endl '\n'
#define int long long
#define IOS ios::sync_with_stdio(false),cin.tie(0),cout.tie(0)
using namespace std;
const int N = 2e7 + 10;
string a,b;
int n;// 初始化转化字符串为 奇数 长度
inline void Init()
{b = "$#"; // 添加边界分隔符// 更新字符串for(int i = 0;i < n;++i){b += a[i];b += '#';}b += "^"; // 添加边界分隔符n = b.size(); // 更新字符串新长度
}// 马拉车算法求最长回文字符串长度
inline int Manacher()
{Init(); // 初始化vector<int>p(N,0); // 记录最长回文串int mr = 0,  mid; // 定义 mid 以及  mr 暂时还没发现最长回文串,所以 mr  == 0, mid  未知for(int i = 1;i < n;++i){// 如果维护的最长回文串超出了 当前 i 字符串的 边界if(mr > i) p[i] = min(p[(mid << 1) - i],mr - i);  // 那么取相应的 对称情况的 有效的回文串长度else p[i] = 1;  // 如果没超出, 当前最长回文串长度就是 一个字符  长度就是 1while(b[i - p[i]] == b[i + p[i]]) ++p[i]; // 暴力扩张对称的回文串长度,并记录好// 如果 扩张后 存在当前回文串长度比之前维护的回文串长度还长// 那么更新最长回文串   mr   以及  mid if(i + p[i] > mr){mr = i + p[i];mid = i;}}int ans = 0;// 遍历记录数组 p 查找我们推导过程中出现的    最右最长回文串长度// 原字符串最长回文串长度  =   新字符串最右最长回文串长度 - 1    所以  p[i] - 1for(int i = 0;i < n;++i) ans = max(ans,p[i] - 1); return ans;
}signed main()
{IOS;cin >> a;n = a.size();int ans = Manacher();cout << ans << endl;return 0;
}

最后提交:

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/641752.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

SpringBoot分片上传、断点续传、大文件极速秒传功能

前言 文件上传是一个老生常谈的话题了&#xff0c;在文件相对比较小的情况下&#xff0c;可以直接把文件转化为字节流上传到服务器&#xff0c;但在文件比较大的情况下&#xff0c;用普通的方式进行上传&#xff0c;这可不是一个好的办法&#xff0c;毕竟很少有人会忍受&#…

华为机试真题实战应用【算法代码篇】-图谱新词挖掘1(附python、C++和JAVA代码实现)

目录 题目描述 代码实现 Java python C++ 代码2 题目描述 小华负责公司知识图谱产品,现在要通过新词挖掘完善知识图谱新词挖掘: 给出一个待挖掘问题内容字符串Content和一人词的字符串word,找到content中所有word的新词。新词: 使用词word的字符排列形成的字符串。 请帮…

Hive 数仓及数仓设计方案

数仓(Data Warehouse) 数据仓库存在的意义在于对企业的所有数据进行汇总&#xff0c;为企业各个部门提供一个统一、规范的出口。做数仓就是做方案&#xff0c;是用数据治理企业的方案。 数据仓库的特点 面向主题集成 公司中不同的部门都会去数据仓库中拿数据&#xff0c;把独…

框架概述和MyBatis环境搭建

学习视频&#xff1a;1001 框架概述_哔哩哔哩_bilibili 目录 框架概述 1.1为什么要学 1.2框架的优点 1.3 当前主流框架 Spring框架 Spring MVC框架 MyBatis框架 ​编辑 Spring Boot框架 Spring Cloud框架 1.4 传统JDBC的劣势 MyBatis 2.1 MyBatis概述 ORM框架工作原…

分布式日志

1 日志管理 1.1 日志管理方案 服务器数量较少时 直接登录到目标服务器捞日志查看 → 通过 rsyslog 或shell/python 等脚本实现日志搜集并集中保存到统一的日志服务器 服务器数量较多时 ELK 大型的日志系统&#xff0c;实现日志收集、日志存储、日志检索和分析 容器环境 …

编程语言MoonBit新增矩阵函数的语法糖

MoonBit更新 1. 新增矩阵函数的语法糖 新增矩阵函数的语法糖&#xff0c;用于方便地定义局部函数和具有模式匹配的匿名函数&#xff1a; fn init {fn boolean_or { // 带有模式匹配的局部函数true, _ > true_, true > true_, _ > false}fn apply(f, x) {f(x)}le…

【分布式技术】注册中心zookeeper

目录 一、ZooKeeper是什么 二、ZooKeeper的工作机制 三、ZooKeeper特点 四、ZooKeeper数据结构 五、ZooKeeper应用场景 ●统一命名服务 ●统一配置管理 ●统一集群管理 ●服务器动态上下线 ●软负载均衡 六、ZooKeeper的选举机制 七、实操部署ZooKeeper集群 步骤一…

compose部署tomcat

1.部署tomcat 1.1.下载相关镜像tomcat8.5.20 $ docker pull tomcat:8.5.20 1.2 在/data目录下创建tomcat/webapps目录 mkdir -p /data/tomcat/webapps 注意&#xff1a;这里是准备将宿主机的/data/tomcat/webapps映射到容器的 /usr/…

LC 2765. 最长交替子数组

2765. 最长交替子数组 难度: 简单 题目大意: 给你一个下标从 0 开始的整数数组 nums 。如果 nums 中长度为 m 的子数组 s 满足以下条件&#xff0c;我们称它是一个 交替子数组 &#xff1a; m 大于 1 。s1 s0 1 。下标从 0 开始的子数组 s 与数组 [s0, s1, s0, s1,...,s(…

HarmonyOS鸿蒙学习笔记(22)@Builder实战

Builder标签是一种更轻量的UI元素复用机制&#xff0c;下面通过一个简单的例子来具体说明&#xff1a; 比如如下布局效果&#xff1a;上面是一个轮播的Swiper,下面是一个Grid 布局代码如下&#xff1a; build() {Navigation() {Scroll() {Column({ space: CommonConstants.CO…

测试老司机聊聊测试设计都包含什么?

一、数据组合测试设计 数据组合测试设计&#xff08;Combinatorial Test Design&#xff0c;CTD&#xff09;是一种优化测试用例的方法&#xff0c;它通过系统地组合不同的测试数据输入&#xff0c;以确保全面覆盖各种可能的测试情况。这种方法主要应用于软件测试领域&#xff…

性能优化-HVX 开发环境介绍

「发表于知乎专栏《移动端算法优化》」 本篇以 HVX 的开发环境配置以及应用实例编译测试为主进行讲述。 &#x1f3ac;个人简介&#xff1a;一个全栈工程师的升级之路&#xff01; &#x1f4cb;个人专栏&#xff1a;高性能&#xff08;HPC&#xff09;开发基础教程 &#x1f3…

scanf解决遇到空格停止问题

scanf解决遇到空格停止问题 gets修改scanf的停止符 我们经常输入字符串的时候&#xff0c;遇到空格&#xff0c;scanf就会停止&#xff1a; 比如这时候我想输入一个句子&#xff1a;“My Love”&#xff1a; char* s (char*)malloc(sizeof(char)*100);scanf("%s", s…

Python数据分析:基础

数据分析是使用统计和计算机科学技术来解决现实世界问题的过程。Python是一种流行的编程语言&#xff0c;也是数据分析的常用工具之一。在Python中&#xff0c;有许多库和工具可用于数据分析&#xff0c;比如NumPy、Pandas、Matplotlib等。 NumPy是Python中用于科学计算的一个…

python 32位图转换8位图并压缩

安装 python -m pip install opencv-python 脚本 from PIL import Image import numpy as np import cv2 oldName "Down32.png" tempName "Down8.png" newName "Down.png" img Image.open(oldName).convert(P) print(img.getbands()) #…

摄像头电机马达驱动芯片LV8548/LV8549/ONSEMI替代料GC8548

摄像头电机马达驱动芯片GC8548&#xff0c;兼容替代 ON的LV8548 无需更改外围 . 下图为其参数分析&#xff1a; GC8548 是一款双通道 12V 直流电机驱动芯片&#xff0c;为摄像机、消费类产品、玩具和其他低压或者电池供电的运动控制类应用提供了集成的电机驱动解决方案。芯片…

hpa自动伸缩

1、定义&#xff1a;hpa全称horizontal pod autoscaling&#xff08;pod的水平自动伸缩&#xff09;&#xff0c;这是k8s自带的模块。pod占用CPU的比率到达一定阀值会触发伸缩机制&#xff08;根据CPU使用率自动伸缩&#xff09; replication controller副本控制器&#xff0c…

带头 + 双向 + 循环链表增删查改实现

目录 源码&#xff1a; List.c文件&#xff1a; List.h文件&#xff1a; 简单的测试&#xff1a; 很简单&#xff0c;没什么好说的&#xff0c;直接上源码。 源码&#xff1a; List.c文件&#xff1a; #include"DLList.h"ListNode* creadNode(LTDataType x) {L…

力扣!30天60道(第2天)

第1题(1.22) &#xff1a;两数之和 解法一&#xff1a;暴力破解 #include <iostream> #include <vector> #include <map> using namespace std;class Solution { public:vector<int> twoSum1(vector<int>& nums, int target) {for (int i …

Java8 安装

> 新版任你发&#xff0c;我用JAVA8 1. 下载jdk包 2. 解压到自己要的位置 3. 配置环境变量 export JAVA_HOME/usr/lib/jvm/jdk1.8.0_202/ export JRE_HOME${JAVA_HOME}/jre export CLASSPATH.:${JAVA_HOME}/lib:${JRE_HOME}/lib export PATH${JAVA_HOME}/bin:$PATH