[USACO06DEC]Milk Patterns G

我们想要得到一个：至少出现k次的子串的最长长度。

后缀排序后可得：含有这个"出现k次的子串"的后缀，必定是连续的，且前若共有k个连续的后缀拥有这个“出现k次的子串”，那么他们连续的区间一段lcp的最小值，一定就是这个“出现k次的子串”的最大值了。

所以，求出height以后，找出：区间长度为k-1的区间的最大的最小值。

单调队列或二分求解。

由于本菜菜单调队列还是会的（单调栈实在恶心），所以就偷懒写这个码量较小的了…

#include <bits/stdc++.h>
using namespace std;
const int N=1e6+5;
int n,d,m,ans;
int s[N],sum[N],rk[N],rk2[N],tp[N],sa[N],height[N];

inline void qsort()
{
	for (register int i=0; i<=m; ++i) sum[i]=0;
	for (register int i=1; i<=n; ++i) sum[rk[i]]++;
	for (register int i=1; i<=m; ++i) sum[i]+=sum[i-1];
	for (register int i=n; i>=1; --i) sa[sum[rk[tp[i]]]]=tp[i],sum[rk[tp[i]]]--;
}

inline void SA()
{
	m=1e6;
	for (register int i=1; i<=n; ++i) rk[i]=s[i],tp[i]=i;
	qsort();
	int p=0;
	for (register int len=1; p<n; m=p,len<<=1)
	{
		p=0;
		for (register int i=n-len+1; i<=n; ++i) tp[++p]=i;
		for (register int i=1; i<=n; ++i) if (sa[i]>len) tp[++p]=sa[i]-len;
		qsort();
		memcpy(rk2,rk,sizeof(rk2));
		p=1; rk[sa[1]]=p;
		for (register int i=2; i<=n; ++i)
		{
			if (rk2[sa[i]]==rk2[sa[i-1]] && rk2[sa[i]+len]==rk2[sa[i-1]+len]) p=p; else p++;
			rk[sa[i]]=p;
		}
	}
}

inline void LCP()
{
	for (register int i=1; i<=n; ++i) rk[sa[i]]=i;
	int k=0;
	for (register int i=1; i<=n; ++i)
	{
		if (rk[i]==1) continue;
		if (k) k--;
		int j=sa[rk[i]-1];
		while (i+k<=n && j+k<=n && s[i+k]==s[j+k]) k++;
		height[rk[i]]=k;	
	}
}

int main(){
	scanf("%d%d",&n,&d);
	for (register int i=1; i<=n; ++i) scanf("%d",&s[i]);
	SA();
	LCP();
//	for (register int i=1; i<=n; ++i) printf("%d ",height[i]); puts("");
	deque<int>q;
	d--;
	for (register int i=1; i<=n; ++i)
	{
		while (q.size() && q.front()<i-d+1) q.pop_front();
		while (q.size() && height[q.back()]>=height[i]) q.pop_back();
		q.push_back(i);
		if (i>=d) ans=max(ans,height[q.front()]);
	}
	printf("%d\n",ans);
return 0;	
}