我倒是觉得上面那一问难一点呢。。。。
下面一个SAM就没掉了,上面还要加上倍增和马拉车
size为right集合大小,mx为最大长度,跑出parent树,OK搞定
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define N 2000010
using namespace std;
char str[N]; long long A=0;
int s[N][26],sz[N],f[N],mx[N],v[N],r[N],n,cnt=1,lst=1;
inline int extend(char c){
int p=lst,np=lst=++cnt,q,nq;
c-='a'; mx[np]=mx[p]+1; sz[np]=1;
for(;p&&!s[p][c];p=f[p]) s[p][c]=np;
if(!p) return f[np]=1;
q=s[p][c];
if(mx[q]==mx[p]+1) f[np]=q;
else{
nq=++cnt;
mx[nq]=mx[p]+1;
f[nq]=f[q]; f[q]=f[np]=nq;
memcpy(s[nq],s[q],26<<2);
for(;p&&s[p][c]==q;p=f[p]) s[p][c]=nq;
}
}
int main(){
scanf("%s",str); n=strlen(str);
for(int i=0;i<n;++i) extend(str[i]);
for(int i=1;i<=cnt;++i) ++v[mx[i]];
for(int i=1;i<=n;++i) v[i]+=v[i-1];
for(int i=cnt;i;--i) r[v[mx[i]]--]=i;
for(int p,i=cnt;i;--i){
p=r[i]; sz[f[p]]+=sz[p];
if(sz[p]>1) A=max(A,1ll*sz[p]*mx[p]);
} printf("%lld
",A);
}