HDU 6194 后缀数组

思路:

ans+=max(0,lcp(i,i+k-1)-max(ht[i],ht[i+k]));

特判k=1

//By SiriusRen
#include <bits/stdc++.h>
using namespace std;
const int N=100050;
int cntA[N],cntB[N],A[N],B[N],sa[N],tsa[N],rk[N],ht[N],f[N][20];
int n,cases,k,base[N];
char s[N];
void SA(){
    for(int i=1;i<=256;i++)cntA[i]=0;
    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;
    }base[0]=-1;
    for(int i=1;i<=n;i++)f[i][0]=ht[i],base[i]=base[i>>1]+1;
    for(int j=1;j<=18;j++)
        for(int i=1;i+(1<<(j-1))<=n;i++)
            f[i][j]=min(f[i][j-1],f[i+(1<<(j-1))][j-1]);
}
int lcp(int x,int y){
    if(x==y)return n-x+1;
    x++;
    int t=base[y-x+1];
    return min(f[x][t],f[y-(1<<t)+1][t]);
}
int main(){
    scanf("%d",&cases);
    while(cases--){
        memset(ht,0,sizeof(ht));
        long long ans=0;
        scanf("%d%s",&k,s+1),n=strlen(s+1);SA();
        if(k==1){for(int i=1;i<=n;i++)ans+=max(0,n-sa[i]-max(ht[i],ht[i+1])+1);goto ed;}
        for(int i=1;i+k-1<=n;i++)
            ans+=max(0,lcp(i,i+k-1)-max(ht[i],ht[i+k]));
        ed:printf("%lld\n",ans);
    }
}

 

posted @ 2018-08-03 15:27 SiriusRen 阅读(...) 评论(...) 编辑 收藏