POJ 1743 Musical Theme(二分+后缀数组)

题意：

求不重叠的最长相同变化的子串的最长长度，比如1 2 3 4 5 6 7 8 9 10,最长长度为5，因为子串1 2 3 4 5 和 6 7 8 9 10变化都一样的

思路：

男人八题之一，然而poj上a了都一万多了。。。。

思路在论文里面讲的很清楚了（其实应该是一个区间最大sa减去最小sa大于k，而不是大于等于，dicuss有人说了这个问题，我也是这么认为的，不过因为poj数据太水导致都能过）。这里讲一个坑点，就是表示前后的差值的时候会出现负数，要手动加上88，而因为都是差值且数值都是1-88，所以m最大也不会超过N

错误及反思：

数据很水，第一次交的有个地方其实不对但是也算我ac了

代码：

#include<cstd
4000
io>
#include<cstring>
#include<algorithm>
using namespace std;
const int N = 20100;
int sa
,rak
,height
;
int arr
;
void construct(const int *s,int n,int m) {
static int t1
,t2
,c
;
int *x = t1,*y = t2;
int i,j,k,p,l;
for (i = 0; i < m; ++ i) c[i] = 0;
for (i = 0; i < n; ++ i) c[x[i] = s[i]] ++;
for (i = 1; i < m; ++ i) c[i] += c[i - 1];
for (i = n - 1; i >= 0; -- i) sa[--c[x[i]]] = i;
for (k = 1; k <= n; k <<= 1) {
p = 0;
for (i = n - k; i < n; ++ i) y[p++] = i;
for (i = 0; i < n; ++ i) if (sa[i] >= k) y[p++] = sa[i] - k;
for (i = 0; i < m; ++ i) c[i] = 0;
for (i = 0; i < n; ++ i) c[x[y[i]]] ++;
for (i = 1; i < m; ++ i) c[i] += c[i - 1];
for (i = n - 1; i >= 0; -- i) sa[--c[x[y[i]]]] = y[i];
std::swap(x,y);
p = 1; x[sa[0]] = 0;
for (i = 1; i < n; ++ i)
x[sa[i]] = y[sa[i - 1]] == y[sa[i]]
&& y[sa[i - 1] + k] == y[sa[i] + k] ? p - 1: p ++;
if (p >= n) break;
m = p;
}
for (i = 0; i < n; ++ i) rak[sa[i]] = i;
for (i = 0,l = 0; i < n; ++ i) {
if (rak[i]) {
j = sa[rak[i] - 1];
while (s[i + l] == s[j + l]) l++;
height[rak[i]] = l;
if (l) l--;
}
}
}
bool judge(int x,int n)
{
for(int i=2;i<=n;i++)
{
int minn=n,maxn=-1;
while(height[i]>=x)
{
minn=min(minn,sa[i-1]);
minn=min(minn,sa[i]);
maxn=max(maxn,sa[i-1]);
maxn=max(maxn,sa[i]);
i++;
}
if(maxn-minn>x) return true;
}
return false;
}
int main()
{
int n;
while(scanf("%d",&n)&&n)
{
for(int i=0;i<n;i++)
scanf("%d",&arr[i]);
for(int i=0;i<n-1;i++)
arr[i]=arr[i+1]-arr[i]+88;
construct(arr,n,20100);
int l=0,r=n-1;
while(l+1<r)
{
int m=(l+r)>>1;
if(judge(m,n-1)) l=m;
else r=m;
}
if(l+1<5) printf("0\n");
else printf("%d\n",l+1);
}
}