学习主要是参考 罗穗骞 的《后缀数组——处理字符串的有力工具》论文
然后自己结合其他书上的后缀数组的实现的方法,选择了一个比较简单的实现版本。
我的SA模板
*
我的后缀数组和论文的后缀数组
计算出来的lcp高度数组有一个下标的偏移
我的高度数组是lcp数组,也就是[0,n)是有效值
论文的后缀数组在计算sa数组的时候是在字符串的末尾加上一个0
我的这个模板不用这样操作(因为在计算的时候特判了)
*/
int n,k;
int rank[maxn],tmp[maxn],sa[maxn],lcp[maxn];
char s[maxn];
int cmp(int i,int j){
if(rank[i]!=rank[j])return rank[i]<rank[j];
int ri=i+k<=n?rank[i+k]:-1;
int rj=j+k<=n?rank[j+k]:-1;
return ri<rj;
}
void getSa(){
n=strlen(s);
for(int i=0;i<=n;i++){
sa[i]=i;
rank[i]=i<n?s[i]:-1;
}
for(k=1;k<=n;k<<=1){
sort(sa,sa+n+1,cmp);
tmp[sa[0]]=0;
for(int i=1;i<=n;i++){
tmp[sa[i]]=tmp[sa[i-1]]+cmp(sa[i-1],sa[i]);
}
for(int i=0;i<=n;i++)rank[i]=tmp[i];
}
}
void getLcp(){
// n=strlen(s);
for(int i=0;i<=n;i++)rank[sa[i]]=i;
int h=0;
lcp[0]=0;
for(int i=0;i<n;i++){
int j=sa[rank[i]-1];
if(h>0)h--;
while(i+h<n&&j+h<n&&s[i+h]==s[j+h])h++;
lcp[rank[i]-1]=h;
}
}
后缀数组应用:
1,字符串的匹配,利用二分进行后缀的匹配 ,时间O(TlogS)
2,找回文串(也可用简单的manacher)
3,可以重叠的最长重复子串,按照lcp分组,二分答案
4,无重叠的最长重复子串,做法基本如上
5,不同的子串的个数
6,可重叠的出现k次的最长重复子串
多个字符串一般是用不同的分隔符,把他们连接起来
然后求后缀数组,与高度数组,在按照高度数组分组等
详细解法,看论文。