題目鏈接
有N(1<=N<=20000)個音符的序列來表示一首樂曲,每個音符都是1..88範圍內的整數,現在要找一個重複的子串,它需要滿足如下條件:1.長度至少爲5個音符。 2.在樂曲中重複出現(就是出現過至少兩次)。(可能經過轉調,“轉調”的意思是主題序列中每個音符都被加上或減去了同一個整數值) 3.重複出現的同一主題不能有公共部分。
既然可以同時減去或者加上一個數,那麼其實就是換作算差值,使得差值是一個相同的子串就是可以的了。
於是,問題變成了求最長的兩個相同且不重疊的子串,這裏的不重疊是保證兩者的距離大於等於len,並且相同要指的是兩者相同部分大於等於len-1即可,因爲算的是差值,所以無法避免的是第一個就被略過了。
所以,我們可以去二分答案這個長度len,然後因爲height的傳遞性,我們可以來判斷height小於len-1之後就看前面的區間的最大的sa和最小的sa是否差值大於等於了len。
#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <limits>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <bitset>
//#include <unordered_map>
//#include <unordered_set>
#define lowbit(x) ( x&(-x) )
#define pi 3.141592653589793
#define e 2.718281828459045
#define INF 0x3f3f3f3f
#define HalF (l + r)>>1
#define lsn rt<<1
#define rsn rt<<1|1
#define Lson lsn, l, mid
#define Rson rsn, mid+1, r
#define QL Lson, ql, qr
#define QR Rson, ql, qr
#define myself rt, l, r
using namespace std;
typedef unsigned long long ull;
typedef unsigned int uit;
typedef long long ll;
const int maxN = 2e4 + 7;
struct SA
{
int n, m;
int s[maxN];
int y[maxN], x[maxN], c[maxN], sa[maxN], rk[maxN], height[maxN];
inline void get_SA()
{
for(int i=1; i<=m; i++) c[i] = 0; //桶的初始化
for(int i=1; i<=n; i++) ++c[x[i] = s[i]];
for(int i=2; i<=m; i++) c[i] += c[i - 1]; //利用差分前綴和的思想知道每個關鍵字最多是在第幾名
for(int i=n; i>=1; i--) sa[c[x[i]]--] = i;
for(int k=1; k<=n; k<<=1)
{
int num = 0;
for(int i=n - k + 1; i<=n; i++) y[++num] = i;
for(int i=1; i<=n; i++) if(sa[i] > k) y[++num] = sa[i] - k; //是否可以作爲第二關鍵字
for(int i=1; i<=m; i++) c[i] = 0;
for(int i=1; i<=n; i++) c[x[i]]++; //因爲上一次循環已經求出這次的第一關鍵字了
for(int i=2; i<=m; i++) c[i] += c[i - 1];
for(int i=n; i>=1; i--) //在同一第一關鍵字下,按第二關鍵字來排
{
sa[c[x[y[i]]]--] = y[i];
y[i] = 0;
}
swap(x, y);
x[sa[1]] = 1; num = 1;
for(int i=2; i<=n; i++)
{
x[sa[i]] = (y[sa[i]] == y[sa[i - 1]] && y[sa[i] + k] == y[sa[i - 1] + k]) ? num : ++num;
}
if(num == n) break;
m = num;
}
}
inline void get_height()
{
int k = 0;
for(int i=1; i<=n; i++) rk[sa[i]] = i;
for(int i=1; i<=n; i++)
{
if(rk[i] == 1) continue; //第一名的height爲0
if(k) k--; //height[i] >= height[i - 1] - 1
int j = sa[rk[i] - 1];
while(j + k <= n && i + k <= n && s[i + k] == s[j + k]) k++;
height[rk[i]] = k;
}
}
inline void clear()
{
n = 0; m = 200;
}
} sa;
int N, a[maxN];
bool check(int len)
{
int mx, mn;
mx = mn = sa.sa[1];
for(int i=2; i<=N; i++)
{
if(sa.height[i] < len - 1)
{
if(mx - mn >= len) return true;
mx = mn = sa.sa[i];
}
else
{
mx = max(mx, sa.sa[i]);
mn = min(mn, sa.sa[i]);
}
}
return mx - mn >= len;
}
signed main()
{
while(scanf("%d", &N) && N)
{
sa.clear();
for(int i=1; i<=N; i++) scanf("%d", &a[i]);
for(int i=1; i<=N; i++) sa.s[++sa.n] = a[i] - a[i - 1] + 100;
sa.get_SA();
sa.get_height();
// for(int i=1; i<=N; i++) printf("%3d%c", sa.s[i], i == N ? '\n' : ' ');
// for(int i=1; i<=N; i++) printf("%3d%c", sa.sa[i], i == N ? '\n' : ' ');
// for(int i=1; i<=N; i++) printf("%3d%c", sa.height[i], i == N ? '\n' : ' ');
int l = 5, r = N / 2, mid, ans = 0;
while(l <= r)
{
mid = HalF;
if(check(mid))
{
ans = mid;
l = mid + 1;
}
else r = mid - 1;
}
printf("%d\n", ans);
}
return 0;
}