标签:
求字符串内的非空子串的数量
后缀数组!!!
#include<algorithm> #include<cstdlib> #include<cstring> #include<cstdio> #include<cmath> using namespace std; typedef long long LL; #define N 1510 int wa[N],wb[N],ws[N],wv[N]; int rank[N],height[N],sa[N]; int r[N]; char s[N]; int ans; 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 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<<=1,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 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() { scanf("%s",s); int n=strlen(s); for (int i=0;i<n;i++) r[i]=s[i]-‘A‘+2; r[n]=0; da(r,n+1,300); calheight(r,n); for(int i=0;i<n;i++) ans+=n-i-height[rank[i]]; printf("%d",ans); return 0; }
还可以用Trie,短多了。。。
#include <cstdio> #include <cstring> #include <algorithm> #include<iostream> #include<bitset> #define mem(a,b) memset(a,b,sizeof(a)) using namespace std; char s[1505]; int sum,i,j,len,u,ch[1200000][26]; int main() { scanf("%s",s); len=strlen(s); for(i=0;i<len;i++) { u=0; for(j=i;j<len;j++) { int c=s[j]-‘A‘; if(!ch[u][c]) ch[u][c]=++sum; u=ch[u][c]; } } cout<<sum<<endl; return 0; }
标签:
原文地址:http://www.cnblogs.com/yangjiyuan/p/5376305.html