博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
真-后缀数组模板
阅读量:5143 次
发布时间:2019-06-13

本文共 1411 字,大约阅读时间需要 4 分钟。

#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 10000000//往大的开也不会死的
using namespace std;
int wa[maxn],wb[maxn],wv[maxn],ws[maxn];
int rank[maxn],height[maxn];
int r[maxn],sa[maxn],n;
int cmp(int *r,int a,int b,int l)
{return r[a]==r[b]&&r[a+l]==r[b+l];} 
void da(int *r,int *sa,int n,int m)
{
    int i,j,p,*x=wa,*y=wb,*t;
    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;  
    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;  
        for(i=0;i<n;i++) wv[i]=x[y[i]];  //这里已经不知道写错了多少遍了
        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];
        for(t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1;i<n;i++)
        x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++; 
    }
    return;
 }
void calheight(int *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()
{
char rr[1000000];
scanf("%s",rr);
n=strlen(rr);
for(int i=0;i<n;i++)r[i]=rr[i]-'a'+1;
da(r,sa,n+1,128);//相当于传了strlen()+1过去,这里非常迷,不知道为什么要传strlen()+1过去??!
calheight(r,sa,n);//这里有一个非常坑的地方这里是n而不是n+1啊,我已经调了一个下午了QAQ,
for(int i=1;i<=n;i++)printf("%d ",sa[i]+1);
printf("\n");
for(int i=2;i<=n;i++)printf("%d ",height[i]);
}

转载于:https://www.cnblogs.com/OcahIBye/p/6706960.html

你可能感兴趣的文章
luogu1328 [NOIp2014]生活大爆炸版石头剪刀布 (模拟)
查看>>
BZOJ 1003: [ZJOI2006]物流运输trans(最短路+dp)
查看>>
求解两个字符串的最长公共子序列
查看>>
(转) 垂直同步、绘制效率、显示器刷新频率与帧率
查看>>
Typescript error
查看>>
flask记录
查看>>
学习进度条12
查看>>
spotlight_监控Linux_无需修改用户权限
查看>>
数据结构学习记录_2019.02.09
查看>>
深入分析LInux内核链表
查看>>
关于Gvim中textwidth被自动设置成78造成输入时自动换行的问题
查看>>
MATLAB绘制向量图
查看>>
10款有趣创意的LOADING等待体验动画作品
查看>>
任意的四个点,判断是不是矩形
查看>>
Java中3DES加密解密与其他语言(如C/C++)通信
查看>>
log4j分级输出日志文件
查看>>
Palindrome Number
查看>>
测试用例覆盖率converage
查看>>
extjs xtype 类型
查看>>
FileUpload控件的使用
查看>>