回文自动机:从入门到只会打板

写在前面

如果你会SAMSAMSAM,相信回文自动机不会难懂。

如果你不会,你可以参考我的上一篇文章。

至少回文自动机是治愈系的吧。

作用

回文自动机,也叫回文树,简称PAMPAMPAM实际上它既不是自动机也不是树

处理回文串的有力工具。可完全代替ManacherManacherManacher,就是会多个字符集。

算法流程

首先一个结论:往一个串末尾加入一个字符,最多会增加一个没出现过的回文串。

证明:

假设出现了两个之前没有出现的,它们都是当前串的后缀。

根据对称性,把短的关于长的的中心作对称,得到一个相同的回文串,说明已经出现过了。

矛盾,所以最多出现一个。

这么说,对于字符串SSS本质不同的回文串最多有∣S∣|S|S个。

于是可以每个状态表示一个本质相同的回文串。用lenilen_ileni表示状态iii表示的回文串的长度。

考虑到回文串两边一样的,定义转移ccc表示当前串左右各加上一个ccc

因为分奇回文和偶回文,所以定义两个初始节点0,10,10,1,其中len[0]=0,len[1]=−1len[0]=0,len[1]=-1len[0]=0,len[1]=1(让len[1]=−1len[1]=-1len[1]=1可以说是这个算法最妙的地方)

000表示偶回文,111表示奇回文。−1-11可以看做添加时吞掉一个字符。

在这里插入图片描述
类比ACACAC自动机,定义fail[i]fail[i]fail[i]这个回文串的最长非自己的回文后缀(以下简称回文真后缀)。

强行让fail[0]=fail[1]=1fail[0]=fail[1]=1fail[0]=fail[1]=1

实际上有个隐含条件:1的所有转移都到0
在这里插入图片描述
接下来同样考虑增量。

当我们插入S[i]S[i]S[i]时,现在我们知道S[i−1]S[i-1]S[i1]的所有回文后缀,要求S[i]S[i]S[i]的回文后缀。

我们可以从前一个点跳failfailfail,如果到某个点ppp刚好可以前后接上S[i]S[i]S[i],即S[i−len[p]−1]=S[i]S[i-len[p]-1]=S[i]S[ilen[p]1]=S[i],说明有一个ppp到当前点的转移。

如果已经有转移了,说明这个回文串出现过,直接退出。

接下来维护failfailfail。我们发现最长回文真后缀和它本身具有相同的性质。

在之前的基础上继续跳就可以了。如果跳到某个pppS[i]S[i]S[i]的转移,说明ch[p][S[i]]ch[p][S[i]]ch[p][S[i]]是个回文后缀,连过去就可以了。

栗子:AABAAABAAABA

插入AAA

此时len[1]=−1len[1]=-1len[1]=1的优势就体现出来了,因为刚好是S[i]=S[i]S[i]=S[i]S[i]=S[i]

在这里插入图片描述
插入AAA,依次跳到0,10,10,1
在这里插入图片描述
插入BBB

在这里插入图片描述插入AAA

在这里插入图片描述

实现

第一步插入的时候由于一些玄学问题,pppfailfailfail可能接到自己身上

解决策略是先把fail算出来,再接到之前的节点后面

剩下的就很容易了

代码有点古怪,仅供参考

char s[MAXN];
int n;
int ch[MAXN][26],fail[MAXN];
int len[MAXN];
int las=1,tot=1;
void init()
{len[1]=-1;fail[0]=fail[1]=1;
}
void insert(int i)
{int p=las;while (s[i-len[p]-1]!=s[i]) p=fail[p];if (ch[p][s[i]-'a']) return (void)(las=ch[p][s[i]-'a']);int q=fail[p];while (s[i-len[q]-1]!=s[i]) q=fail[q];las=++tot;fail[las]=ch[q][s[i]-'a'];	len[ch[p][s[i]-'a']=las]=len[p]+2;
}

运用

①每个点结尾的回文个数

failfailfail树的深度

②本质不同的回文个数

就是状态数

好像只有这些……

想到再补吧

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

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

相关文章

2021牛客第一场 I. Increasing Subsequence-前缀和优化dp

https://ac.nowcoder.com/acm/contest/11166/I 思路:dp[i][j] 是表示上上步走在i点,上一步走在j点的期望。首先我们很容易想到n^3的做法,那我们必须考虑去优化一维的时间复杂度。我们可以考虑使用前缀和优化dp转移。 我们枚举i点&#xff0c…

Codeforces Round #620 (Div. 2) E. 1-Trees and Queries 思维 + LCA

传送门 文章目录题意思路&#xff1a;题意 思路&#xff1a; 照例&#xff0c;先考虑不加边怎么做。由于可以经过重复的边或点&#xff0c;设aaa与bbb之间长度为lenlenlen&#xff0c;那么需要len<klen<klen<k并且还需要(k−len)mod20(k-len) \bmod 20(k−len)mod20&…

.NET Core 微服务之Polly熔断策略

紧接着上一篇说&#xff0c;咱们继续介绍Polly这个类库熔断策略&#xff08;Circuit-breaker&#xff09;如果调用某个目标服务出现过多超时、异常等情况&#xff0c;可以采取一定时间内熔断该服务的调用&#xff0c;熔断期间的请求将不再继续调用目标服务&#xff0c;而是直接…

【洛谷P4169】天使玩偶/SJY摆棋子【CDQ分治】

传送门 题意&#xff1a;动态加点&#xff0c;给定点询问曼哈顿距离最近的点 N,M≤3e5,x,y≤1e6N,M \leq 3e5,x,y \leq 1e6N,M≤3e5,x,y≤1e6 经(kan)过(le)分(ti)析(jie),这是一道cdqcdqcdq分治 考虑当前区间左半边修改对右半边的询问的影响 设左边某个修改为(x1,y1)(x_1,…

牛客第二场 G.League of Legends-单调队列优化dp

https://ac.nowcoder.com/acm/contest/11253/G 上面出题人给的题解&#xff1a; 思路基本差不多&#xff0c;这里主要说一下合并小区间的dp&#xff0c; dp[i][j]代表前i个分成j组最大的时间max 我们首先将区间排好序&#xff0c;如果满足a[k]>b[i] ,则有 j都是由j-1转…

.NET中扩展方法和Enumerable(System.Linq)

LINQ是我最喜欢的功能之一&#xff0c;程序中到处是data.Where(xx>5).Select(x)等等的代码&#xff0c;她使代码看起来更好&#xff0c;更容易编写&#xff0c;使用起来也超级方便&#xff0c;foreach使循环更加容易&#xff0c;而不用for int..&#xff0c;linq用起来那么爽…

Peaks加强版 黑暗爆炸 - 3551 Kruskal重构树 + 主席树

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一张图&#xff0c;有nnn个山峰&#xff0c;每个山峰高度为hih_ihi​&#xff0c;有mmm条边&#xff0c;每条边有个难度值wiw_iwi​&#xff0c;现在有qqq个询问&#xff0c;每次询问给定一个山峰vvv&am…

【BJOI2017】树的难题【点分治】【线段树】

传送门 传送门 题意&#xff1a;给一棵树&#xff0c;树上有颜色&#xff0c;每种颜色有权值&#xff0c;定义一条路径的权值为所有颜色相同段的权值之和&#xff0c;求长度在[L,R][L,R][L,R]中的路径的最大权值。 数据范围&#xff1a;暴力过不了 显然是个点分治 对于分治…

Planes, Trains, but not Automobiles-求最小路径覆盖的起点终点

https://vjudge.net/problem/Kattis-planestrainsbutnotautomobiles 题意&#xff1a;给一个有向图&#xff0c;火车可以由任意一个起点开始&#xff0c;每一个点只能经过一次&#xff0c;在坐火车的时候你可以选择坐飞机到另外一个点&#xff0c;求坐飞机的最小次数&#xff…

译 | .NET Core 3.0 对诊断的改进

原文&#xff1a;Sourabh Shirhatti翻译&#xff1a;Edi Wang在 .NET Core 3.0 中&#xff0c;我们将引入一套工具&#xff0c;这些工具利用 .NET 运行时中的新功能&#xff0c;使诊断和解决性能问题变得更加容易。这些运行时功能可帮助您回答一些常见的诊断问题&#xff1a;我…

P4768 [NOI2018] 归程 Kruskal重构树 + 倍增 + 最短路

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个联通无向图&#xff0c;每条边有一个长度lll和海拔aaa&#xff0c;当海拔≤\le≤水位线的时候&#xff0c;说明这个道有积水。在起始点有一辆车&#xff0c;车可以走没有积水的路&#xff0c;下车后…

【AT987】高橋君【组合数】【莫队】

传送门 题意&#xff1a;TTT组询问NNN个相同物品选不超过KKK个的方案数&#xff0c;T,N≤1e5T,N \leq 1e5T,N≤1e5 设f(x,y)∑i0yCxif(x,y)\sum_{i0}^{y}C_x^if(x,y)∑i0y​Cxi​即所求 直接求并没有很好的性质 但我们发现:f(x,y)2f(x−1,y)−Cx−1yf(x,y)2f(x-1,y)-C_{x-1…

2021牛客第一场H.Hash Function—FFT求差值的卷

https://ac.nowcoder.com/acm/contest/11166/H 官方题解。 比赛时&#xff0c;我们都是用暴力写的&#xff0c;数据太弱了&#xff0c;今天突然想起来&#xff0c;用fft写了一下。 主要使用fft求差值的卷&#xff0c;乘法我们在多项式的系数上做&#xff0c;加减法 我们在多项…

设计模式之策略者模式

策略者模式简介策略者模式定义一个算法接口&#xff0c;并由其实现类去实现&#xff0c;使得每一个算法都得到封装&#xff0c;并让他们可以相互替换。这是一种行为型模式。策略者模式降低了算法行为和环境角色的耦合度&#xff0c;使得算法可以独立发生变化。策略者模式在现实…

Codeforces Round #622 (Div. 2) D. Happy New Year 状压dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; n≤1e5,m≤1e9,k≤8.n\le 1e5,m\le 1e9,k\le 8.n≤1e5,m≤1e9,k≤8. 思路&#xff1a; 注意到题目中保证了每个孩子至多收到kkk个&#xff0c;且k≤8k\le 8k≤8&#xff0c;注意到这是题目保证的&#xff0…

【NOI2015】品酒大会【后缀数组】【并查集】

传送门 传送门 题意&#xff1a;给一个长度为NNN的字符串和一个长度为NNN的序列AAA&#xff0c;对于所有的k∈[0,N−1]k \in [0,N-1]k∈[0,N−1]&#xff0c;求选出两个数i,ji,ji,j满足lcp(suffix(i),suffix(j))≥klcp(suffix(i),suffix(j))\geq klcp(suffix(i),suffix(j))≥…

2021牛客第四场-E Tree Xor-线段树区间异或

https://ac.nowcoder.com/acm/contest/11255/E 题目大意&#xff1a;有n个点&#xff0c;每一点有一个区间[l,r]范围的权值可能&#xff0c;有n-1条边&#xff0c;边上的权值为w,ww(u)^w(v),问你满足权值范围和边权值的方案数。 思路&#xff1a;我们把树根的值设为0&#xff…

Codeforces Round #624 (Div. 3) D. Three Integers 数论

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给定A,B,CA,B,CA,B,C&#xff0c;找到一个三元组(a,b,c)(a,b,c)(a,b,c)&#xff0c;使得aaa是bbb的因子&#xff0c;bbb是ccc的因子&#xff0c;且abs(A−a)abs(B−b)abs(C−c)abs(A-a)abs(B-b)abs(C-c)abs(…

Gym - 100917F Find the Length-用最小路径树求最小环

https://codeforces.com/gym/100917/problem/F 题目大意&#xff1a;就给你一个无向图&#xff0c;求出所有点属于的最小环大小。 思路&#xff1a;我们给每一个点用dji求出最小路径树&#xff0c;然后对于一颗最小路径树&#xff0c;求环。 和树根构成环只有两种情况&#x…

【SPOJ2666】QTree4【链分治】

传送门 题意&#xff1a;给一棵带边权的树&#xff0c;每个点开始时为白色&#xff0c;维护两种操作&#xff1a; 1.改变一个点的颜色&#xff08;白变黑&#xff0c;黑变白&#xff09; 2.询问最远的两个白点之间的距离 树分治国集论文 链分治的本质其实就是树链剖分。它们…