\(1\) Manacher

挺短,背是挺好背的

Manacher用于求回文串长度。思想大概就是:

1、加入字符集之外的识别字符(比如#)分隔开原来相邻的字母,这样所有的回文串都变成了以某个字符为中心的(否则如果是偶数长度的回文串还要特判)。

2、考虑借由以前的信息求出新的回文串长度。记录到现在为止最靠右的回文串中最右侧的字符下标&其对称轴的下标,不妨记这个最靠右的串为\(\rm S\)。那么考虑以当前位置作为对称轴的答案,一定至少是\(\min\){隔着\(\rm S\)的对称轴与其对称的另一个位置ans,\(|S|-i+1\)} 。然后就不断扩展即可。

3、关于复杂度证明。我们记一次帅气的操作的意义是成功让\(ans_i\)的初始值继承了与之对称的点的答案和边界的取\(\min\),记以当前点为轴的最长回文子串为\(\rm T\),\(T\)的右端点为\(q\)。可以知道

  • (1)\(\rm S\)的右端点是单增的;
  • (2)如果当前旧的\(maxlen<i\),即未成功进行一次帅气的操作,那么显然while1次,\(maxlen\)增大一次;
  • (3)如果当前的串经过了一次帅气的操作,那么当\(q<maxlen\)时,直接跳出while;当\(q\geq maxlen\)时,\(q\)增大\(maxlen\)必增大。所以得出结论,进行一次帅气的操作和\(maxlen\)的增大次数是严格同阶的。

So,最终复杂度就是\(\Theta(n)\)的。

void Manacher(char *s){
int id, fars, i ;
id = 0, fars = 0 ;
//id : 最靠右的回文串的中心位置
//fars : 迄今为止最靠右的回文串的最右侧
for (i = 1 ; i <= N ; ++ i)
ns[++ L] = (int)In[i], ns[++ L] = '#' ;
for (i = 1 ; i <= L ; ++ i){
if (fars <= i) base[i] = 1 ;
else base[i] = min(fars - i + 1, base[id * 2 - i]) ;
while (ns[i + base[i]] == ns[i - base[i]]) base[i] ++ ;
if (i + base[i] > fars) id = i, fars = i + base[i] - 1 ;
}
}
int main(){
scanf("%s", In + 1),
L = -1, N = strlen(In + 1) ;
ns[++ L] = '$', ns[++ L] = '#' ; Manacher(In) ;
for (int i = 1 ; i <= 2 * N + 2 ; ++ i) ans = max(ans, base[i] - 1) ;
cout << ans << endl ; return 0 ;
}

\(2\) PAM

学了PAM,不知道为啥感觉比SAM简单?参考的资料会放在最后。

其实就是一种自动机,以回文串为状态,左右各添加一个字符为转移的自动机。要点如下:

0、一个串的回文子串至多有\(O(n)\)个。

1、首先每个节点需要保存这个节点中回文串的长度。

2、显然始状态需要有两个,即奇数长度的\(s\)和偶数长度的\(s\),称作“奇根”“偶根”。那么为了方便呢,奇根的长度设置为\(-1\),偶根长度设置为\(0\)。

3、考虑要从\(last\)指针扩展当前状态,假设当前需要insert的字母是\(c\),是这个串里面的第\(p\)个字符,那我们需要找到一个后缀\(s[j...p-1]\quad s.t.\quad s[j...p-1]\)本身回文且\(s[j-1]=c\),那么就可以向下扩展。

4、考虑怎么找这个后缀,显然对于一个串\(S\),他的所有回文后缀都是其最长回文后缀的回文后缀。所以考虑\(fail\)指针,应当从当前状态连向它的最长回文后缀

5、插入新节点时,考虑跳完\(fail\)后如果没有相应的转移边,就要新建一个状态然后连\(fail\).

然后是代码和一点注意:

struct PAM{
int trie[MAXN][Sigma] ;
int rt0, rt1, last, sz ;
int len[MAXN], fail[MAXN] ;
}P ;
void _init(PAM &p){
p.sz = -1,
p.rt0 = ++ p.sz, p.rt1 = ++ p.sz ;
p.fail[p.rt0] = p.fail[p.rt1] = p.rt1 ;
p.last = p.rt0, p.len[p.rt0] = 0, p.len[p.rt1] = -1 ;
}
void _insert(PAM &p, int x, int pos, char *s){
int u = p.last ;
while (s[pos - p.len[u] - 1] != s[pos]) u = p.fail[u] ;
if (!p.trie[u][x]){
int fa = p.fail[u] ;
int newn = ++ p.sz ;
p.len[newn] = p.len[u] + 2 ;
while (s[pos - p.len[fa] - 1] != s[pos]) fa = p.fail[fa] ;
p.fail[newn] = p.trie[fa][x], p.trie[u][x] = newn,
}
p.last = p.trie[u][x] ;
}

6、\(\rm \color{red}{WARNING}\),以下两句顺序不要写反:

p.fail[newn] = p.trie[fa][x], p.trie[u][x] = newn,

原因是当\(fa=u\)时就出现环了。

\(3\) 闲扯

学完才知道,\(\rm PAM\)又简单又好背功能又多……Manacher被打爆了啊喂qwq。

[学习笔记] Manacher与PAM的更多相关文章

  1. 学习笔记 - Manacher算法

    Manacher算法 - 学习笔记 是从最近Codeforces的一场比赛了解到这个算法的~ 非常新奇,毕竟是第一次听说 \(O(n)\) 的回文串算法 我在 vjudge 上开了一个[练习],有兴趣 ...

  2. Manacher算法学习笔记 | LeetCode#5

    Manacher算法学习笔记 DECLARATION 引用来源:https://www.cnblogs.com/grandyang/p/4475985.html CONTENT 用途:寻找一个字符串的 ...

  3. 【学习笔记】字符串—马拉车(Manacher)

    [学习笔记]字符串-马拉车(Manacher) 一:[前言] 马拉车用于求解连续回文子串问题,效率极高. 其核心思想与 \(kmp\) 类似:继承. --引自 \(yyx\) 学姐 二:[算法原理] ...

  4. C++学习笔记(2)

    本学习笔记是C++ primer plus(第六版)学习笔记.是C++学习笔记(1)的后续.复习C++基础知识的可以瞄瞄. 转载请注明出处http://www.cnblogs.com/zrtqsk/p ...

  5. Beego学习笔记——Logs

    日志处理 这是一个用来处理日志的库,它的设计思路来自于database/sql,目前支持的引擎有file.console.net.smtp,可以通过如下方式进行安装: go get github.co ...

  6. OI知识点|NOIP考点|省选考点|教程与学习笔记合集

    点亮技能树行动-- 本篇blog按照分类将网上写的OI知识点归纳了一下,然后会附上蒟蒻我的学习笔记或者是我认为写的不错的专题博客qwqwqwq(好吧,其实已经咕咕咕了...) 基础算法 贪心 枚举 分 ...

  7. 后缀自动机&回文自动机学习笔记

    在学了一天其实是边学边摆之后我终于大概$get$后缀自动机了,,,就很感动,于是时隔多年我终于决定再写篇学习笔记辽$QwQ$ $umm$和$FFT$学习笔记一样,这是一篇单纯的$gql$的知识总结博, ...

  8. js学习笔记:webpack基础入门(一)

    之前听说过webpack,今天想正式的接触一下,先跟着webpack的官方用户指南走: 在这里有: 如何安装webpack 如何使用webpack 如何使用loader 如何使用webpack的开发者 ...

  9. PHP-自定义模板-学习笔记

    1.  开始 这几天,看了李炎恢老师的<PHP第二季度视频>中的“章节7:创建TPL自定义模板”,做一个学习笔记,通过绘制架构图.UML类图和思维导图,来对加深理解. 2.  整体架构图 ...

随机推荐

  1. 修改so库中的依赖名

    修改so库中的依赖名 在ArchLinuxArm上有一些针对aarch64, arm, armeabi-v7a等Android常用架构的so库可以下载,有时候可以省去很多编译时间,且都是编译optim ...

  2. VS2017创建的单元测试不支持顺序测试

      问题:使用IDE创建的单元测试项目,标准引用是,导致不能添加顺序测试,复制其它项目的顺序测试文件进行编辑时,也会提示基于MSTest V2的测试不能用于顺序测试     解决办法: 移除自带的NU ...

  3. leetcode动态规划--基础题

    跳跃游戏 给定一个非负整数数组,你最初位于数组的第一个位置. 数组中的每个元素代表你在该位置可以跳跃的最大长度. 判断你是否能够到达最后一个位置. 思路 根据题目意思,最大跳跃距离,说明可以跳0--n ...

  4. 【cf932E】E. Team Work(第二类斯特林数)

    传送门 题意: 求\(\displaystyle \sum_{i=0}^n{n\choose i}i^k,n\leq 10^9,k\leq 5000\). 思路: 将\(i^k\)用第二类斯特林数展开 ...

  5. BZOJ3894/LG4313 文理分科 新建点最小割

    问题描述 BZOJ3894 LG4313 题解 显然一个人只能选文/理 -> 一个人只能属于文(S).理(T)集合中的一个 可以把选择文得到 \(art\) 的收益看做选择文失去 \(scien ...

  6. Web安全测试学习笔记 - DVWA+PHP环境搭建

    DVWA(Damn Vulnerable Web Application),是一个用PHP编写的,作为Web安全测试练习平台的合法环境(毕竟咱不能为了练习就随便找个网站去攻击...),也就是俗称的靶场 ...

  7. Map 集合 和 String 字符串相互转换工具类

    package com.skynet.rimp.common.utils.util; import java.util.Arrays; import java.util.HashMap; import ...

  8. Python中7个不一样的代码写法

    打印index 对于一个列表,或者说一个序列我们经常需要打印它的index,一般传统的做法或者说比较low的写法: 更优雅的写法是多用enumerate 两个序列的循环 我们会经常对两个序列进行计算或 ...

  9. Selenium模块的安装

    Selenium模块 1.安装selenium python2:pip install selenium python3:pip install selenium 2.设置浏览器驱动 解压后必须与浏览 ...

  10. 【SDUT】【链表】2120 - 数据结构实验之链表五:单链表的拆分

    Problem Description 输入N个整数顺序建立一个单链表,将该单链表拆分成两个子链表,第一个子链表存放了所有的偶数,第二个子链表存放了所有的奇数.两个子链表中数据的相对次序与原链表一致. ...