[后缀数组]【学习笔记】

飒飒飒飒飒飒飒飒飒飒

---

研究了好长时间....(诶好像莫比乌斯反演时也说过这句话)

参考资料：

1.http://wenku.baidu.com/link?url=Beh6Asxvtm7M2QY5kiPyKKaP87xvBrNBKW9LXOeGKm-WM4GoUM3opnHZ8z-DahF7TRaLZZ4cpUe6jfFF064XUEmAiIDF7t90CpgNfSC3_Pq

2.http://www.cnblogs.com/staginner/archive/2012/02/02/2335600.html

3.http://www.xymlife.com/2016/01/16/后缀数组小结/#more

【概念】：（论文里很清楚了，这里说自己的理解）

1. 后缀：后缀是指从某个位置i开始到结束的子串。Suffix(i)=r[i..len(r)]

后缀i是从i开始的子串

2. 后缀数组： 后缀数组 SA 是一个一维数组，它保存1..n 的某个排列 SA[1]，SA[2]， ……，SA[n]，并且保证Suffix(SA[i]) < Suffix(SA[i+1])，1 ≤ i < n。也就是将 S 的 n 个后缀从小到大进行排序之后把排好序的后缀的开头位置顺次放入SA 中。

后缀排序之后的结果

sa[i]是排名为i的后缀的起点下标(开始位置)

3. 名次数组：名次数组 Rank[i] 保存的是 Suffix(i) 在所有后缀中从小到大排列的 “ 名次 ” 。
rank[i]是后缀i的名次(从i开始的子串 位置i的子串)

---

【倍增算法】：

见论文描述

一张图，每个位置开始长度为j的排序，然后用两个j组成2*j排序.........

---

【实现】：

排序使用计数排序，

n是字符串长度，m是字符大小上限

a是要排序数组 sa后缀数组

c是计数排序用的 

r为“第一关键字key1”的名次数组(类似于给每个位置开始的子串分配一个名次/权值 一开始一个字符r[i]=a[i])

k为“第二关键字key2”的排序结果 

k[i]表示的是按key2排序后第i名整个子串(key1+key2)开始的位置，不是key2开始的位置

可以发现 key1开始位置+j=key2开始位置

 

先给单个字符计数排序

然后从j=1开始(每次循环是把两个j合成j<<1)

p最后是当前不同的子串(不一定是后缀)的个数，p==n就没必要继续排序了(已经可以区分每个位置开始的字符串的大小关系了)

1.给key2排序：先处理没有key2的(n-j+1...n)，再通过上一次排序后的sa得到k(上次sa就是长度j的结果，这次的一个关键字长度也是j，只要sa[i]-j变成key1开始的位置就行了)

2.这时r[k[i]]就是key2排第i名的key1的权值，k[i]就是i的开始位置，计数排序，结果用sa保存

3.更新r为长度j<<1时的名次，把r复制给k，这个名次会有相同的，应该算一个，用离散化去重的写法....如何判断相同:先判断是不是都有key2，有一个没有key2的话两个一定不相同(长度就不可能相同了啊)，再分别判断两个关键字是否相同(这时的k是长度j时的名次，判断就行了)

「这里和论文里不太一样，他的做法是加了一个0....」

---

【求LCP】

height数组:定义 height[i]=suffix(sa[i-1])和 suffix(sa[i])的最长公共前缀

排名i和排名i-1后缀的LCP

对于第i名和第j名，他们的LCP就是RMQ{height[i+1],height[i+2],...,height[j]}

定义h[i]=height[rank[i]]

第i个后缀和在它的前一名的后缀的LCP。 

性质：h[i]>=h[i-1]-1

设 suffix(k)是排在 suffix(i-1)前一名的后缀，则它们的最长公共前缀是 h[i-1]。那么 suffix(k+1)将排在 suffix(i)的前面(这里要求 h[i-1]>1，如果 h[i-1]≤1，原式显然成立 )并且 suffix(k+1)和 suffix(i)的最长公共前缀是 h[i-1]-1，所以 suffix(i)和在它前一名的后缀的最长公共前缀至少是 h[i-1]- 1。按照h[1],h[2],......,h[n]的顺序计算,并利用 h 数组的性质，时间复杂度可 以降为 O(n)。 

也就是说第i个和他的前一名的LCP至少是i-1和它的前一名的LCP-1，直观上理解，i-1往后一位(长度-1)就是i啊

实现上，并没有必要弄一个h[]，只需要按照h的顺序(也就是原来1..n的顺序)算就行了

---

注意：m的大小最后会为n，所以c数组至少为n 

    void getHeight(){
        int k=0;
        for(int i=1;i<=n;i++) rnk[sa[i]]=i;
        for(int i=1;i<=n;i++){
            if(k) k--;
            if(rnk[i]==1) continue;
            int j=sa[rnk[i]-1];
            while(i+k<=n&&j+k<=n&&s[i+k]==s[j+k]) k++;
            height[rnk[i]]=k;
        }
    }
    inline bool cmp(int *r,int a,int b,int j){
        return a+j<=n&&b+j<=n&&r[a]==r[b]&&r[a+j]==r[b+j];
    }
    void getSA(){ m=300
        int *r=t1,*k=t2;
        for(int i=0;i<=m;i++) c[i]=0;
        for(int i=1;i<=n;i++) c[r[i]=s[i]]++;
        for(int i=1;i<=m;i++) c[i]+=c[i-1];
        for(int i=n;i>=1;i--) sa[c[r[i]]--]=i;
        
        for(int j=1;j<=n;j<<=1){
            int p=0;
            for(int i=n-j+1;i<=n;i++) k[++p]=i;
            for(int i=1;i<=n;i++) if(sa[i]>j) k[++p]=sa[i]-j;
            
            for(int i=0;i<=m;i++) c[i]=0;
            for(int i=1;i<=n;i++) c[r[k[i]]]++;
            for(int i=1;i<=m;i++) c[i]+=c[i-1];
            for(int i=n;i>=1;i--) sa[c[r[k[i]]]--]=k[i];
            
            swap(r,k);p=0;r[sa[1]]=++p;
            for(int i=2;i<=n;i++) r[sa[i]]=cmp(k,sa[i],sa[i-1],j)?p:++p;
            if(p>=n) break;m=p;
        }
    }
    
    int lcp(int x,int y){
        x=rnk[x];y=rnk[y];
        if(x>y) swap(x,y);x++;
        int t=Log[y-x+1];
        return min(mn[x][t],mn[y-Pow[t]+1][t]);
    }