[POJ1743] Musical Theme [SA]

Link
POJ - http://poj.org/problem?id=1743
Luogu - https://www.luogu.org/problemnew/show/P2743

---

SA。
经典题目：
本质不同子串数量
[L,R]内本质不同子串的数量
快速求两个后缀的LCP
可重叠最长重复子串
不可重叠最长重复子串

---

首先要想到所谓的字符串算法并不一定要处理 “字符”，处理较小的数字也是没有问题的。。
做一下差分就变成求不重叠最长重复子串
显然原问题可以表示为求 $Ans=\underset{i,j}{\max }\min (abs(i-j-1),\mathrm{L}\mathrm{C}\mathrm{P}(\mathrm{i},\mathrm{j}))+1$
其中 $\mathrm{L}\mathrm{C}\mathrm{P}(\mathrm{i},\mathrm{j})$ 表示对差分后的后缀 $i$ 和 $j$ 取 $\mathrm{L}\mathrm{C}\mathrm{P}$ 的长度
（注意那个差分产生的 $-1$ 和 $+1$ ！！！不要忽略掉呀）
另外题目限定了 $Ans\ge 5$
问题有点复杂啊？不过你看到了一个 $\max$ 和一个 $\min$ 二分？
……虽然这个看起来并不是常规的最小值最大，不过套了两层好像是不太好做。
那么，二分 枚举 什么 呢
$x=\mathrm{L}\mathrm{C}\mathrm{P}(\mathrm{i},\mathrm{j})$ 如何 ？
所以现在的问题就是判断是否存在
$\mathrm{m}\mathrm{i}\mathrm{n}(\mathrm{a}\mathrm{b}\mathrm{s}(\mathrm{i}-\mathrm{j}-1),\mathrm{L}\mathrm{C}\mathrm{P}(\mathrm{i},\mathrm{j}))\ge x$
爽多了
接着把所有 $\mathrm{L}\mathrm{C}\mathrm{P}$ 满足条件的抓出来逐个判断 $abs(i-j-1)\ge x$ 是否成立即可。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<cctype>
#include<ctime>
#include<cmath>
using namespace std;
#define R register
const int INF = 0x3f3f3f3f;
const int MAXN = 4e4 + 5;
int n, S[MAXN], m, Rk[MAXN], Tp[MAXN], SA[MAXN], Tax[MAXN], Height[MAXN], Ans;
void Radix_Sort()
{
	fill(Tax, Tax + 1 + m, 0);
	for (R int i = 1; i <= n; ++i) ++Tax[Rk[i]];
	for (R int i = 1; i <= m; ++i) Tax[i] += Tax[i-1];
	for (R int i = n; i >= 1; --i) SA[Tax[Rk[Tp[i]]]--] = Tp[i];
}
void Suffix_Sort()
{
	m = 176;
	for (R int i = 1; i <= n; ++i) Rk[i] = S[i], Tp[i] = i;
	Radix_Sort();
	for (R int Len = 1, tot = 0; tot < n; Len <<= 1, m = tot)
	{
		tot = 0;
		for (R int i = n - Len + 1; i <= n; ++i) Tp[++tot] = i;
		for (R int i = 1; i <= n; ++i) if (SA[i] > Len) Tp[++tot] = SA[i] - Len;
		Radix_Sort();
		swap(Rk, Tp);
		Rk[SA[1]] = tot = 1;
		for (R int i = 2; i <= n; ++i)
		{
			if (! (Tp[SA[i]] == Tp[SA[i-1]] && Tp[SA[i]+Len] == Tp[SA[i-1]+Len]) ) ++tot;
			Rk[SA[i]] = tot;
		}
	}
}
void Get_Height()
{
	for (R int i = 1, k = 0, j; i <= n; ++i)
	{
		if (k) --k;
		j = SA[Rk[i]-1];
		while (S[i+k]==S[j+k]) ++k;
		Height[Rk[i]] = k;
	}
}
bool Check(const int& x)
{
	for (R int St = INF, Ed = -INF, i = 2; i <= n; ++i)
	{
		if (Height[i] >= x)
		{
			St = min(St, min(SA[i], SA[i-1]));
			Ed = max(Ed, max(SA[i], SA[i-1]));
			if (Ed - St - 1 >= x) return true;
		}
		else St = INF, Ed = -INF;
	}
	return false;
}
int Divide(const int& DLL, const int& DRR)
{
	if (DLL > DRR) return 0;
	const int DMM = DLL + DRR >> 1;
	if (Check(DMM)) return max(DMM, Divide(DMM + 1, DRR));
	return Divide(DLL, DMM - 1);
}
int main()
{
	while (~scanf("%d", &n))
	{
		if (!n) return 0;
		for (R int i = 1; i <= n; ++i) scanf("%d", &S[i]);
		for (R int i = n; i >= 2; --i) S[i] = S[i] - S[i-1] + 88;
		--n;
		for (R int i = 1; i <= n; ++i) S[i] = S[i+1];
		Suffix_Sort();
		Get_Height();
		Ans = Divide(1, n >> 1) + 1;
		printf("%d\n", (Ans > 4) ? Ans : 0);
	}
	return 0;
}