比較簡(jiǎn)單的應(yīng)用。

#include <stdio.h> #include <string.h> #define maxn 200002int wa[maxn],wb[maxn],wv[maxn],ws[maxn]; int cmp(int *r,int a,int b,int l) {return r[a]==r[b]&&r[a+l]==r[b+l];} int max(int x,int y) {return x>y?x:y;} int min(int x,int y) {return x<y?x:y;} 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; } int rank[maxn],height[maxn]; 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 r[maxn],sa[maxn]; char s[maxn]; int main() {int i,j,len;while(scanf("%s",s)!=EOF){memset(r,0,sizeof(r));len=strlen(s);s[len]=1;scanf("%s",s+len+1);int x=strlen(s);for(i=0;i<x;i++)r[i]=s[i];da(r,sa,x+1,123);calheight(r,sa,x);int ans=0;for(i=1;i<x;i++){if(height[i]>ans){if((sa[i-1]+height[i]<len&&sa[i]>len)||(sa[i-1]+height[i]>len&&sa[i]<len))ans=height[i];}}printf("%d\n",ans);} }

?

轉(zhuǎn)載于:https://www.cnblogs.com/sweat123/p/4792741.html

總結(jié)

以上是生活随笔為你收集整理的hdu1403 后缀数组的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。