poj1745(后缀数组+二分)

最新推荐文章于 2018-06-14 12:27:06 发布

原创最新推荐文章于 2018-06-14 12:27:06 发布 · 402 阅读

1 ·

CC 4.0 BY-SA版权

后缀数组专栏收录该内容

6 篇文章

订阅专栏

博客介绍了如何利用后缀数组和二分搜索算法来找出字符串中最长的不重叠重复子串的最大长度。在解决过程中，需注意输入的细节，并通过将高度大于等于L的元素分组来验证条件。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目大意：求最长的不重叠重复字串的最大值

思路：二分最长的长度L，然后按照height>=L分组，可以得出每一个L都会在同一组内，而这组sa[i]max-sa[i]min>=L 则成立
这题输入的细节比较多需要注意（~~我会说我T了~~）

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fod(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=20005,inf=1e9;
int n,k,p,q=1;
int a[N],v[N],h[N],sa[2][N],rk[2][N];
void calsa(int sa[N],int rk[N],int SA[N],int RK[N])
{
    for(int i=1;i<=n;i++)v[rk[sa[i]]]=i;
    for(int i=n;i;i--)
        if(sa[i]>k)
            SA[v[rk[sa[i]-k]]--]=sa[i]-k;
    for(int i=n-k+1;i<=n;i++)SA[v[rk[i]]--]=i;
    for(int i=1;i<=n;i++)
        RK[SA[i]]=RK[SA[i-1]]+(rk[SA[i-1]]!=rk[SA[i]]||rk[SA[i]+k]!=rk[SA[i-1]+k]);
}
void getsa()
{
    memset(v,0,sizeof(v));
    p=0,q=1;
    for(int i=1;i<=n;i++)v[a[i]]++;
    for(int i=1;i<=200;i++)v[i]+=v[i-1];
    for(int i=1;i<=n;i++)
        sa[p][v[a[i]]--]=i;
    for(int i=1;i<=n;i++)
        rk[p][sa[p][i]]=rk[p][sa[p][i-1]]+(a[sa[p][i-1]]!=a[sa[p][i]]);
    for(k=1;k<n;k<<=1,swap(p,q))
        calsa(sa[p],rk[p],sa[q],rk[q]);
}
void geth()
{
    k=0;
    for(int i=1;i<=n;i++)
        if(rk[p][i]==1)h[rk[p][i]]=0;
        else 
        {
            int j=sa[p][rk[p][i]-1];
            while(a[i+k]==a[j+k])k++;
            h[rk[p][i]]=k;if(k>0)k--;
        }
}
bool check(int x)
{
    int mn=sa[p][1],mx=sa[p][1];
    for(int i=2;i<=n;i++)
        if(h[i]<x)
            mn=mx=sa[p][i];
        else 
        {
            mn=min(mn,sa[p][i]);
            mx=max(mx,sa[p][i]);
            if(mx-mn>=x)return 1;
        }
    return 0;
}
void solve()
{
    int l=0,r=n/2,ans;
    while(l<=r)
    {
        int mid=(l+r)>>1;
        if(check(mid))ans=mid,l=mid+1;
        else r=mid-1;
    }
    if(ans<4)puts("0");
    else printf("%d\n",ans+1);
}
int main()
{
    while(cin>>n&&n) {
        fo(i,1,n) cin>>a[i];
        fo(i,1,n-1) a[i]=a[i+1]-a[i]+100;
        n--;
        getsa();geth();
        solve();
    }
    return 0;
}