回文树/回文自动机

放链接:

状态数的线性证明

并没有看懂上面的证明,所以自己脑补了一个...

引理: 每一个回文串都是字符串某个前缀的最长回文后缀.

证明. 考虑一个回文串在字符串中第一次出现的位置, 记为 \(S_{p_1 ... p_2}\), 它一定是 \(S_{1 ... p_2}\)的最长回文后缀.

否则, 如果有 \(S_{p_3 ... p_2} (p_3<p_1)\) 也为回文串, 那么由于回文, \(S_{p_3 ... p_3-p_2+p_1} = S_{p_1 ... p_2}\), \(S_{p_1 ... p_2}\)并不是它第一次出现的位置.矛盾.

因而命题得证.

而每个点的最长回文后缀是唯一的, 因此\(S\)最多只有\(|S|\)个不同的回文子串.

引理的推论. 一个回文串 \(\leftrightarrow\) 某个串的最长回文子串 && 某个串的最长回文子串的回文后缀.

关于fail指针

fail指针指向的是一个节点代表的回文串的最长回文后缀.

在build时, 它也可以理解为以某个点为结尾的次长回文后缀.

Code

const int ssz=300050;
ll n;
char s[ssz]; struct te{int l,fail,cnt,ch[27];}tree[ssz]{{0,1},{-1,1}};
int pt=1,rt0=0,rt1=1;
#define ch(p,c) tree[p].ch[c]
#define fail(p) tree[p].fail
int newnd(){return ++pt;} int getfail(int p,int i){
while(s[i-1-tree[p].l]!=s[i])p=fail(p);
return p;
}
void build(){
int p,q,last=0;
rep(i,1,n){
p=getfail(last,i);
if(ch(p,s[i])==0){
q=newnd();
tree[q].l=tree[p].l+2,fail(q)=ch(getfail(fail(p),i),s[i]);
ch(p,s[i])=q;
}
last=ch(p,s[i]);
++tree[last].cnt;
}
}

应用

枚举所有回文子串

dfs即可.

拓扑序

显然拓扑序就是 ${ 1, 2, \cdots, n } $.

求字符串出现次数

加入每个字符后, ++cnt[last];;

然后逆拓扑序dp, cnt[fa(p)] += cnt[p].

cnt[p] 即为回文串 \(p\) 出现次数.

详见下面的题.

每个节点长度 \(\le \frac {len}2\) 的回文后缀

和维护fail指针大体类似, 加上限制条件即可.

详见代码.

其中tree[p].tr表示的是 \(p\) 节点长度 \(\le \frac {len}2\) 的回文后缀

struct tnd{int l,fi,ch[csz],tr;}tree[ssz];
#define ch(p,c) tree[p].ch[c]
#define fail(p) tree[p].fi
#define trl(p) tree[p].l
#define trtr(p) tree[p].tr
int rt0=0,rt1=1,pt=1;
int getfail(int p,int i){
while(s[i-1-trl(p)]!=s[i])p=fail(p);
return p;
}
void build(){
int p,q,last=0;
rep(i,1,n){
p=getfail(last,i);
if(ch(p,s[i])==0){
q=++pt;
trl(q)=trl(p)+2,fail(q)=ch(getfail(fail(p),i),s[i]);
//get tr(p) start
if(trl(q)<=1)trtr(q)=fail(q);
else{
int z=trtr(p);
while(s[i-1-trl(z)]!=s[i]||(trl(z)+2)*2>trl(q))z=fail(z);
trtr(q)=ch(z,s[i]);
}
//end
ch(p,s[i])=q;
}
last=ch(p,s[i]);
}
}

例题

BZOJ3676:[Apio2014]回文串

求回文串长度*出现次数的最大值.

板子题.

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<set>
#include<map>
using namespace std;
#define rep(i,l,r) for(register int i=(l);i<=(r);++i)
#define repdo(i,l,r) for(register int i=(l);i>=(r);--i)
#define il inline
typedef double db;
typedef long long ll; //---------------------------------------
const int ssz=300050;
ll n,ans=0;
char s[ssz]; struct te{int l,fail,cnt,ch[27];}tree[ssz]{{0,1},{-1,1}};
int pt=1,rt0=0,rt1=1;
#define ch(p,c) tree[p].ch[c]
#define fail(p) tree[p].fail
int newnd(){return ++pt;} int getfail(int p,int i){
while(s[i-1-tree[p].l]!=s[i])p=fail(p);
return p;
}
void build(){
int p,q,last=0;
rep(i,1,n){
p=getfail(last,i);
if(ch(p,s[i])==0){
q=newnd();
tree[q].l=tree[p].l+2,fail(q)=ch(getfail(fail(p),i),s[i]);
ch(p,s[i])=q;
}
last=ch(p,s[i]);
++tree[last].cnt;
}
} int main(){
ios::sync_with_stdio(0),cin.tie(0);
cin>>(s+1);
n=strlen(s+1);
rep(i,1,n)s[i]-='a'-1;
build();
repdo(i,pt,2){
tree[fail(i)].cnt+=tree[i].cnt;
ans=max(ans,(ll)tree[i].cnt*tree[i].l);
}
cout<<ans<<'\n';
return 0;
}

[模板] 回文树/回文自动机 && BZOJ3676:[Apio2014]回文串的更多相关文章

  1. bzoj3676 [Apio2014]回文串 卡常+SAM+树上倍增

    bzoj3676 [Apio2014]回文串 SAM+树上倍增 链接 bzoj luogu 思路 根据manacher可以知道,每次暴力扩展才有可能出现新的回文串. 所以推出本质不同的回文串个数是O( ...

  2. [Bzoj3676][Apio2014]回文串(后缀自动机)(parent树)(倍增)

    3676: [Apio2014]回文串 Time Limit: 20 Sec  Memory Limit: 128 MBSubmit: 3396  Solved: 1568[Submit][Statu ...

  3. [BZOJ3676][APIO2014]回文串(Manacher+SAM)

    3676: [Apio2014]回文串 Time Limit: 20 Sec  Memory Limit: 128 MBSubmit: 3097  Solved: 1408[Submit][Statu ...

  4. HDU - 5157 :Harry and magic string (回文树,求多少对不相交的回文串)

    Sample Input aca aaaa Sample Output 3 15 题意: 多组输入,每次给定字符串S(|S|<1e5),求多少对不相交的回文串. 思路:可以用回文树求出以每个位置 ...

  5. HDU5658:CA Loves Palindromic (回文树,求区间本质不同的回文串数)

    CA loves strings, especially loves the palindrome strings. One day he gets a string, he wants to kno ...

  6. bzoj千题计划304:bzoj3676: [Apio2014]回文串(回文自动机)

    https://www.lydsy.com/JudgeOnline/problem.php?id=3676 回文自动机模板题 4年前的APIO如今竟沦为模板,,,╮(╯▽╰)╭,唉 #include& ...

  7. 【回文自动机】bzoj3676 [Apio2014]回文串

    回文自动机讲解!http://blog.csdn.net/u013368721/article/details/42100363 pam上每个点代表本质不同的回文子串.len(i)代表长度,cnt(i ...

  8. BZOJ3676[Apio2014]回文串——回文自动机

    题目描述 考虑一个只包含小写拉丁字母的字符串s.我们定义s的一个子串t的“出 现值”为t在s中的出现次数乘以t的长度.请你求出s的所有回文子串中的最 大出现值. 输入 输入只有一行,为一个只包含小写字 ...

  9. BZOJ3676 APIO2014回文串(manacher+后缀自动机)

    由于本质不同的回文子串数量是O(n)的,考虑在对于每个回文子串在第一次找到它时对其暴力统计.可以发现manacher时若右端点移动则找到了一个新回文串.注意这样会漏掉串长为1的情况,特判一下. 现在问 ...

随机推荐

  1. 关于Keychain

    1.Keychain 浅析 2.iOS的密码管理系统 Keychain的介绍和使用 3.iOS开发中,唯一标识的解决方案之keyChain+UUID

  2. Android下获取FPS的几种方法

    FPS(Frames Per Second)是关乎Android用户体验最为重要的指标之一,而在VR中更是如此.为了评估VR系统.VR SDK及Unity应用的性能,通常会实时获取FPS并将其显示出来 ...

  3. Git 最佳实践:分支管理

    5月份,为统一团队git分支管理规范,刚开始准备自己写,在网上搜了下,发现不少不错的git分支管理实践.最后我为团队选择了这个git分支管理实践 A successful Git branching ...

  4. (二)版本控制管理器之CVS(上)

    在前一篇<(一)版本控制管理器之发展史>的介绍中,有提到古典时期的CVS,那什么是CVS?CVS特点是什么?怎么个用法?等一系列的问题,虽然这个版本控制管理器早已过时,但大家了解下也不妨, ...

  5. GenericServlet 、Servlet和httpServler

    -------[转] 1.GenericServlet类是所有Servlet类的祖先类. 2.HttpServlet类继承了GenericServlet类. 3.Servlet有两个非常重要的的对象, ...

  6. 高端内存映射之kmap_atomic固定映射--Linux内存管理(二十一)

    1 固定映射 1.1 数据结构 linux高端内存中的临时内存区为固定内存区的一部分, 对于固定内存在linux内核中有下面描述 x86 arm arm64 arch/x86/include/asm/ ...

  7. wxPython树控件

    1.树控件 树(tree)是一种通过层次结构展示信息的控件,如下图所示是树控件示例,左窗口中是树控件,在wxPython中树控件类是wx.TreeCtrl. wx.TreeCtrl中一个常用的方法有: ...

  8. 公共的JS组件-告别CURD

    urls.py urlpatterns = [ url('^asset.html$', views.AssetView.as_view()), url('^asset-json.html$', vie ...

  9. Docker,Docker Compose,Docker Swarm,Kubernetes之间的区别

    Dcoker Docker 这个东西所扮演的角色,容易理解,它是一个容器引擎,也就是说实际上我们的容器最终是由Docker创建,运行在Docker中,其他相关的容器技术都是以Docker为基础,它是我 ...

  10. Ambari与Kerberos 集成

    Kerberos 介绍 Kerberos 是一个网络认证的框架协议,其设计的初衷便是通过密钥系统为 Client 和 Server 应用程序之间提供强大的认证服务.在使用 Kerberos 认证的集群 ...