关于SAM和广义SAM

不是教程

某些思考先记下来

SAM

终于学会了这个东西诶......

一部分重要性质

确定一个重要事情,S构造出的SAM的一个重要性质是当且仅当对于S的任意一个后缀,可以从1号节点走到终止状态。专业的名词叫做有限状态自动机。

trans[st][c]表示的是对于状态st,如果将st中任意串s加一个c,那么会到达的新状态new,显然new是唯一的。假如不唯一那么s一定不属于同一个st。

fa[st]表示的是对于状态st,如果慢慢缩小st中后缀长度,会到达的第一个状态。规定一个虚点\(1\)全集,那么fa能够构成一个树形结构(parent tree)。

很显然\(endpos[st]\subset endpos[fa[st]]\)。所以虽然\(\sum|endpos|=O(n^2)\)但是可以通过线段树合并得到endpos集合。你也可以发现Parent Tree任何一个到根的链是一个长度递减且连续的某个后缀的前缀。

trans构成的有向图G一定是一个DAG,因为可以建立这个图从1到任意点的路径构成的串对应的原串,所以有环就出大事了。

G上任意一条从1开始走的路径都代表一个后缀的前缀——也就是一个子段。

而对于G任意一个i节点开始走的路径,表示的是对于后缀i的任意一个子段。

可以发现一个确定的Trans[][]可以唯一确定一个串和SAM的形态。因为DAG的最长路是原串。

构建和复杂度证明

先说构建,构建的目的是要让这个东西能够识别新的后缀并且还要维持好Parent Tree的各种性质。

设原来是S[1...r]现在我们要变成S[1...r+1],S[r+1]=c:

首先可以发现的事情是,新状态一定包括\(S[1....r+1]\),设这个新状态为u,那么首先我们要让所有\(S[i...r+1]\)都被表示出来,我们可以一直跳\(v=r\)的fa,根据定义如果没有trans[v][c]我们要更新trans[v][c]。然后一直跳如果直到1才停下来,说明之前\(S[(i \in[1,r])...r]\)都没出现过,那没事了。

但是如果到中间一个\(v\)而trans[v][c]=x,说明\(v\)代表的\(S[i...r+1]\)之前被表示出来过,但按照parent tree的定义我们要更新endspos。此时有一个问题:

x代表的串最长的那个串可能不是我们当前要使得自动机以后识别的串S[i...r+1]。根据parent tree那个一条链的性质,这个集合一定要假如x的所有父亲的endpos集合中,但是可以发现,比S[i...r+1]小的串构成的endpos集合和大于等于S[i...r+1]的串endpos集合现在多了一个r+1,所以我们把x分裂成上下两半y和x',然后让x'和u共同当y的儿子。此时有一个问题是我们要更新一下trans,也就是把之前DAG图上到x的边全部转到y上,根据trans还有fa的性质这些边一定是v的连续一段fa,我们暴力修改一下(复杂度等下分析)。更新后的DAG只能只能识别原来的串+c这些串了(trans全更新了),并且我们成功更新了endpos。


分析一下复杂度:

首先如果trans[v][c]=NULL我们修改成r+1的这一部分复杂度是没问题的。因为你修改一次意味着一条DAG上的边,而DAG上的边是\(O(n)\)的。

唯一的问题是到了v之后,你还要跳连续的v的trans[v'][c]=x的一段链。这里复杂度是多少?

根据trans和Parent Tree分析一下,可以发现,任何一段暴力跳都可以对应上一个另一条树链上的一段,且另一段大小一定不小于这一段。那么一直这样对应下去一定可以找到一个最长的链。由于parent tree总点数是\(O(n)\)的,所以显然暴力跳的次数远远小于\(O(n)\) 。

广义SAM

普通SAM是针对一个串而言的,而广义SAM是针对一个trie树而言的。此时后缀和endpos的定义一一扩大,后缀变成了包括叶子的一段缀,endpos变成了trie树上的节点编号。

考虑一下对于一个串我们弄SAM出来利用的性质,我们只是利用了后缀的后缀也是后缀以及一个节点加上一个字符c后是不是仍然是后缀是确定的,在trie树上,叶子的缀的缀也是叶子的缀是显然的。现在SAM表达的含义是trie树上的节点,trans[u][c]表示的是在u这个点加上一个c字符后会到trie树上哪个点。

考虑以后我们之前的增量构造的含义,意义其实就是更新了到根的一条链的后缀,在这些后缀的后面加上一个字符。

一个构造方法是按照trie树bfs序构造。

关于SAM和广义SAM的更多相关文章

  1. 后缀自动机(SAM)+广义后缀自动机(GSA)

    经过一顿操作之后竟然疑似没退役0 0 你是XCPC选手吗?我觉得我是! 稍微补一点之前丢给队友的知识吧,除了数论以外都可以看看,为Dhaka和新队伍做点准备... 不错的零基础教程见 IO WIKI ...

  2. 【HDU 4436】 str2int (广义SAM)

    str2int Problem Description In this problem, you are given several strings that contain only digits ...

  3. 【BZOJ 3926】 [Zjoi2015]诸神眷顾的幻想乡 (广义SAM)

    3926: [Zjoi2015]诸神眷顾的幻想乡 Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 974  Solved: 573 Descriptio ...

  4. 【BZOJ 3473】 字符串 (后缀数组+RMQ+二分 | 广义SAM)

    3473: 字符串 Description 给定n个字符串,询问每个字符串有多少子串(不包括空串)是所有n个字符串中至少k个字符串的子串? Input 第一行两个整数n,k. 接下来n行每行一个字符串 ...

  5. luogu3346 诸神眷顾的幻想乡 (广义SAM)

    首先,让每一个叶节点做一次树根的话,每个路径一定至少有一次会变成直上直下的 于是对于每个叶节点作为根产生的20个trie树,把它们建到同一个广义SAM里 建法是对每个trie dfs去建,last就是 ...

  6. loj#6031. 「雅礼集训 2017 Day1」字符串(SAM 广义SAM 数据分治)

    题意 链接 Sol \(10^5\)次询问每次询问\(10^5\)个区间..这种题第一感觉就是根号/数据分治的模型. \(K\)是个定值这个很关键. 考虑\(K\)比较小的情况,可以直接暴力建SAM, ...

  7. Luogu P3181 [HAOI2016]找相同字符 广义$SAM$

    题目链接 \(Click\) \(Here\) 设一个串\(s\)在\(A\)中出现\(cnt[s][1]\)次,在\(B\)中出现\(cnt[s][2]\)次,我们要求的就是: \[\sum cnt ...

  8. CF666E Forensic Examination 广义SAM、线段树合并、倍增、扫描线

    传送门 朴素想法:对\(M\)个匹配串\(T_1,...,T_M\)建立广义SAM,对于每一次询问,找到这个SAM上\(S[pl...pr]\)对应的状态,然后计算出对于每一个\(i \in [l,r ...

  9. Luogu4022 CTSC2012 熟悉的文章 广义SAM、二分答案、单调队列

    传送门 先将所有模板串扔进广义SAM.发现作文的\(L0\)具有单调性,即\(L0\)更小不会影响答案,所以二分答案. 假设当前二分的值为\(mid\),将当前的作文放到广义SAM上匹配. 设对于第\ ...

随机推荐

  1. Android学习笔记之 SimpleAdapter 中添加按钮响应事件,getView的重写

    Andriod 里面的ListView是一个显示列表数据的控件,常用适配器SimpleAdapter进行绑定,绑定代码如下: ListView lstView = (ListView) this.fi ...

  2. Android GDI 图形渲染

      发布于2011-07-26   导读:对于Android开发者来说,成系列的技术文章对他们的技术成长帮助最大.如下是我们向您强烈推荐的主题为Android开发的第一个系列文章. <Andro ...

  3. 15-8 pymysql的使用

    一 安装pymysql模块 1 pycharm安装 file-setting如图:然后点加号,搜索pymsql即可,点击安装 2 pip 安装 pip3 install pymysql 二  连接数据 ...

  4. 解锁当前XXX用户

    pam_tally2 查看当前锁账户 pam_tally2 --user=XXX用户 --reset 解锁当前XXX用户

  5. 第三期 预测——Frenet 坐标

    Frenet坐标 在讨论过程模型之前,我们应该提到“Frenet Coordinates”,它是一种以比传统x,y笛卡尔坐标更直观的方式表示道路位置的方式. 用Frenet坐标,我们使用变量 s和d描 ...

  6. @noi.ac - 491@ explore

    目录 @description@ @solution@ @accepted code@ @details@ @description@ 最近有一个巨大的古代地下遗迹在比特镇被发现.这个地下遗迹的俯视图 ...

  7. Ubuntu 开机自动挂载NTFS分区

    先安装两个软件,ntfs-3g 和 ntfs-config,前面那个最新版本的Ubuntu已经自带了,所以只需要安装ntfs-config. 输入 sudo apt-get install ntfs- ...

  8. H3C DCC概念

  9. 洛谷P2672 推销员 题解 贪心

    题目链接:https://www.luogu.org/problem/P2672 这道题目是贪心,贪心的思想是: 选择 \(m\) 户人家的最大疲劳值应该是以下两种方案中的较大值: 方案一:选择 \( ...

  10. jQuery中动态创建、添加元素的方法总结

    <input type="button" value="创建元素" id="btn"> <div id="box ...