BZOJ 3473

Posted 2020-09-26 SiriusRen

思路：

CF原题

ZYF有题解

O(nlog^2n)

//By SiriusRen
#include <bits/stdc++.h>
using namespace std;
const int N=400050;
int n,m,q,cntA[N],cntB[N],A[N],B[N],rk[N],ht[N],sa[N],tsa[N],f[N][20];
int from[N],cnt[N],rec[N],tl[N],ans[N];
char ch[N],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+1]: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;
    }
    for(int i=1;i<=n;i++)f[i][0]=ht[i];
    for(int j=1;j<=19;j++)
        for(int i=1;i<=n;i++)
            f[i][j]=min(f[i][j-1],f[i+(1<<(j-1))][j-1]);
}
bool check(int pos,int len){
    int l=pos,r=pos;
    for(int j=19;~j;j--){
        if(l+1>=(1<<j)&&f[l+1-(1<<j)][j]>=len)l-=(1<<j);
        if(f[r+1][j]>=len)r+=(1<<j);
    }return rec[r]>=l;
}
int main(){
    scanf("%d%d",&m,&q);
    for(int i=1;i<=m;i++){
        scanf("%s",ch);
        int t=strlen(ch);
        s[n++]=‘ ‘;
        for(int j=0;j<t;j++)from[n]=i,s[n++]=ch[j];
        tl[i]=n;
    }n--,s[0]=0;SA();
    int t=1,k=0;
    for(int i=1;i<=n;i++)if(from[sa[i]]){
        if(!cnt[from[sa[i]]])k++;
        cnt[from[sa[i]]]++;
        if(k>=q){
            for(;k-(cnt[from[sa[t]]]==1)>=q;k-=(cnt[from[sa[t]]]==1),--cnt[from[sa[t++]]]);
            rec[i]=t;
        }
    }
    for(int i=1;i<=n;i++)if(from[sa[i]]){
        int l=1,r=tl[from[sa[i]]]-sa[i],dt=0;
        while(l<=r){
            int mid=(l+r)>>1;
            if(check(i,mid))dt=mid,l=mid+1;
            else r=mid-1;
        }ans[from[sa[i]]]+=dt;
    }
    for(int i=1;i<=m;i++)printf("%d ",ans[i]);
}