尾碼數組（不完善）

前天在看尾碼數組相關的內容時，看了羅穗騫的論文，看了2天才基本看懂了論文中的2倍增法。但是還是沒有理解完全。現在先做一個筆記，以便以後參考。(說實話，我覺得trie， p - trie什麼的都沒這麼傷腦筋)，以下是我在網上down的一個樣本程式。

#include<stdio.h> #include<string.h> #include<algorithm> using namespace std; int cmp(int *r,int a,int b,int l){     return r[a]==r[b]&&r[a+l]==r[b+l]; } //為什麼設定為100結果就不正確呢#define M 200int wa[M],wb[M],wv[M],ws[M]; int SA[M];  char s[M]="aabaaaab";  void print(char *label,int *arr,int n){     puts(label);     for(int i=0;i<n;i++)         printf("%d ",arr[i]);     puts(""); } void da(char *r,int *sa,int n,int m){      int i, j, p, *x = wa, *y = wb;      for (i = 0; i < m; i++)          ws[i] = 0;      for (i = 0; i < n; i++)         ws[x[i] = r[i]]++;      for (i = 1; i < m; i++)          ws[i] += ws[i - 1];      for (i = n - 1; i >= 0; i--)          sa[--ws[x[i]]] = i;      print("array sa:",sa,n); //**************print**************       for (j = 1, p = 1; p < n; j *= 2,m=p)      {          //對第二關鍵字進行排序:由上一次的結果 直接算出來          for (p = 0, i = n - j; i < n; i++) y[p++] = i;          for (i = 0; i < n; i++)              if (sa[i] >= j) y[p++] = sa[i] - j;          print("array y:",y,n); print("array x:",x,n);          //對第一關鍵字進行排序          for (i = 0; i < n; i++) wv[i] = x[y[i]]; //根據第二關鍵字的順序將尾碼數組進行排序以便用於第一關鍵字的排序，注意論文中的圖xy， //先進行y排序，然後進行x排序。注意看論文中的圖2，每次都是先根據上一次的結果算出第二個關鍵字的順序放在y中，然後y的排序結果

//對第一個關鍵字即圖中的x進行排序。         print("wv:",wv,n);//**************************************          for (i = 0; i < m; i++) ws[i] = 0;          for (i = 0; i < n; i++) ws[wv[i]]++;          for (i = 1; i < m; i++) ws[i] += ws[i - 1];          for (i = n - 1; i >= 0; i--) sa[--ws[ wv[i]]] = y[i];  print("sa:",SA,n+1);//**************print**************          for ( p = 1, y[sa[0]] = 0, i = 1; i < n; i++)             y[sa[i]] = cmp(x, sa[i - 1], sa[i], j) ? p - 1 : p++;  print("array y:",y,n);     print("array x:",x,n);          swap(x, y);  printf("----------------------------------------\n");     } } #define  maxn 100int rank[maxn],height[maxn];void calheight(char*r,int*sa,int n){int i,j,k=0;for(i=1;i<=n;i++)rank[sa[i]]=i;for(i=0;i<n;height[rank[i++]]=k)for(k?k--:0,j=sa[rank[i]-1];r[i+k]==r[j+k];k++);return;}int main(){     int n=strlen(s);printf("n = %d\n",n);    int m=128;     da(s,SA,n + 1,m);     print("sa:",SA,n + 1);//**************print************** for (int i = 0; i <= n; ++i){printf("%s\n",&s[SA[i]]);}calheight(s,SA,n + 1);print("height:",height,n + 1);    return 0; } /*  * sa:    8 3 4 5 0 6 1 7 2  * */