POJ3261：Milk Patterns——题解
tags:
篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了POJ3261：Milk Patterns——题解相关的知识，希望对你有一定的参考价值。
http://poj.org/problem?id=3261
给一个序列，求至少出现 k 次的最长重复子串，这 k 个子串可以重叠。
论文题+傻逼题。
上一道题（POJ1743）会做即可。
还是二分长度，按照高度数组分组做，统计每个组有多少个后缀即可。
#include<algorithm> #include<iostream> #include<cstring> #include<cctype> #include<cstdio> #include<vector> #include<queue> #include<cmath> using namespace std; const int N=2e4+5; int n,m,rank[N],sa[N],height[N],w[N],s[N]; inline bool pan(int *x,int i,int j,int k){ int ti=i+k<n?x[i+k]:-1; int tj=j+k<n?x[j+k]:-1; return x[i]==x[j]&&ti==tj; } inline void SA_init(){ int *x=rank,*y=height,r=10000; for(int i=0;i<r;i++)w[i]=0; for(int i=0;i<n;i++)w[s[i]]++; for(int i=1;i<r;i++)w[i]+=w[i-1]; for(int i=n-1;i>=0;i--)sa[--w[s[i]]]=i; r=1;x[sa[0]]=0; for(int i=1;i<n;i++) x[sa[i]]=s[sa[i]]==s[sa[i-1]]?r-1:r++; for(int k=1;r<n;k<<=1){ int yn=0; for(int i=n-k;i<n;i++)y[yn++]=i; for(int i=0;i<n;i++) if(sa[i]>=k)y[yn++]=sa[i]-k; for(int i=0;i<r;i++)w[i]=0; for(int i=0;i<n;i++)++w[x[y[i]]]; for(int i=1;i<r;i++)w[i]+=w[i-1]; for(int i=n-1;i>=0;i--)sa[--w[x[y[i]]]]=y[i]; swap(x,y);r=1;x[sa[0]]=0; for(int i=1;i<n;i++) x[sa[i]]=pan(y,sa[i],sa[i-1],k)?r-1:r++; } for(int i=0;i<n;i++)rank[i]=x[i]; } inline void height_init(){ int i,j,k=0; for(i=1;i<=n;i++)rank[sa[i]]=i; for(i=0;i<n;i++){ if(k)k--; else k=0; j=sa[rank[i]-1]; while(s[i+k]==s[j+k])k++; height[rank[i]]=k; } } bool check(int k){ int cnt=1; for(int i=2;i<=n;i++){ if(height[i]>=k){ cnt++; }else cnt=1; if(cnt>=m)return 1; } return 0; } int erfen(int l,int r){ while(l<r){ int mid=(l+r+1)>>1; if(check(mid))l=mid; else r=mid-1; } return l; } int main(){ scanf("%d%d",&n,&m); for(int i=0;i<n;i++)scanf("%d",&s[i]); s[n++]=0; SA_init(); n--; height_init(); printf("%d\\n",erfen(0,n)); return 0; }
+++++++++++++++++++++++++++++++++++++++++++
+本文作者：luyouqi233。 +
+欢迎访问我的博客：http://www.cnblogs.com/luyouqi233/+
+++++++++++++++++++++++++++++++++++++++++++
以上是关于POJ3261：Milk Patterns——题解的主要内容，如果未能解决你的问题，请参考以下文章
POJ 3261 Milk Patterns（后缀数组+二分答案）