记录编号 |
454380 |
评测结果 |
AAAAAAAAAA |
题目名称 |
[SPOJ 705] 不同的子串 |
最终得分 |
100 |
用户昵称 |
Hzoi_Mafia |
是否通过 |
通过 |
代码语言 |
C++ |
运行时间 |
0.000 s |
提交时间 |
2017-09-28 20:33:44 |
内存使用 |
0.00 MiB |
显示代码纯文本
#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
char s[50005];
int n,m;
int t1[50005],t2[50005],t3[50005],buc[50005];
int sa[50005],rank[50005],height[50005];
inline void Suffix(){
int i,j,k(0),p(0),*x(t1),*y(t2),*t;
for(i=0;i<=m;++i)buc[i]=0;
for(i=1;i<=n;++i)++buc[x[i]=s[i]];
for(i=1;i<=m;++i)buc[i]+=buc[i-1];
for(i=n;i>=1;--i)sa[buc[x[i]]--]=i;
for(j=1;p<n;j<<=1,m=p){
for(p=0,i=n-j+1;i<=n;++i)y[++p]=i;
for(i=1;i<=n;++i)
if(sa[i]>j)
y[++p]=sa[i]-j;
for(i=0;i<=m;++i)buc[i]=0;
for(i=1;i<=n;++i)t3[i]=x[y[i]];
for(i=1;i<=n;++i)++buc[t3[i]];
for(i=1;i<=m;++i)buc[i]+=buc[i-1];
for(i=n;i>=1;--i)sa[buc[t3[i]]--]=y[i];
for(t=x,x=y,y=t,x[sa[1]]=1,p=1,i=2;i<=n;++i)
x[sa[i]]=((y[sa[i]]==y[sa[i-1]])&&(y[sa[i]+j]==y[sa[i-1]+j]))?p:++p;
}
for(i=1;i<=n;++i)rank[sa[i]]=i;
for(i=1;i<=n;height[rank[i++]]=k)
for(k?--k:0,j=sa[rank[i]-1];s[i+k]==s[j+k];++k);
}
inline int gg(){
freopen("subst1.in","r",stdin);
freopen("subst1.out","w",stdout);
scanf("%s",s+1);
n=strlen(s+1);
m=130;
Suffix();
int ans(n);
for(int i=1;i<=n;++i)
ans+=n-height[i]-sa[i];
printf("%d\n",ans);
return 0;
}
int K(gg());
int main(){;}