题意:
- 给你一个长度为n的序列,序列中的元素大小为[1,88]。
- 现在让你找重复主题:主题指序列中的一个子串,其中满足以下条件
- 1.长度>=5
- 2.子串元素间的差值相等。
- 3.重复出现的同一主题不能有公共部分。
- 求满足题意的最大主题长度,即就是最长不可重叠子串
做法:
- 预处理的时候,把原输入的数组差分处理成新数组。求得的最大长度为k
- ans就是k+1
- 比如:n = 10
- 1 2 3 4 5 6 7 8 9 10
- 符合题意的就是 1 2 3 4 5 和 5 6 7 8 9 10
AC代码:
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#define IO ios_base::sync_with_stdio(0),cin.tie(0),cout.tie(0)
#define pb(x) push_back(x)
#define sz(x) (int)(x).size()
#define sc(x) scanf("%d",&x)
#define abs(x) ((x)<0 ? -(x) : x)
#define all(x) x.begin(),x.end()
#define mk(x,y) make_pair(x,y)
#define debug printf("!!!!!!\n")
#define fin freopen("in.txt","r",stdin)
#define fout freopen("out.txt","w",stdout)
using namespace std;
typedef long long ll;
typedef pair<int,int> PII;
const int mod = 1e9+7;
const double PI = 4*atan(1.0);
const int maxm = 1e6;
const int maxn = 1e5+10;
const int INF = 0x3f3f3f3f;
const ll LINF = 1ll<<62;
const int m = 176; //
int s[maxn]; //待排序的字符串放在s数组中,从s[0~n-1],长度为n,且最大值小于m,最后一位是0(无效值)
int sa[maxn],t1[maxn],t2[maxn],c[maxn],n; //sa[1~n]为有效值,sa[0]必定为n是无效值
int rk[maxn]; //rank[0~n-1]为有效值,rank[n]必定为0无效值
int height[maxn]; //height[1~n]
//lcp(x,y):字符串x与字符串y的最长公共前缀,在这里指x号后缀与与y号后缀的最长公共前缀
//height[i]:lcp(sa[i],sa[i-1]),即排名为i的后缀与排名为i-1的后缀的最长公共前缀
//H[i]:height[rak[i]],即i号后缀与它前一名的后缀的最长公共前缀
inline int read()
{
char x;
int u,flag = 0;
while(x = getchar(),x<'0' || x>'9') if(x == '-') flag = 1;
u = x-'0';
while(x = getchar(),x>='0' && x<='9') u = (u<<3)+(u<<1)+x-'0';
if(flag) u = -u;
return u;
}
void build_sa(int m)
{
n++;
int *x = t1,*y = t2;
for(int i=0;i<m;i++) c[i] = 0;
for(int i=0;i<n;i++) c[x[i] = s[i]]++;
for(int i=1;i<m;i++) c[i]+=c[i-1];
for(int i=n-1;i>=0;i--) sa[--c[x[i]]] = i;
for(int j=1;j<=n;j<<=1)
{
int p = 0;
for(int i = n-j;i<n;i++) y[p++] = i;
for(int i=0;i<n;i++) if(sa[i]>=j) y[p++] = sa[i]-j;
for(int i=0;i<m;i++) c[i] = 0;
for(int i=0;i<n;i++) c[x[y[i]]]++;
for(int i=1;i<m;i++) c[i]+=c[i-1];
for(int i=n-1;i>=0;i--) sa[--c[x[y[i]]]] = y[i];
swap(x,y);
p = 1;x[sa[0]] = 0;
for(int i=1;i<n;i++)
x[sa[i]] = (y[sa[i-1]] == y[sa[i]] && y[sa[i-1]+j] == y[sa[i]+j])?p-1:p++;
if(p>=n) break;
m = p;
}
n--;
int k = 0;
for(int i=0;i<=n;i++) rk[sa[i]] = i;
for(int i=0;i<n;i++){
if(k) k--;
int j = sa[rk[i]-1];
while(s[i+k] == s[j+k]) k++;
height[rk[i]] = k;
}
}
bool check(int res)
{
int minsa = INF;
int maxsa = 0;
for(int i=1;i<=n;i++)
{
if(height[i]<res){
minsa = sa[i];
maxsa = sa[i];
}
else{
maxsa = max(maxsa,sa[i]);
minsa = min(minsa,sa[i]);
if(maxsa-minsa>=res) return true;
}
}
return false;
}
int main()
{
// fin;
// fout;
while(true)
{
n = read();
if(n == 0) break;
for(int i=0;i<n;i++) s[i] = read();
// for(int i=0;i<n;i++) cout<<s[i]<<" ";
// cout<<endl;
for(int i=0;i<n-1;i++) s[i] = s[i]-s[i+1]+88; //m(max) = 88-1+88 = 175
s[n-1] = 0;
n--;
build_sa(m);
int l = 1,r = n,ans = 0;
while(l<=r)
{
int mid = (r-l)/2+l;
if(check(mid)){
ans = mid;
l = mid+1;
}else r = mid-1;
}
int res = ans+1;
if(res>=5) printf("%d\n",res);
else printf("0\n");
}
return 0;
}