后缀数组详解

这是一部分，以后还会更新的

 

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#define maxn 10000
using namespace std;

char s[maxn],s1[maxn];
int sa[maxn];/***sa[第几大]=第几位开始的后缀***/
int Rank[maxn];/***rank[第几位开始的后缀]=第几大***/
int t[maxn],t2[maxn],c[maxn],n;

void build_sa(int m)
{
	int i,*x=t,*y=t2;
	//基数排序
	for(i=0;i<m;++i)c[i]=0;
	for(i=0;i<n;++i)c[x[i]=s[i]]++;
	for(i=1;i<m;++i)c[i]+=c[i-1];/**计算出每种字母排第几**/
	for(i=n-1;i>=0;--i)sa[--c[x[i]]]=i;
	for(int k=1;k<=n;k<<=1)
	{
		int p=0;
		/**直接利用sa数组排序第二关键字**/
		for(i=n-k;i<n;++i)y[p++]=i;
		for(i=0;i<n;++i)
			if(sa[i]>=k)y[p++]=sa[i]-k;
		/**基数排序第一关键字**/
		for(i=0;i<m;++i)c[i]=0;
		for(i=0;i<n;++i)c[x[y[i]]]++;
		for(i=0;i<m;++i)c[i]+=c[i-1];
		for(i=n-1;i>=0;i--)sa[--c[x[y[i]]]]=y[i];
		/**根据sa和y数组计算新的x数组**/
		swap(x,y);
		p=1;x[sa[0]]=0;
		for(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;/**以后即使继续倍增，sa也不会改变，退出**/
		m=p;/**下次基数排序的最大值**/
	}
}

int m;
int cmp_suffix(char *pattern,int p)
{
	return strncmp(pattern,s+sa[p],m);
	/***************************
	用 法: int strncmp(char *str1, char *str2, int maxlen);
说明:此函数功能即比较字符串str1和str2的前maxlen个字符。
如果前maxlen字节完全相等，返回值就=0；
在前maxlen字节比较过程中，如果出现str1[n]与str2[n]不等，
则返回（str1[n]-str2[n]）。
	****************************/
}

int find(char *p)
{
	m=strlen(p);
	if(cmp_suffix(p,0)<0)return -1;/**下界为第0个后缀**/
	if(cmp_suffix(p,n-1)>0)return -1;/**上界为第n-1个后缀**/
	int l=0,r=n-1;
	while(r>=l)/**二分判断是否有该后缀与P的公共前缀=p的长度**/
	{
		int m=l+(r-l)/2;
		int res=cmp_suffix(p,m);
		if(!res)return m;
		if(res<0)r=m-1;else l=m+1;
	}
	return -1;
}


int main()
{
    scanf("%s",s);
    n=strlen(s);
    build_sa(300);/**char的最大字符值**/
    scanf("%s",s1);
    cout<<find(s1);/**从s中查找s1**/
	return 0;
}

 

 

 

下面这个是基础版，常数会大些，复杂度不变

 1 #include <iostream>
 2 #include <string>
 3 #include <cstring>
 4 #include <algorithm>
 5 #include <cstdio>
 6 #define maxlen 1000
 7 #define maxn 10000
 8 
 9 using namespace std;
10 
11 char s[maxlen],s1[maxlen];
12 int sa[maxn],tsa[maxn],Rank[maxn],trank[maxn],sum[maxn];
13 int n,m;
14 
15 void sorting(int j)
16 {
17     memset(sum,0,sizeof(sum));
18     for(int i=1;i<=n;i++)sum[Rank[i+j]]++;
19     for(int i=1;i<=maxlen;i++)sum[i]+=sum[i-1];
20     for(int i=n;i!=0;i--)tsa[sum[Rank[i+j]]--]=i;
21 
22     memset(sum,0,sizeof(sum));
23     for(int i=1;i<=n;++i)sum[Rank[i]]++;
24     for(int i=1;i<=maxlen;++i)sum[i]+=sum[i-1];
25     for(int i=n;i!=0;i--)sa[sum[Rank[tsa[i]]]--]=tsa[i];
26 }
27 
28 void build_sa()
29 {
30     int p;
31     for(int i=0;i!=n;++i)trank[i+1]=s[i];
32     for(int i=1;i<=n;++i)sum[trank[i]]++;
33     for(int i=1;i<=maxlen;++i)sum[i]+=sum[i-1];
34     for(int i=n;i!=0;i--)sa[sum[trank[i]]--]=i;
35     Rank[sa[1]]=1;
36     for(int i=2,p=1;i<=n;++i)
37     {
38         if(trank[sa[i]]!=trank[sa[i-1]])p++;
39         Rank[sa[i]]=p;
40     }
41     for(int j=1;j<=n;j*=2)
42     {
43         sorting(j);
44         trank[sa[1]]=1;p=1;
45         for(int i=2;i<=n;++i)
46         {
47             if((Rank[sa[i]]!=Rank[sa[i-1]]) || (Rank[sa[i]+j]!=Rank[sa[i-1]+j]))p++;
48             trank[sa[i]]=p;
49         }
50         for(int i=1;i<=n;++i)Rank[i]=trank[i];
51     }
52 }
53 
54 void init()
55 {
56     memset(sum,0,sizeof(sum));
57     memset(Rank,0,sizeof(Rank));
58     memset(trank,0,sizeof(trank));
59     memset(sa,0,sizeof(sa));
60     memset(tsa,0,sizeof(tsa));
61 }
62 
63 int cmp_suffix(char *pattern,int p){return strncmp(pattern,s+sa[p]-1,m);}
64 
65 int finds(char *p)
66 {
67     m=strlen(p);
68     if(cmp_suffix(p,1)<0)return -1;
69     if(cmp_suffix(p,n)>0)return -1;
70     int l=1,r=n;
71     while(r>=l)
72     {
73         int m=l+(r-l)/2;
74         int res=cmp_suffix(p,m);
75         if(!res)return m;
76         if(res<0)r=m-1;else l=m+1;
77     }
78     return -1;
79 }
80 
81 int main()
82 {
83     init();
84     scanf("%s",s);
85     n=strlen(s);
86     build_sa();
87     scanf("%s",s1);
88     cout<<finds(s1);
89     return 0;
90 }