hihoCoder 1403 后缀数组 重复旋律
思路:
后缀数组 第一次写 留个模板吧
先求出后缀数组,问题转换为询问height数组中连续k-1个数的最小值的最大值,单调队列扫描一遍即可。—yousiki
手懒用得STL
//By SiriusRen
#include <deque>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define N 22222
int n,k,s[N],cntA[N],cntB[N],sa[N],tsa[N],A[N],B[N],rk[N],ht[N],ans;
deque<int>q;
int main(){
scanf("%d%d",&n,&k);
for(int i=1;i<=n;i++)scanf("%d",&s[i]);
for(int i=1;i<=n;i++)cntA[s[i]]++;
for(int i=1;i<=100;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>=1;i--)tsa[cntB[B[i]]--]=i;
for(int i=n;i>=1;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<k;i++){
while(!q.empty()&&ht[i]<ht[q.back()])q.pop_back();
q.push_back(i);
}ans=ht[q.front()];
for(int i=k;i<=n;i++){
while(!q.empty()&&ht[i]<ht[q.back()])q.pop_back();
q.push_back(i);
while(!q.empty()&&q.front()<=i-k+1)q.pop_front();
if(!q.empty())ans=max(ans,ht[q.front()]);
}printf("%d\n",ans);
}