后缀数组学习小结

学习主要是参考 罗穗骞 的《后缀数组——处理字符串的有力工具》论文
然后自己结合其他书上的后缀数组的实现的方法，选择了一个比较简单的实现版本。
我的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次的最长重复子串

多个字符串一般是用不同的分隔符，把他们连接起来
然后求后缀数组，与高度数组，在按照高度数组分组等
详细解法，看论文。