http://poj.org/problem?id=1743

給一段數，求最大相似子串長度，如果沒有輸出0。

相似子串定義：

1.兩個不重疊的子串，其中一個是另一個加/減一個數得來的。

2.長度>=5

二分答案，然後想如何表示相似。

實際上我們對原數進行差分然後找相同不重疊子串即可，答案為長度+1。

但是我們考慮對於1 2 3這組數據答案為1（雖然你應該輸出0），而用上面的方法會得到2，顯然是不對的。

實際上上面差分得1 1 -3，1和1雖然為相同不重疊子串，但在原數組上他們是重疊的，所以我們規定兩個子串之間必須差一個數。

具體的方法就是處理高度函數，找每個後綴的最長公共前綴即可。

細節：

處理高度函數和二分答案的時候要把末尾的數去掉。

數字不可為負，為此可以同時加上一個數。

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cctype>
#include<cstdio>
#include<vector>
#include<queue>
#include<cmath>
using namespace std;
const int N=2e4+5;
int n,rank[N],sa[N],height[N],w[N],s[N];
inline 
 
bool pan(int *x,int i,int j,int k){
    int ti=i+k<n?x[i+k]:-1;
    int tj=j+k<n?x[j+k]:-1;
    return x[i]==x[j]&&ti==tj;
}
inline void SA_init(){
    int *x=rank,*y=height,r=200;
    for(int i=0;i<r;i++)w[i]=0;
    for(int i=0;i<n;i++)w[s[i]]++;
    for(int i=1;i<r;i++)w[i]+=w[i-1
 
];
    for(int i=n-1;i>=0;i--)sa[--w[s[i]]]=i;
    r=1;x[sa[0]]=0;
    for(int i=1;i<n;i++)
    x[sa[i]]=s[sa[i]]==s[sa[i-1]]?r-1:r++;
    for(int k=1;r<n;k<<=1){
    int yn=0;
    for(int i=n-k;i<n;i++)y[yn++]=i;
    for(int i=0;i<n;i++)
        if(sa[i]>=k)y[yn++]=sa[i]-k;
    for(int i=0;i<r;i++)w[i]=0;
    for(int i=0;i<n;i++)++w[x[y[i]]];
    for(int i=1;i<r;i++)w[i]+=w[i-1];
    for(int i=n-1;i>=0;i--)sa[--w[x[y[i]]]]=y[i];
    swap(x,y);r=1;x[sa[0]]=0;
    for(int i=1;i<n;i++)
        x[sa[i]]=pan(y,sa[i],sa[i-1],k)?r-1:r++;
    }
    for(int i=0;i<n;i++)rank[i]=x[i];
}
inline void height_init(){
    int i,j,k=0;
    for(i=1;i<=n;i++)rank[sa[i]]=i;
    for(i=0;i<n;i++){
    if(k)k--;
    else k=0;
    j=sa[rank[i]-1];
    while(s[i+k]==s[j+k])k++;
    height[rank[i]]=k;
    }
}
bool check(int k){
    int maxn,minn;
    maxn=minn=sa[1];
    for(int i=2;i<=n;i++){
    if(height[i]>=k&&i<n){
        minn=min(minn,sa[i]);
        maxn=max(maxn,sa[i]);
        continue;
    }
    if(maxn-minn>=k)return 1;
    maxn=minn=sa[i];
    }
    return 0;
}
int erfen(int l,int r){
    while(l<r){
    int mid=(l+r+1)>>1;
    if(check(mid))l=mid;
    else r=mid-1;
    }
    if(l<4)return 0;
    return l+1;
}
int main(){
    while(scanf("%d",&n)!=EOF&&n){
    for(int i=0;i<n;i++)scanf("%d",&s[i]);
    for(int i=0;i<n-1;i++)s[i]=s[i+1]-s[i]+100;
    s[n-1]=0;
    SA_init();
    n--;
    height_init();
    printf("%d\n",erfen(0,n));
    }
    return 0;
}

+++++++++++++++++++++++++++++++++++++++++++

+本文作者：luyouqi233。 　　　　　　　　　　　　　　+

+歡迎訪問我的博客：http://www.cnblogs.com/luyouqi233/+

+++++++++++++++++++++++++++++++++++++++++++

POJ1743：Musical Theme——題解