后缀数组-倍增发(模板)
来源:互联网 发布:微波仿真软件下载 编辑:程序博客网 时间:2024/05/01 15:10
//后缀数组,倍增法实现int n,k;int rank[MAX],tmp[MAX];bool compare_sa(int i,int j){ if(rank[i]!=rank[j]) return rank[i]<rank[j]; int ri=i+k<=n?rank[i+k]:-1; int rj=j+k<=n?rank[j+k]:-1; return ri<rj;}void construct_sa(string S,int *sa){ for(int i=0;i<=n;i++) { sa[i]=i; rank[i]=i<n?S[i]:-1; } for(int k=1;k<=n;k*=2) { sort(sa,sa+n+1,compare_sa); tmp[sa[0]]=0; for(int i=1;i<=n;i++) tmp[sa[i]]=tmp[sa[i-1]]+(compare_sa(sa[i-1],sa[i])?1:0); for(int i=0;i<=n;i++) rank[i]=tmp[i]; }}/*后缀数组倍增法,基数排序实现,需要注意很多细节,还有几个地方不太理解*/char s[MAXN];int sa[MAXN],t[MAXN],t2[MAXN],c[MAXN],n;void build_sa(int m){ int i,*x=t,*y=t2; for(int i=0;i<m;i++)c[i]=0; for(int i=0;i<n;i++)c[x[i]=s[i]]++; for(int i=1;i<m;i++)c[i]+=c[i-1]; for(int i=n-1;i>=0;i--)sa[--c[x[i]]]=i; for(int k=1;k<=n;k<<=1) { 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; for(int i=0;i<m;i++)c[i]=0; for(int i=0;i<n;i++)c[x[y[i]]]++; for(int i=1;i<m;i++)c[i]+=c[i-1]; for(int i=n-1;i>=0;i--)sa[--c[x[yi]]]=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; }}
发现一个注释的很详细的网站http://www.cnblogs.com/staginner/archive/2012/02/02/2335600.html
- 后缀数组-倍增发(模板)
- 后缀数组 倍增算法模板
- 后缀数组 倍增法模板
- 后缀数组模板(倍增法)
- 后缀数组二倍增模板解析
- 后缀数组倍增算法模板详解
- 后缀数组的倍增模板与DC3模板
- SA(后缀数组)——倍增法模板
- 后缀数组(倍增)
- 后缀数组倍增法
- 后缀数组(倍增)
- 后缀数组 倍增算法
- 后缀数组,倍增算法
- 后缀数组--倍增法实现
- 倍增算法实现后缀数组
- 后缀数组之倍增法
- 后缀数组--倍增法注解
- 后缀数组之倍增算法
- Win32编程---实现点,线,面的绘制
- 递归的使用
- 索引及存储结构随想【收集中】
- 《深入理解计算机系统》第八章 异常控制流 读书笔记
- linux下touch
- 后缀数组-倍增发(模板)
- jQuery 还未生成元素时 如何绑定事件?
- 多线程模拟时钟
- 细胞
- HDU 1081 To The Max
- 统计学概要
- Insertion Sort List | leetcode
- 一些网页
- 字符串分割(C++)