BZOJ 3238 后缀数组+单调栈

来源：互联网发布：免费交换机网管软件编辑：程序博客网时间：2024/06/16 05:30

单调栈跑两遍求出来 ht[i]为最小值的那段区间

//By SiriusRen#include <cstdio>#include <cstring>#include <algorithm>using namespace std;const int N=500050;int n,cntA[N],cntB[N],A[N],B[N],sa[N],tsa[N],rk[N],ht[N],stka[N],stkb[N],top,l[N],r[N];char s[N];void SA(){    for(int i=1;i<=n;i++)cntA[s[i]]++;    for(int i=1;i<=256;i++)cntA[i]+=cntA[i-1];    for(int i=n;i;i--)sa[cntA[s[i]]--]=i;    rk[sa[1]]=1;    for(int i=2;i<=n;i++)rk[sa[i]]=rk[sa[i-1]]+(s[sa[i]]!=s[sa[i-1]]);    for(int l=1;rk[sa[n]]<n;l<<=1){        memset(cntA,0,sizeof(cntA));        memset(cntB,0,sizeof(cntB));        for(int i=1;i<=n;i++)            cntA[A[i]=rk[i]]++,            cntB[B[i]=(i+l<=n?rk[i+l]:0)]++;        for(int i=1;i<=n;i++)cntA[i]+=cntA[i-1],cntB[i]+=cntB[i-1];        for(int i=n;i;i--)tsa[cntB[B[i]]--]=i;        for(int i=n;i;i--)sa[cntA[A[tsa[i]]]--]=tsa[i];        rk[sa[1]]=1;        for(int i=2;i<=n;i++)            rk[sa[i]]=rk[sa[i-1]]+(A[sa[i]]!=A[sa[i-1]]||B[sa[i]]!=B[sa[i-1]]);    }    for(int i=1,j=0;i<=n;i++){        j=j?j-1:0;        while(s[i+j]==s[sa[rk[i]-1]+j])j++;        ht[rk[i]]=j;    }}int main(){    scanf("%s",s+1);    n=strlen(s+1),SA();    for(int i=1;i<=n;i++){        while(top&&stka[top]>ht[i])top--;        stka[++top]=ht[i],stkb[top]=i;        l[i]=stkb[top-1]+1;    }top=0;stkb[0]=n+1;    for(int i=n;i;i--){        while(top&&stka[top]>=ht[i])top--;        stka[++top]=ht[i],stkb[top]=i;        r[i]=stkb[top-1]-1;    }    long long ans=(long long)n*(n-1)*(n+1)>>1;    for(int i=2;i<=n;i++)        ans-=(long long)(i-l[i]+1)*(r[i]-i+1)*2*ht[i];    printf("%lld\n",ans);}

0 0