传送门：后缀数组一·重复旋律
最长不可重叠重复子串问题
二分答案，转化成判定问题
bool check(int K)
{
for(int i=1;i<=n;i++)
if(height[i]< K)
{
minsa=sa[i];
maxsa=sa[i];
}
else
{
minsa=min(minsa,sa[i]);
maxsa=max(maxsa,sa[i]);
if(maxsa-minsa>=K)return true;
}
return false;
}最长公共子串问题
求排名相邻，原来不在同一个字符串的 height 值的最大值
重复次数最多的连续字串
逃。。。
题意：最长可重叠重复K次子串问题
分析：
难~，求sa、rank、height数组跪了，好难理解啊。
suffix(i)表示从原字符串第i个字符开始到字符串结尾的后缀。我们把它所有的后缀拿出来按字典序排序
把排好序的数组记做sa。比如sa[1]=7,sa[4]=2
名次数组 Rank， Rank[i] 保存的是后缀 i 在所有后缀中从小到大排列的“名次”
height[i] 是 suffix(sa[i-1]) 和 suffix(sa[i]) 的最长公共前缀长度，即排名相邻的两个后缀的最长公共前缀长度。
有几个结论：
有了height，求最长可重叠重复K次子串就方便了。重复子串即两后缀的公共前缀，最长重复子串，等价于两后缀的最长公共前缀的最大值。问题就转化成了，求height 数组中最大的长度为 K的子序列的最小值。
后缀数组的求法有很多，最有名的是两种倍增算法和DC算法。时间复杂度上DC算法更优，但是很复杂。我们这里只介绍相对容易的倍增算法。
有空再仔细研究这个算法，先溜了。
Online AC Code
#include <iostream>
#include <string>
#include <cstring>
using namespace std;
//给定一个字符串，求至少出现K次的最长重复子串，这k个子串可以重叠。
//所谓出现K次就是在任意K个地方出现了这个子串，并不要求这K个是连续的。
int N, K;
int array[20001];
//产生后缀字符串, 并求出在全部后缀字符串字典排序后, 每个后缀字符串与其前面一个后缀字符串的最长公共子串长度
void solve()
{
int cntA[20001], cntB[20001];//表示每个数字的计数, 以解决如果有相同数字的排名问题
int sa[N + 1];//cntA的严格递增排名映射:排名->下标
int tsa[N + 1];//cntB的严格递增排名映射:排名->下标
int rank[N + 1];//以下标的为第一个字符的子串的真实排名映射: 下标->排名
int height[N + 1];//在全部后缀字符串字典排序后, 每个后缀字符串与其前面一个后缀字符串的最长公共子串长度
memset(cntA, 0, sizeof(cntA));
for(int i = 1; i <= N; ++i) ++cntA[array[i]];
for(int i = 1; i < 256; ++i) cntA[i] += cntA[i - 1]; //将小数字的计数累加到大数字的计数上
for(int i = N; i >= 1; --i) sa[cntA[array[i]]--] = i; //从N到1的循环的意思是: 如果有相同的数字, 则它们之间的排名顺序是按照它们出现的先后顺序严格增大的
//for(int i = 1; i <= N; ++i) cout << i << ": " << sa[i] << endl;
rank[sa[1]] = 1;
for(int i = 2; i <= N; ++i)
{
rank[sa[i]] = rank[sa[i - 1]]; //对于数字sa[i], 它的真正排名要么和sa[i-1]的排名一样(数字一样), 要么比sa[i-1]真的大1
if(array[sa[i]] != array[sa[i - 1]]) ++rank[sa[i]];
}
//开始基数排序
int A[N + 1], B[N + 1]; //A[i]表示array中以第i个数字为开始的子串的rank值, B[i]表示array中以第i数字后面的第i+l个数字为开始的子串的rank值
for(int l = 1; rank[sa[N]] < N; l <<= 1)//此循环只要出现ranke[sa[N]] == N, 就可以结束了
{
memset(cntA, 0, sizeof(cntA));
memset(cntB, 0, sizeof(cntB));
for(int i = 1; i <= N; ++i)
{
++cntA[A[i] = rank[i]];//算出此次循环的cntA值, 此处cntA的下标值是rank值不再是数字值
++cntB[B[i] = (i + l <= N)? rank[i + l]: 0]; //算出此次循环的cntB值, c此处的cntB的下标值是rank值不再是数字值
}
for(int i = 1; i <= N; ++i) cntB[i] += cntB[i - 1];
for(int i = N; i >= 1; --i) tsa[cntB[B[i]]--] = i;//从N到1的循环的意思是: 如果有相同的rank:B[i]值, 则它们之间的排名顺序是按照它们出现的先后顺序严格增大的
for(int i = 1; i <= N; ++i) cntA[i] += cntA[i - 1];
for(int i = N; i >= 1; --i) sa[cntA[A[tsa[i]]]--] = tsa[i];//从N到1的循环的意思是: 如果有相同的rank:A[i]值, 则它们之间的排名顺序是按照它们sa[i]的值严格增大的
rank[sa[1]] = 1;
for(int i = 2; i <= N; ++i)
{
rank[sa[i]] = rank[sa[i - 1]];
if(A[sa[i]] != A[sa[i - 1]] || B[sa[i]] != B[sa[i - 1]]) ++rank[sa[i]];
}
}
//算出height
for (int i = 1, j = 0; i <= N; ++i)
{
if(j) j--;
while (array[i + j] == array[sa[rank[i] - 1] + j]) ++j;
height[rank[i]] = j;
}
//cout << endl;
//for (int i = 1; i <= N; ++i) cout << sa[i] << ": " << height[i] << endl;
int ret = 0;
int left = 1, right = N, mid;
while(left <= right)
{
mid = (left + right) >> 1;
int cnt = 1;
for(int i = 1; i <= N; ++i)
{
if(height[i] >= mid) ++cnt;
else cnt = 1;
if(cnt >= K) break;
}
if(cnt >= K)
{
left = mid + 1;
ret = max(ret, mid);
}
else right = mid - 1;
}
cout << ret << endl;
return;
}
int main()
{
cin >> N >> K;
for(int i = 1; i <= N; ++i) cin >> array[i];
solve();
return 0;
}