#include<iostream>
#include<stdio.h>
#include<string.h>
using namespace std;
const int maxn=200005;
char text[maxn];
int SA[maxn];
int wa[maxn],wb[maxn],wv[maxn],cnt[maxn];
int height[maxn],h[maxn];
int rank[maxn];
int R[maxn];
void DA (int *r,int *sa,int n,int m)
{
int *x=wa,*y=wb,p;
for(int i=0;i<m;i++) cnt[i]=0;
for(int i=0;i<n;i++) cnt[x[i]=r[i]]++;
for(int i=1;i<m;i++) cnt[i]+=cnt[i-1];
for(int i=n-1;i>=0;i--) sa[--cnt[x[i]]]=i;
for(int j=1,p=1;p<n;j*=2,m=p)
{
p=0;
for(int i=n-j;i<n;i++) y[p++]=i;
for(int i=0;i<n;i++) if(sa[i]>=j) y[p++]=sa[i]-j;
for(int i=0;i<n;i++) wv[i]=x[y[i]];
for(int i=0;i<m;i++) cnt[i]=0;
for(int i=0;i<n;i++) cnt[wv[i]]++;
for(int i=1;i<m;i++) cnt[i]+=cnt[i-1];
for(int i=n-1;i>=0;i--) sa[--cnt[wv[i]]]=y[i];
swap(x,y);
p=1;
x[sa[0]]=0;
for(int i=1;i<n;i++)
x[sa[i]]=((y[sa[i]]==y[sa[i-1]])&&(y[sa[i]+j]==y[sa[i-1]+j])) ? p-1:p++;
}
}
void get_H(int *r,int n)
{
for(int i=1;i<=n;i++) rank[SA[i]]=i;
int j=0,k=0;
for(int 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",text);
int nn=strlen(text);
for(int i=0;i<nn;i++)
R[i]=text[i]-'0'+1;
R[nn]=0;
DA(R,SA,nn+1,130);
get_H(R,nn);
for(int i=1;i<=nn;i++) printf("%d ",SA[i]+1);
puts("");
for(int i=1;i<=nn;i++) printf("%d ",height[i]);
return 0;
}
单调栈维护,求特征值
for (int i=1;i<=n;i++){
while(t&&height[i]<=height[q[t]])t--;
L[i]=q[t];
q[++t]=i;
}
q[0]=n+1;t=0;
for (int i=n;i;i--){
while(t&&height[i]<=height[q[t]])t--;
R[i]=q[t]-1;
q[++t]=i;
}
// for(int i=1;i<=n;i++) printf("%d %d\n ",L[i],R[i]);
for (int i=1;i<=n;i++)if(R[i]>L[i])
ans=max(ans,(ll)(R[i]-L[i]+1)*height[i]);