求不重复最大子序列，后缀数组。
#include<stdio.h>#include<iostream>#include<algorithm>#include<string.h>using namespace std;const int nMax=20005;int r[nMax],N,ret;int sa[nMax], rank[nMax], height[nMax];int wa[nMax], wb[nMax], wv[nMax], wd[nMax];int cmp(int *r, int a, int b, int l){ return r[a] == r[b] && r[a+l] == r[b+l];}void da(int *r, int n, int m){ int i, j, p, *x = wa, *y = wb, *t; for(i = 0; i < m; i ++) wd[i] = 0; for(i = 0; i < n; i ++) wd[x[i]=r[i]] ++; for(i = 1; i < m; i ++) wd[i] += wd[i-1]; for(i = n-1; i >= 0; i --) sa[-- wd[x[i]]] = i; for(j = 1, p = 1; p < n; j *= 2, m = p) { for(p = 0, i = n-j; i < n; i ++) y[p ++] = i; for(i = 0; i < n; i ++) if(sa[i] >= j) y[p ++] = sa[i] - j; for(i = 0; i < n; i ++) wv[i] = x[y[i]]; for(i = 0; i < m; i ++) wd[i] = 0; for(i = 0; i < n; i ++) wd[wv[i]] ++; for(i = 1; i < m; i ++) wd[i] += wd[i-1]; for(i = n-1; i >= 0; i --) sa[-- wd[wv[i]]] = y[i]; for(t = x, x = y, y = t, p = 1, x[sa[0]] = 0, i = 1; i < n; i ++) { x[sa[i]] = cmp(y, sa[i-1], sa[i], j) ? p - 1: p ++; } }}void calHeight(int *r, int n){ // 求height数组。 int i, j, k = 0; for(i = 1; i <= n; i ++) rank[sa[i]] = i; // 1->n for(i = 0; i < n; i++){ for(k ? k -- : 0, j = sa[rank[i]-1]; r[i+k] == r[j+k]; k ++); height[rank[i]] = k; }}int ok(int lim){ int mi=sa[1],ma=sa[1]; for(int i=2;i<=N;i++) { if(height[i]<lim) { if(ma-mi>=lim) return 1; ma=mi=sa[i]; } else { mi=min(sa[i],mi); ma=max(sa[i],ma); } } if(ma-mi>=lim) return 1; return 0;}int main(){ // freopen("test.txt","r",stdin); while(scanf("%d",&N) && N) { for(int i=0;i<N;i++) scanf("%d",&r[i]); for(int i=0;i<=N-2;i++) r[i]=r[i+1]-r[i]+88; r[N-1]=0; da(r,N,88*2); calHeight(r,N-1); int l=0,r=N/2,mid; while(l<=r) { mid=(l+r)>>1; if(ok(mid)) { ret=mid; l=mid+1; } else r=mid-1; } ret++; if(ret<5) ret=0; printf("%d\n",ret); } return 0;}
