-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathSuffix_array
More file actions
35 lines (34 loc) · 993 Bytes
/
Suffix_array
File metadata and controls
35 lines (34 loc) · 993 Bytes
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
//Suffix array in O(nlog^2n)
vector<int> suffix_array(string s){
s += '$';
int n = s.size();
vector<int> p(n),c(n); // p stores suffix array and c stroes its equivalence class
{
// k = 0 phase
vector<pair<char,int>> a(n);
for(int i=0;i<n;++i) a[i] = {s[i],i};
sort(all(a));
for(int i=0;i<n;++i) p[i] = a[i].S;
c[p[0]] = 0;
for(int i=1;i<n;++i){
if(a[i].F == a[i-1].F) c[p[i]] = c[p[i-1]];
else c[p[i]] = c[p[i-1]] + 1;
}
}
int k=0;
while((1<<k) < n ){
// k -> k + 1
vector<pair<pii,int>> a(n);
for(int i=0;i<n;++i) a[i] = {{c[i],c[( i+(1<<k) )%n]} ,i};
sort(all(a));
for(int i=0;i<n;++i) p[i] = a[i].S;
c[p[0]] = 0;
for(int i=1;i<n;++i){
if(a[i].F == a[i-1].F) c[p[i]] = c[p[i-1]];
else c[p[i]] = c[p[i-1]]+1;
}
if(c[p.back()] == n-1) break;
++k;
}
return p;
}