題目鏈接：點擊查看

題目大意：給出n個連續的數字組成的序列，現在要求出其中兩個不重疊的字序列，滿足兩個子序列“相似”，相似的定義是兩個子序列當且僅當長度相等并且每一位的數字差都相等，現在需要求出最大的兩個相似子序列

題目分析：因為是要求兩個不重疊的重復子序列，顯然是后綴數組的模板題目之一，不過有點不同的是我們需要先預處理轉換一下，根據題目“相似”的條件，我們很容易可以想到將整個數組轉換為差分數組，之后就可以通過二分答案，按照height數組的值將整個后綴分為不同的幾個組，在每個組中尋找答案即可，為什么要分組呢，因為height數組可以傳遞的原因，如果相鄰的height都滿足大于等于當前二分的答案時，那么這幾個height數組所代表的的sa數組，也就是后綴，他們的前綴都滿足了長度的要求，且都出現次數必定大于等于二，現在只需要對于不重疊這個因素進行判斷就好了

注意，因為求height數組的時候，我們的rk[n+1]=0，會出現rk[n+1]-1=0-1=-1的情況造成數組越界，所以我們將區間限制為0~n就好了

代碼：

#include<iostream> #include<cstdio> #include<string> #include<ctime> #include<cstring> #include<algorithm> #include<stack> #include<queue> #include<map> #include<set> #include<sstream> using namespace std;typedef long long LL;const int inf=0x3f3f3f3f;const int N=2e4+100;int sa[N]; //SA數組，表示將S的n個后綴從小到大排序后把排好序的 //的后綴的開頭位置順次放入SA中 int t1[N],t2[N],c[N];int rk[N],height[N];int s[N],n;void build_sa(int s[], int n, int m) {int i,j,p,*x=t1,*y=t2;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(j=1;j<=n;j<<=1) {p=0;for(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<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];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]+j]==y[sa[i]+j]?p-1:p++;if(p>=n) break;m=p;} }void get_height(int s[],int n) {int i,j,k=0;for(i=0;i<=n;i++)rk[sa[i]]=i;for(i=0;i<n;i++) {if(k) k--;j=sa[rk[i]-1];while(s[i+k]==s[j+k]) k++;height[rk[i]]=k;} }bool check(int x) {int mmax=-inf,mmin=inf;for(int i=2;i<=n;i++){if(height[i]>=x){mmax=max(mmax,max(sa[i],sa[i-1]));mmin=min(mmin,min(sa[i],sa[i-1]));if(mmax-mmin>x)return true;}elsemmax=-inf,mmin=inf;}return false; }int main() { // freopen("input.txt","r",stdin); // ios::sync_with_stdio(false);while(scanf("%d",&n)!=EOF&&n){int pre;scanf("%d",&pre);n--;for(int i=0;i<n;i++){int num;scanf("%d",&num);s[i]=num-pre+88;pre=num;}s[n]=0;build_sa(s, n+1, 200);getHeight(s,n);int l=1,r=inf,ans=0;while(l<=r){int mid=l+r>>1;if(check(mid)){ans=mid;l=mid+1;}elser=mid-1;}if(ans<4)printf("0\n");elseprintf("%d\n",ans+1);}return 0; }

?

總結

以上是生活随笔為你收集整理的POJ - 1743 Musical Theme(二分+后缀数组+差分数组)的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。