[Tyvj 1860] 后缀数组 倍增O(nlogn) 求sa[],height[]模板，单调栈维护，求特征值

#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]);