【*】后缀数组（dc3算法构造）

最新推荐文章于 2023-04-10 11:25:04 发布

famousDT

最新推荐文章于 2023-04-10 11:25:04 发布

阅读量1.7k

点赞数

CC 4.0 BY-SA版权

分类专栏：数据结构文章标签：算法 ini c

本文链接：https://blog.youkuaiyun.com/famousDT/article/details/7373529

数据结构专栏收录该内容

31 篇文章

订阅专栏

本文介绍了一种基于三阶符号线性时间的字符串排序算法——DC3算法。该算法通过使用特定的数据结构和优化技巧实现高效的字符串排序过程。文章详细解释了关键函数如`c0`, `c12`, `sort`, 和 `dc3` 的工作原理及其在算法中的作用。

const int Max = 200001;
int num[Max];
int r[Max * 3], sa[Max * 3];
int rank[Max], height[Max];
int wa[Max], wb[Max], wv[Max], wd[Max];
#define F(x) ((x) / 3 + ((x) % 3 == 1 ? 0 : tb))
#define G(x) ((x) < tb ? (x) * 3 + 1 : ((x) - tb) * 3 + 2)
int c0(int *r, int a, int b)
{ return r[a] == r[b] && r[a + 1] == r[b + 1] && r[a + 2] == r[b + 2]; }
int c12(int k, int *r, int a, int b)
{ if (k == 2) return r[a] < r[b] || r[a] == r[b] && c12(1, r, a + 1, b + 1);
 else return r[a]<r[b]||r[a]==r[b]&&wv[a+1]<wv[b+1];}
void sort(int *r,int *a,int *b,int n,int m)
{
     int i;
     for(i=0;i<n;i++) wv[i]=r[a[i]];
     for(i=0;i<m;i++) wd[i]=0;
     for(i=0;i<n;i++) wd[wv[i]]++;
     for(i=1;i<m;i++) wd[i]+=wd[i-1];
     for(i=n-1;i>=0;i--) b[--wd[wv[i]]]=a[i];
     return;
}
void dc3(int *r,int *sa,int n,int m)
{
     int i,j,*rn=r+n,*san=sa+n,ta=0,tb=(n+1)/3,tbc=0,p;
     r[n]=r[n+1]=0;
     for(i=0;i<n;i++) if(i%3!=0) wa[tbc++]=i;
     sort(r+2,wa,wb,tbc,m);
     sort(r+1,wb,wa,tbc,m);
     sort(r,wa,wb,tbc,m);
     for(p=1,rn[F(wb[0])]=0,i=1;i<tbc;i++)
     rn[F(wb[i])]=c0(r,wb[i-1],wb[i])?p-1:p++;
     if(p<tbc) dc3(rn,san,tbc,p);
     else for(i=0;i<tbc;i++) san[rn[i]]=i;
     for(i=0;i<tbc;i++) if(san[i]<tb) wb[ta++]=san[i]*3;
     if(n%3==1) wb[ta++]=n-1;
     sort(r,wb,wa,ta,m);
     for(i=0;i<tbc;i++) wv[wb[i]=G(san[i])]=i;
     for(i=0,j=0,p=0;i<ta && j<tbc;p++)
     sa[p]=c12(wb[j]%3,r,wa[i],wb[j])?wa[i++]:wb[j++];
     for(;i<ta;p++) sa[p]=wa[i++];
     for(;j<tbc;p++) sa[p]=wb[j++];
     return;
}