首页 > 代码库 > BZOJ 3473
BZOJ 3473
思路:
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]);}
BZOJ 3473
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。