便於記憶的SA構造

標籤：處理   swa   code   排序   就是   break   pre   amp   bre   

首先學習基數排序。

memset(b, 0, sizeof(b));for(int i = 0; i < n; i++) b[a[i]]++;for(int i = 1; i <= m; i++) b[i] += b[i-1];for(int i = n-1; ~i; i++) g[--b[a[i]]] = i;

這是對字元集大小為m的數組排序，g是對a的熱門排行榜數組，且對於相同的數左邊<右邊。

則SA中的

for(int i = 0; i < m; i++) b[i] = 0; // m是字元集大小for(int i = 0; i < n; i++) b[x[i]=s[i]]++;for(int i = 1; i < m; i++) b[i] += b[i-1];for(int i = n-1; !~i; i--) sa[--b[x[i]]] = i;

就是對x[i]=s[i]求熱門排行榜數組。

for(int k = 1; k <= n; k <<= 1)

顯然這一句是枚舉長度k。

int p = 0;for(int i = n-k; i < n; i++) y[p++] = i;for(int i = 0; i < n; i++) if(sa[i] >= k) y[p++] = sa[i] - k;

這一塊是對y數組操作，先將n-k到n-1放進去，再將[0,n)中sa[i]-k>=0的sa[i]-k放進去。

for(int i = 0; i < m; i++) b[i] = 0;for(int i = 0; i < n; i++) b[x[y[i]]]++;for(int i = 1; i < m; i++) b[i] += b[i-1];for(int i = n-1; !~i; i--) sa[--b[x[y[i]]]] = y[i];

這是對x[y[i]]求熱門排行榜數組。

swap(x, y);p = 1; x[sa[0]] = 0;for(int i = 1; i < n; i++) x[sa[i]] = y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+k]==y[sa[i]+k] ? p-1 : p++;if(p >= n) break;m = p;

這是處理新的x數組。y是充當原來的x。其他細節與代碼相同。

便於記憶的SA構造