后缀数组--许智磊

后缀数组芜湖一中许智磊后缀数组——字符串处理中的有力武器后缀树的一个简单而高效的替代品当今字符串处理研究中的热门让我们一同揭开她神秘的面纱后缀数组——定义和符号字符集、字符、字符串都按照惯常的定义字符串S的长度表示为len(S)字符串的下标从1开始到len(S)结束字符串S的第i个字符表示为S[i]从i到j这一段的子串表示为S[i..j]后缀是一种特殊的子串从某个位置i开始到整个串的末尾结束S的从i开头的后缀等价于S[i..len(S)]后缀数组——定义和符号约定一个字符集Σ待处理的字符串约定为S，约定len(S)=n规定S以字符“$”结尾，即S[n]=“$”“$”小于Σ中所有的字符除了S[n]=“$”之外，S的其他字符都属于Σ对于约定的字符串S，其i开头的后缀表示为Suffix(i)后缀数组——定义和符号字符串的大小关系按照通常所说的“字典顺序”进行比较我们对S的n个后缀按照字典顺序从小到大排序将排序后的后缀的开头位置顺次放入数组SA中，称为后缀数组令Rank[i]保存Suffix(i)在排序中的名次，称数组Rank为名次数组后缀数组——构造方法把n个后缀当作n个字符串，按照普通的方法进行排序——O(n2)低效的原因——把后缀仅仅当作普通的、独立的字符串，忽略了后缀之间存在的有机联系。如何构造后缀数组?后缀数组——构造方法倍增算法(DoublingAlgorithm)对字符串u，定义uk=u[1..k],len(u)≥ku,len(u)k定义k-前缀比较关系k，=k和≤k对两个字符串u,v，ukv当且仅当ukvku=kv当且仅当uk=vku≤kv当且仅当uk≤vk后缀数组——构造方法设u=Suffix(i)，v=Suffix(j)后缀u，以i开头后缀v，以i开头对u、v在2k-前缀意义下进行比较kkkk比较红色字符相当于在k-前缀意义下比较Suffix(i)和Suffix(j)比较绿色字符相当于在k-前缀意义下比较Suffix(i+k)和Suffix(j+k)在2k-前缀意义下比较两个后缀可以转化成在k-前缀意义下比较两个后缀i+kj+k后缀数组——构造方法把n个后缀按照k-前缀意义下的大小关系从小到大排序将排序后的后缀的开头位置顺次放入数组SAk中，称为k-后缀数组用Rankk[i]保存Suffix(i)在排序中的名次，称数组Rankk为k-名次数组后缀数组——构造方法利用SAk可以在O(n)时间内求出Rankk利用Rankk可以在常数时间内对两个后缀进行k-前缀意义下的大小比较后缀数组——构造方法如果已经求出Rankk可以在常数时间内对两个后缀进行k-前缀意义下的比较可以在常数时间内对两个后缀进行2k-前缀意义下的比较可以很方便地对所有的后缀在2k-前缀意义下排序采用快速排序O(nlogn)采用基数排序O(n)可以在O(n)时间内由Rankk求出SA2k也就可以在O(n)时间内求出Rank2k后缀数组——构造方法1-前缀比较关系实际上是对字符串的第一个字符进行比较可以直接根据开头字符对所有后缀进行排序求出SA1采用快速排序，复杂度为O(nlogn)然后根据SA1在O(n)时间内求出Rank1可以在O(nlogn)时间内求出SA1和Rank1后缀数组——构造方法直接根据首字符排序SA1Rank1O(nlogn)SA2Rank2O(n)SA4Rank4O(n)SA8Rank8O(n)……O(n)SAmRankmO(n)m=2t且m≥nO(nlogn)的操作一次O(n)的操作[logn]次总复杂度为O(nlogn)后缀数组——构造方法m≥n，SAm=SARankm=Rank我们已经在O(nlogn)的时间内构造出了后缀数组SA和名次数组Rank后缀数组——方法总结利用到后缀之间的联系用k-前缀比较关系来表达2k-前缀比较关系每次可以将参与比较的前缀长度加倍根据SAk、Rankk求出SA2k、Rank2k参与比较的前缀长度达到n以上时结束倍增思想后缀数组——辅助工具仅仅靠后缀数组和名次数组有时候还不能很好地处理问题后缀数组的最佳搭档——LCP定义两个字符串的最长公共前缀LongestCommonPrefixlcp(u,v)=max{i|u=iv}也就是从头开始比较u和v的对应字符持续相等的最远值后缀数组——辅助工具定义LCP(i,j)=lcp(Suffix(SA[i]),Suffix(SA[j]))也就是SA数组中第i个和第j个后缀的最长公共前缀LCPTheorem对任何1≤ij≤nLCP(i,j)=min{LCP(k-1,k)|i+1≤k≤j}称j-i为LCP(i,j)的“跨度”，LCPTheorem意义为：跨度大于1的LCP值可以表示成一段跨度等于1的LCP值的最小值后缀数组——辅助工具4=LCP(i,i+1)3=LCP(i+1,i+2)5=LCP(i+2,j)LCP(i,j)=3Suffix(SA[i])Suffix(SA[j])Suffix(SA[i+1])Suffix(SA[i+2])后缀数组——辅助工具设height[i]=LCP(i-1,i)根据LCPTheoremLCP(i,j)=min{height[k]|i+1≤k≤j}计算LCP(i,j)等价于询问数组height中下标从i+1到j范围内所有元素的最小值经典的RMQ(RangeMinimumQuery)问题!!!线段树、排序树——O(nlogn)预处理,O(logn)每次询问标准RMQ方法——O(n)预处理,O(1)每次询问后缀数组——辅助工具采用一种“神奇的”方法，可以在O(n)时间内计算出height数组采用标准RMQ方法在O(n)时间内进行预处理之后就可以在常数时间内算出任何的LCP(i,j)可以在常数时间内计算出任何两个后缀的最长公共前缀这是后缀数组最常用以及最强大的功能之一后缀数组——应用举例怎样使用后缀数组？后缀数组——应用举例回文串——顺读和倒读完全一样的字符串奇回文串字符串u满足：1.len(u)=p为奇数2.对任何1≤i≤(p-1)/2，u[i]=u[p-i+1]偶回文串字符串v满足：1.len(v)=q为奇数2.对任何1≤i≤q/2，v[i]=v[q-i+1]回文串后缀数组——应用举例字符串T的回文子串——T的子串，并且是回文串字符串T的最长回文子串——T的回文子串中长度最大的给出一个字符串T，求它的最长回文子串给出最大长度即可后缀数组——应用举例分析求最长奇回文子串的算法最长偶回文子串可以类似求出后缀数组——应用举例枚举奇回文串中间一个字符的位置尽量向两边扩展后缀数组——应用举例rrii-r-1i+r+1反射相等以某个位置为中心向两边扩展的复杂度为O(len(T))整个算法的复杂度为O(len(T)2)后缀数组——应用举例求以一个位置i为中心向两边扩展的最远值是算法的核心部分需要降低这一步的复杂度后缀数组——应用举例$中心ii'=2n-i+2i-ri+ri'+r#求以i为中心向两边扩展的最远值，等价于求Suffix(i)和Suffix(i')的最长公共前缀后缀数组！！！同时和粉红串反射相等T串T'串Suffix(i)和Suffix(i')的公共前缀后缀数组——应用举例解法：1.初始化答案为0。按照前述方法修改串T，得到串S2.求出后缀数组SA、名次数组Rank3.计算height数组并进行标准RMQ方法预处理复杂度：O(n)+O(nlogn)+n*O(1)=O(nlogn)4.枚举i，计算以i为对称中心的极长回文串并更新答案+2*O(n)后缀数组VS后缀树后缀树也可以做到类似的事情后缀数组有什么优势呢？后缀数组VS后缀树后缀数组在信息学竞赛中最大的优势：易于理解，易于编程，易于调试后缀数组比后缀树占用的空间少——处理长字符串，如DNA分析后缀数组VS后缀树时间复杂度的比较按照字符总数|Σ|把字符集Σ分为三种类型：ConstantAlphabet——|Σ|是一个常数IntegerAlphabet——|Σ|和字符串长度n规模相当GeneralAlphabet——对|Σ|没有任何限制ConstantAlphabetIntegerAlphabetGeneralAlphabet后缀数组VS后缀树后缀数组是直接针对GeneralAlphabet设计的算法复杂度跟字符集的类型没有关系后缀树则对不同字符集有不同的表现如果采用儿子-兄弟方式来表达后缀树：构造的复杂度为O(n*|Σ|)——显然不适合Integer和GeneralAlphabet，对于|Σ|稍大的ConstantAlphabet也无法胜任解决方法：每个节点建立一棵红黑树来保存儿子，复杂度为O(n*log|Σ|)。——竞赛的时候有时间编吗？结论对于Integer和General以及|Σ|较大的ConstantAlphabet，后缀树甚至在时间复杂度上都无法胜过后缀数组。但是对于|Σ|较小的ConstantAlphabet，后缀树还是有着速度上的优势的。——我们要根据实际情况，因“题”制宜选择合适的数据结构后缀数组——最后的话研究后缀数组，不是因为害怕后缀树的繁琐也没有贬低后缀树，抬高后缀数组的意思对于功能相似的两个数据结构，我们应该灵活地掌握，有比较有选择地使用构造后缀数组用到的倍增思想对我们的思考也是有帮助的后缀数组